1 2 34

46024阅读
36回复

我matlab学习之路

离线守夜人

UID ：118005

注册：2015-09-16
登录：2025-04-27
发帖：92
等级：论坛版主

30楼发表于: 2018-03-17 08:23:45

回澳。。的帖子

澳。。:膜拜大佬[表情] (2018-03-16 16:32) `K^j:fE7n

二大爷

离线守夜人

UID ：118005

注册：2015-09-16
登录：2025-04-27
发帖：92
等级：论坛版主

31楼发表于: 2018-03-17 08:23:59

回澳。。的帖子

澳。。:膜拜大佬[表情] (2018-03-16 16:32) q1Ja*=r

是你吗？

离线yecat20

高频微波线路板生产13823677869叶少

UID ：126335

注册：2018-02-27
登录：2025-09-27
发帖：3313
等级：论坛版主

32楼发表于: 2018-03-17 13:08:51

感谢分享

专注高频线路板pcb生产（板材有：Rogers5880 4350B.4003, TACONIC ，ARLON, F4B,TP-2等），13823677869，QQ：116425136

离线守夜人

UID ：118005

注册：2015-09-16
登录：2025-04-27
发帖：92
等级：论坛版主

33楼发表于: 2019-01-15 18:59:18

根据仿真试验结果绘制出听者接收到信号频率变化曲线以及绘制时间频率图，并生成相应的声音，是二维多普勒处理。

%根据仿真试验结果绘制出听者接收到信号频率变化曲线以及绘制时间频率图，并生成相应的声音
% z=linspace(0,8);
x0=150; %m
v=50; %车速m/s
y0=30; %m
c=330; %声音在空气中的速度m/s
f0=2000; %声音的頻率Hz
fs=8000; %采祥頻率
t=0:1/fs:6; %规定t的范围和步进
x_t=-x0+v.*t;%火车与观察者的水平距禽
y=y0; %火车与观察者的垂直距窩
r=sqrt(x_t.^2+y.^2); %火车与观察者同距高
costheta=x_t./r; %声速矢量与火车和听者同達銭来角的余弦
f=f0./(1+v.*costheta./c); %听者接收到的頻率
fmin=min(f);
fmax=max(f);
voltage=(f-fmin)./(fmax- fmin)*2-1; %日一-化調頻申圧在-1 ~ +1之同
signal=0.5*vco(voltage,[fmin fmax],fs); %圧控震蕩器，由輸入申圧控制輸出信号頻率
u=sin(2*pi*f0*t); %声源友出的信号
sound(u,fs);%播放声源发出的信号
pause(5);%暂停
sound(signal,fs);%播放听者接收到的信号
%audiowrite(signal,fs,'dopp.wav');%音频数据存盘 %老版本可以用wavwrite,新版本建议用audiowrite
figure(1);
plot(t,f);
xlabel('Time/s');
ylabel('Freq/Hz');%做出听者接收到的信号的频率变化曲线
figure(2);
spectrogram(signal,kaiser(256,5),220,512,fs,'yaxis');%做出时间频率图
axis([0 6,1000,3000])

离线守夜人

UID ：118005

注册：2015-09-16
登录：2025-04-27
发帖：92
等级：论坛版主

34楼发表于: 2019-01-15 19:01:00

基于上次二维的多普勒处理，做三维的，算法的思路是将物体在空间移动的路线看做空间曲线，然后将曲线对三个面投影，对三个投影做二维的多普勒处理，再合成。

%基本参数
x0=150; %m
y0=30; %m 起始位置
z0=30; %m 起始位置
fs=8000; %采祥頻率
t=0:1/fs:50;%规定t的范围和步进
v=30; %车速m/s
c=330; %声音在空气中的速度m/s
f0=2000; %声音的頻率Hz
%飞机移动路线，等同输入信号
figure(1);
x=t.*sin(t)+x0;
y=t.*cos(t)+y0;
z=t+z0; %飞机路线公式
%做多普勒处理
%对xoy平面投影的处理
xoy_x=x;%xoy中x轴
xoy_y=y;%xoy中y轴
r=sqrt(xoy_x.^2+xoy_y.^2); %飞机与观察者同距高
costheta=xoy_x./r; %声速矢量与火车和听者同達銭来角的余弦
fx=f0./(1+v.*costheta./c); %听者接收到的頻率
%对yoz平面投影的处理
yoz_y=y;%yoz中y轴
yoz_z=z;%yoz中z轴
r=sqrt(yoz_y.^2+yoz_z.^2); %飞机与观察者同距高
costheta=yoz_y./r; %声速矢量与火车和听者同達銭来角的余弦
fy=f0./(1+v.*costheta./c); %听者接收到的頻率
%对zox轴投影的处理
zox_z=z;%zox中z轴
zox_x=x;%zox中x轴
r=sqrt(zox_z.^2+zox_x.^2); %飞机与观察者同距高
costheta=zox_z./r; %声速矢量与火车和听者同達銭来角的余弦
fz=f0./(1+v.*costheta./c); %听者接收到的頻率
%生成图形
%飞机移动路线
plot3(x,y,z);
xlabel('横向距离');
ylabel('纵向距离');
zlabel('高度');
grid on;
title('飞机移动路线');%给坐标轴加注释
%多普勒处理
figure(2);
plot3(fx,fy,fz);
xlabel('X/T');
ylabel('Y/T');
zlabel('Z/T');
grid on;
title('飞机移动路线的多普勒处理');%给坐标轴加注释
%三维移动路线对平面的投影
figure(3);
subplot(3,1,1);
plot(x,y);
grid on;
title('飞机移动路线对平面的投影');%给坐标轴加注释
xlabel('XOY投影');
subplot(3,1,2);
plot(y,z);
xlabel('YOZ投影');
grid on;
subplot(3,1,3);
plot(z,x);
xlabel('ZOX投影');
grid on;
%投影的多普勒处理
figure(4)
subplot(3,1,1);
plot(fx,fy);
title('投影的多普勒处理');%给坐标轴加注释
xlabel('XOY投影');
subplot(3,1,2);
plot(fy,fz);
xlabel('YOZ投影');
subplot(3,1,3);
plot(fz,fx);
xlabel('ZOX投影');

离线且惜

学而时习之

UID ：134512

注册：2019-11-26
登录：2021-05-19
发帖：31
等级：仿真新人

35楼发表于: 2020-04-18 22:35:34

优秀的人啊坚持下去一定会成功的

加油

离线caijun1234

UID ：115829

注册：2015-04-23
登录：2025-08-20
发帖：541
等级：积极交流四级

36楼发表于: 2022-03-08 22:40:16

发帖回复