登 录 註 冊
  • 1205阅读
  • 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 .. jK6dI 7h  
jK6dI 7h  
亲爱的朋友您仅能浏览部分内容,查看全部内容及附件请先 登录注册

离线守夜人

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

离线yecat20

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

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

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

离线守夜人

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

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