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

离线守夜人

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

离线yecat20

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

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

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

离线守夜人

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

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