具体 *** 作步骤如下:
1、首先,打开Matlab软件,单击[打开]或创建一个新的调试文件,如下图所示,然后进入下一步。
2、其次,单击代码前面的横杠设置断点,如下图所示,然后进入下一步。
3、接着,单击[run],如下图所示,然后进入下一步。
4、然后,光标停在该版本的断点处,然后按[F10]逐步执行,如下图所示,然后进入下一步。
5、最后,将光标放在要观察的变量上,停留一段时间,以查看变量值,再按F5执行完成,如下图所示。这样,问题就解决了。
1、打开Matlab软件,【打开】或新建调试文件。
2、点击代码前的横杠【设置断点】。
3、点击【运行】。
4、光标停到断点处,按【F10】单步执行。
5、将鼠标放在要观察的变量上【停留片刻】,查看变量值(当矩阵太大时,只显示矩阵的维数)。最后F5执行完成。
%2FSK仿真程序clear all
fs=2000
dt=1/fs
f1=20
f2=120
a=round(rand(1,10))
g1=a
g2= ~a
g11=(ones(1,2000))'*g1
g1a=g11(:)'
g21=(ones(1,2000))'*g2
g2a=g21(:)'
t=0:dt:10-dt
t1=length(t)
fsk1 =g1a.*cos(2*pi*f1.*t)%%%%%%%%%%%%%%%%%%%%%%%%%
fsk2 =g2a.*cos(2*pi*f2.*t)%%%%%%%%%%%%%%%%%%%%%%%%%
fsk=fsk1+fsk2 %产生的信号
no=0.0*randn
sn=fsk+no
subplot(311)plot(t,no)
subplot(312)plot(t,fsk)
% title(‘信号’)ylabel(‘’信号幅度)
subplot(313)plot(t,sn)
% title(‘叠加噪声后的信号’);ylabel(‘幅度A’)xlabel(‘时间t’)
figure(2)
b1=fir1(101,[10/80020/800])
b2=fir1(101,[90/800110/800])
H1=filter(b1,1,sn)H2=filter(b2,1,sn)
subplot(211)plot(t,H1)
% title(‘经过带通滤波器后f1的波形’)ylabel(‘幅度’)
subplot(212)plot(t,H2)
% title(‘经过带通滤波器后f2的波形’)ylabel(‘幅度’)xlabel(‘t’)
sw1=H1.*H1
sw2=H2.*H2
figure(3)
subplot(211)plot(t,sw1)
% title(‘经过相乘器h1后的波形’)ylabel(‘幅度’)
subplot(212)plot(t,sw2)
% title(‘经过相乘器h2后的波形’)ylabel(‘?幅度’)xlabel(‘t’)
bn=fir1(101,[2/80010/800])
figure(4)
st1=filter(bn,1,sw1)st2=filter(bn,1,sw2)
subplot(211)plot(t,st1)
% title(‘经过低通滤波器sw1后的波形’)ylabel(‘幅度’)
subplot(212)plot(t,st2)
% title(‘经过低通滤波器sw2后的波形’)ylabel(‘幅度’)xlabel(‘t’)
%判决
for(i=1:length(t))
if(st1(i)>=st2(i))st(i)=0
else st(i)=st2(i)
end
end
figure(5)
subplot(211)plot(t,st)
% title(‘经过抽样判决器后的波形’)ylabel(‘幅度’)
subplot(212)plot(t,sn)
% title(‘原始的波形’)ylabel(‘幅度’)xlabel(‘t’)
%加入噪声后
no=0.3*randn(1,t1)
sn=fsk+no
figure(6)
subplot(311)plot(t,no)
title('噪声波形')ylabel('噪声幅度') %噪声波形
subplot(312)plot(t,fsk)
title('信号')ylabel('信号幅度')
subplot(313)plot(t,sn)
title('叠加有噪声的信号')ylabel('幅度A')xlabel('时间t')
figure(7)
b1=fir1(101,[10/800 20/800])
b2=fir1(101,[90/800 110/800]) %设置带通参数
H1=filter(b1,1,sn)H2=filter(b2,1,sn) %经过带通滤波器后的信号
subplot(211)plot(t,H1)
title('经过带通滤波器后f1的波形')ylabel('幅度')
subplot(212)plot(t,H2)
title('经过带通滤波器后f2的波形')ylabel('幅度')xlabel('t')
sw1=H1.*H1sw2=H2.*H2 %经过相乘器
figure(8)
plot(t,sw2)
title('经过相乘器后h2的波形')ylabel('?幅度')xlabel('t')
bn=fir1(101,[2/800 10/800]) %经过低通滤波器
figure(9)
st1=filter(bn,1,sw1)st2=filter(bn,1,sw2)
plot(t,st1)
title('经过低通滤波器sw1后的波形')ylabel('幅度')xlabel('t')
%判决
for i=1:length(t)
if(st1(i)>=st2(i))
st1(i)=0
else st1(i)=st2(i)
end
end
figure(10)
subplot(211)plot(t,st)title('经过抽样判决后的波形')ylabel('幅度')
subplot(212)plot(t,sn)title('原始波形')ylabel('幅度')xlabel('t')
%程序完
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)