帮我看看matlab程序!%干涉信号构造t=0:0.0000001:0.0040;dt=sin(2*300*pi*t);d1=0.4;d2=0.4+dt;%微小振动d(t)=5e-3*sin(2*10*pi*t)d3=0.2;d4=0.2;bc=5e-3;e0=50*cos(2*pi*8e5*t);subplot(411)plot(t,e0);xlabel('t');ylabel('原始光信号E0');subpl
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/29 02:56:54
![帮我看看matlab程序!%干涉信号构造t=0:0.0000001:0.0040;dt=sin(2*300*pi*t);d1=0.4;d2=0.4+dt;%微小振动d(t)=5e-3*sin(2*10*pi*t)d3=0.2;d4=0.2;bc=5e-3;e0=50*cos(2*pi*8e5*t);subplot(411)plot(t,e0);xlabel('t');ylabel('原始光信号E0');subpl](/uploads/image/z/12540967-7-7.jpg?t=%E5%B8%AE%E6%88%91%E7%9C%8B%E7%9C%8Bmatlab%E7%A8%8B%E5%BA%8F%21%25%E5%B9%B2%E6%B6%89%E4%BF%A1%E5%8F%B7%E6%9E%84%E9%80%A0t%3D0%3A0.0000001%3A0.0040%3Bdt%3Dsin%282%2A300%2Api%2At%29%3Bd1%3D0.4%3Bd2%3D0.4%2Bdt%3B%25%E5%BE%AE%E5%B0%8F%E6%8C%AF%E5%8A%A8d%28t%29%3D5e-3%2Asin%282%2A10%2Api%2At%29d3%3D0.2%3Bd4%3D0.2%3Bbc%3D5e-3%3Be0%3D50%2Acos%282%2Api%2A8e5%2At%29%3Bsubplot%28411%29plot%28t%2Ce0%29%3Bxlabel%28%27t%27%29%3Bylabel%28%27%E5%8E%9F%E5%A7%8B%E5%85%89%E4%BF%A1%E5%8F%B7E0%27%29%3Bsubpl)
帮我看看matlab程序!%干涉信号构造t=0:0.0000001:0.0040;dt=sin(2*300*pi*t);d1=0.4;d2=0.4+dt;%微小振动d(t)=5e-3*sin(2*10*pi*t)d3=0.2;d4=0.2;bc=5e-3;e0=50*cos(2*pi*8e5*t);subplot(411)plot(t,e0);xlabel('t');ylabel('原始光信号E0');subpl
帮我看看matlab程序!
%干涉信号构造
t=0:0.0000001:0.0040;
dt=sin(2*300*pi*t);
d1=0.4;
d2=0.4+dt;%微小振动d(t)=5e-3*sin(2*10*pi*t)
d3=0.2;
d4=0.2;
bc=5e-3;
e0=50*cos(2*pi*8e5*t);
subplot(411)
plot(t,e0);
xlabel('t');ylabel('原始光信号E0');
subplot(412);
plot(t,dt);
xlabel('t');ylabel('振动信号')
em=50*cos(2*pi*8e5*t+2*pi*(d1+2*d3+d4)/bc);
ep=50*cos(2*pi*8e5*t+2*pi*(d1+2*d2+d4)/bc);
ef=em+ep;
subplot(413)
plot(t,ef);
xlabel('t');ylabel('干涉信号');
%低通滤波器
wp=300;ws=350;fs=1000;rp=0.5;rs=20;
[n,wn]=buttord(wp/(fs/2),ws/(fs/2),rp,rs);
[b,a]=butter(n,wn);
[h,w]=freqz(b,a);
y=filtfilt(b,a,ef);
subplot(414)
plot(t,y);
xlabel('t');ylabel('滤波所得振动信号');
我低通滤波的波形不对,怎么改呀(低通滤波后得的波形要和第二个波形差不多),【把它的这个程序过程能不能帮我解释下】
帮我看看matlab程序!%干涉信号构造t=0:0.0000001:0.0040;dt=sin(2*300*pi*t);d1=0.4;d2=0.4+dt;%微小振动d(t)=5e-3*sin(2*10*pi*t)d3=0.2;d4=0.2;bc=5e-3;e0=50*cos(2*pi*8e5*t);subplot(411)plot(t,e0);xlabel('t');ylabel('原始光信号E0');subpl
你的原始信号有问题
t=0:99;
e0=0.3*sin(0.8*pi*t);
subplot(411)
plot(t,e0);
xlabel('t');ylabel('原始光信号E0');
subplot(412);
dt=0.6*sin(0.1*pi*t);
plot(t,dt);
xlabel('t');ylabel('振动信号')
subplot(412)
ef=e0+dt;
subplot(413)
plot(t,ef);
xlabel('t');ylabel('干涉信号');
%低通滤波器
wp=300;ws=350;fs=1000;rp=0.5;rs=20;
[n,wn]=buttord(wp/(fs/2),ws/(fs/2),rp,rs);
[b,a]=butter(n,wn);
[h,w]=freqz(b,a);
y=filtfilt(b,a,ef);
subplot(414)
plot(t,y);
xlabel('t');ylabel('滤波所得振动信号');