登 录 註 冊
    
  • 2174阅读
  • 34回复

我matlab学习之路 [复制链接]

上一主题 下一主题
离线守夜人

只看该作者 30楼 发表于: 2018-03-17
回 澳。。 的帖子
[quote]澳。。:膜拜大佬  (2018-03-16 16:32) [url=http://bbs.rfeda.cn/job.php?action=topost&tid=110052&pid=1874643][/url .. H2iC? cSR  
H2iC? cSR  
亲爱的朋友您仅能浏览部分内容,查看全部内容及附件请先 登录注册

离线守夜人

只看该作者 31楼 发表于: 2018-03-17
回 澳。。 的帖子
[quote]澳。。:膜拜大佬  (2018-03-16 16:32) < .. H2iC? cSR  
H2iC? cSR  
亲爱的朋友您仅能浏览部分内容,查看全部内容及附件请先 登录注册

离线yecat20

只看该作者 32楼 发表于: 2018-03-17
感谢分享[s:9 .. H2iC? cSR  
H2iC? cSR  
亲爱的朋友您仅能浏览部分内容,查看全部内容及附件请先 登录注册

专注高频线路板pcb生产(板材有:Rogers5880 4350B.4003, TACONIC ,ARLON, F4B,TP-2等),13823677869,QQ:116425136
离线守夜人

只看该作者 33楼 发表于: 01-15
根据仿真试验结果绘制出听者接收到信号频率变化曲线以及绘制时间频率图,并生成相应的声音,是二维多普勒处理。[code]%根据仿真试验结果绘制出听者接收到信号频率变化曲线以及绘制时间频率图,并生成相应的声音 H2iC? cSR  
% z=linspace(0,8); H2iC? cSR  
x0=150;  %m H2iC? cSR  
v=50;    %车速m/s H2iC? cSR  
y0=30;   %m H2iC? cSR  
c=330;   %声音在空气中的速度m/s H2iC? cSR  
f0=2000;  %声音的頻率Hz H2iC? cSR  
fs=8000;  %采祥頻率 H2iC? cSR  
t=0:1/fs:6;  %规定t的范围和步进 H2iC? cSR  
x_t=-x0+v.*t;%火车与观察者的水平距禽 H2iC? cSR  
y=y0;     %火车与观察者的垂直距窩 H2iC? cSR  
r=sqrt(x_t.^2+y.^2);    %火车与观察者同距高 H2iC? cSR  
costheta=x_t./r;  %声速矢量与火车和听者同達銭来角的余弦 H2iC? cSR  
f=f0./(1+v.*costheta./c);   %听者接收到的頻率 H2iC? cSR  
fmin=min(f); H2iC? cSR  
fmax=max(f); H2iC? cSR  
voltage=(f-fmin)./(fmax- fmin)*2-1; %日一-化調頻申圧在-1 ~ +1之同 H2iC? cSR  
signal=0.5*vco(voltage,[fmin fmax],fs);   %圧控震蕩器,由輸入申圧控制輸出信号頻率 H2iC? cSR  
u=sin(2*pi*f0*t);   %声源友出的信号 H2iC? cSR  
sound(u,fs);%播放声源发出的信号 H2iC? cSR  
pause(5);%暂停 H2iC? cSR  
sound(signal,fs);%播放听者接收到的信号 H2iC? cSR  
%audiowrite(signal,fs,'dopp.wav');%音频数据存盘 %老版本可以用wavwrite,新版本建议用audiowrite H2iC? cSR  
figure(1); H2iC? cSR  
plot(t,f); H2iC? cSR  
xlabel('Time/s'); H2iC? cSR  
ylabel('Freq/Hz');%做出听者接收到的信号的频率变化曲 .. H2iC? cSR  
H2iC? cSR  
亲爱的朋友您仅能浏览部分内容,查看全部内容及附件请先 登录注册

离线守夜人

只看该作者 34楼 发表于: 01-15
基于上次二维的多普勒处理,做三维的,算法的思路是将物体在空间移动的路线看做空间曲线,然后将曲线对三个面投影,对三个投影做二维的多普勒处理,再合成。[code]%基本参数 `z`=!1  
x0=150;  %m `z`=!1  
y0=30;   %m 起始位置 `z`=!1  
z0=30;   %m 起始位置 `z`=!1  
fs=8000;  %采祥頻率 `z`=!1  
t=0:1/fs:50;%规定t的范围和步进 `z`=!1  
v=30;    %车速m/s `z`=!1  
c=330;   %声音在空气中的速度m/s `z`=!1  
f0=2000;  %声音的頻率Hz `z`=!1  
`z`=!1  
`z`=!1  
%飞机移动路线,等同输入信号 `z`=!1  
figure(1); `z`=!1  
x=t.*sin(t)+x0; `z`=!1  
y=t.*cos(t)+y0; `z`=!1  
z=t+z0;     %飞机路线公式 `z`=!1  
`z`=!1  
%做多普勒处理 `z`=!1  
`z`=!1  
`z`=!1  
%对xoy平面投影的处理 `z`=!1  
xoy_x=x;%xoy中x轴 `z`=!1  
xoy_y=y;%xoy中y轴 `z`=!1  
r=sqrt(xoy_x.^2+xoy_y.^2);    %飞机与观察者同距高 `z`=!1  
costheta=xoy_x./r;  %声速矢量与火车和听者同達銭来角的余弦 `z`=!1  
fx=f0./(1+v.*costheta./c);   %听者接收到的頻率 `z`=!1  
`z`=!1  
%对yoz平面投影的处理 `z`=!1  
yoz_y=y;%yoz中y轴 `z`=!1  
yoz_z=z;%yoz中z轴 `z`=!1  
r=sqrt(yoz_y.^2+yoz_z.^2);    %飞机与观察者同距高 `z`=!1  
costheta=yoz_y./r;  %声速矢量与火车和听者同達銭来角的余弦 `z`=!1  
fy=f0./(1+v.*costheta./c);   %听者接收到的頻率 `z`=!1  
`z`=!1  
%对zox轴投影的处理 `z`=!1  
zox_z=z;%zox中z轴 `z`=!1  
zox_x=x;%zox中x轴 `z`=!1  
r=sqrt(zox_z.^2+zox_x.^2);    %飞机与观察者同距高 `z`=!1  
costheta=zox_z./r;  %声速矢量与火车和听者同達銭来角的余弦 `z`=!1  
fz=f0./(1+v.*costheta./c);   %听者接收到的頻率 `z`=!1  
`z`=!1  
%生成图形 `z`=!1  
%飞机移动路线 `z`=!1  
plot3(x,y,z); `z`=!1  
xlabel('横向距离'); `z`=!1  
ylabel('纵向距离'); `z`=!1  
zlabel('高度'); `z`=!1  
grid on; `z`=!1  
title('飞机移动路线');%给坐标轴加注释 `z`=!1  
`z`=!1  
%多普勒处理 `z`=!1  
figure(2); `z`=!1  
plot3(fx,fy,fz); `z`=!1  
xlabel('X/T'); `z`=!1  
ylabel('Y/T'); `z`=!1  
zlabel('Z/T'); `z`=!1  
grid on; `z`=!1  
title('飞机移动路线的多普勒处理');%给坐标轴加注释 `z`=!1  
`z`=!1  
%三维移动路线对平面的投影 `z`=!1  
`z`=!1  
figure(3); `z`=!1  
subplot(3,1,1); `z`=!1  
plot(x,y); `z`=!1  
grid on; `z`=!1  
title('飞机移动路线对平面的投影');%给坐标轴加注释 `z`=!1  
xlabel('XOY投影'); `z`=!1  
subplot(3,1,2); `z`=!1  
plot(y,z); `z`=!1  
xlabel('YOZ投影'); `z`=!1  
grid on; `z`=!1  
subplot(3,1,3); `z`=!1  
plot(z,x); `z`=!1  
xlabel('ZOX投影'); `z`=!1  
grid on; `z`=!1  
`z`=!1  
%投影的多普勒处理 `z`=!1  
figure(4) `z`=!1  
subplot(3,1,1); `z`=!1  
plot(fx,fy); `z`=!1  
title('投 .. `z`=!1  
`z`=!1  
亲爱的朋友您仅能浏览部分内容,查看全部内容及附件请先 登录注册

快速回复
限1000 字节
 
上一个 下一个