信号上机实验三.docx

1、信号上机实验三 信号上机实验三一、实验目的1、掌握信号时域运算的MATLAB实现方法。2、掌握信号时域变换的MATLAB实现方法二、实验原理信号的时域运算包括信号的相加、相乘，信号的时域变换包括信号的平移、反折、倒相及信号的尺度变换三、实验内容1、将实验原理中提到的例子在计算机上全部练习一遍；解:例3-2. 求下列两序列的和序列： f1(k)=-2,-1,0,1,2, f2(k)=1,1,1,其MATLAB命令为functionf,k=isxj(f1,f2,k1,k2)k=min(min(k1),min(k2):max(max(k1),max(k2);s1=zeros(1,length(k);

2、s2=s1;s1(find(k=min(k1)&(k=min(k2)&(k=max(k2)=1)=f2;f=s1+s2;stem(k,f,filled)axis(min(min(k1),min(k2)-1),(max(max(k1),max(k2)+1),(min(f)-0.5),(max(f)+0.5)例3-3：己知f(k)=2k -3k3，画出f(-k)波形。functionf,k=isfz(f1,k1)f=fliplr(f1);k=-fliplr(k1);stem(k,f,filled)axis(min(k)-1,max(k)+1,min(f)-0.5,max(f)+0.5)例3-4.

3、己知f(k)=k2 -4k4，画出f(k-2)波形。functionf,k=isyw(ff,kk,k0)k=kk+k0;f=ff;stem(k,f,filled)axis(min(k)-1,max(k)+1,min(f)-0.5,max(f)+0.5)2、已知信号，画出的波形；解;3、已知两个连续信号，用MATLAB绘出下列信号的波形：1） 2）3） 4）解;4、离散序列，用MATLAB绘出下列序列的波形。1） 2）3） 4）解：（1）（2）（3）(4)四、实验思考题在对信号进行平移、反折和尺度变换时，运算顺序对结果是否有影响？在运算中应该注意什么？解：运算顺序对结果没影响。为了画f（-t-t

4、0）或f（-k-k0），最好是先平移，将f(t)平移为f(t-t0)或 f(k)平移为f(k-k0),然后再反转，将变量t和k相应地换为-t或-k。如果反转后再进行平移，由于这时自变量变为-t或-k，故平移方向与前述相反。信号上机实验五一、实验目的1.用MATLAB实现周期信号傅里叶级数分解与综合。2.利用MATLAB对周期信号进行频谱分析二、实验原理任何满足狄里赫里条件的周期信号，都可以表示成三角形式或指数形式的傅里叶级数展开式。 （5.1） （5.2）一般来说，傅里叶级数有无限个非零值，即任何具有有限个间断点的周期信号都一定有一个无限非零系数的傅里叶级数表示。但对数值计算来说，这是无法实现

5、的。在实际应用中，我们可以用有限项的傅里叶级数求和来逼近。即对有限项和： （5.3）当N值取得较大时，上式就是原周期信号的一个很好的近似。上式常称作截断傅里叶级数表示。三、实验内容1.将例5-1在计算机上练习一遍，观察运行结果是否正确；解;function A_sym,B_sym=CTFShchsymsyms t n k x T=5;tao=0.2*T;a=0.5;if nargin4;Nf=6;endif nargin5;Nn=32;endx=time_fun_x(t);A0=2*int(x,t,-a,T-a)/T; As=int(2*x*cos(2*pi*n*t/T)/T,t,-a,T-a

6、); Bs=int(2*x*sin(2*pi*n*t/T)/T,t,-a,T-a); A_sym(1)=double(vpa(A0,Nn); for k=1:NfA_sym(k+1)=double(vpa(subs(As,n,k),Nn);B_sym(k+1)=double(vpa(subs(Bs,n,k),Nn); end if nargout=0 c=A_sym;disp(c) d=B_sym;disp(d) t=-8*a:0.01:T-a;f1=0.4/2+0.3742.*cos(2*pi*1*t/5)+0.*sin(2*pi*1*t/5); f2=0.3027.*cos(2*pi*2*

7、t/5)+0.*sin(2*pi*2*t/5); f3=0.2018.*cos(2*pi*3*t/5)+0.*sin(2*pi*3*t/5); f4=0.0935.*cos(2*pi*4*t/5)+0.*sin(2*pi*4*t/5); f5=-0.0624.*cos(2*pi*6*t/5)+0.*sin(2*pi*6*t/5); f6=f1+f2; f7=f6+f3; f8=f7+f4+f5; subplot(2,2,1)plot(t,f1),hold ony=time_fun_e(t) plot(t,y,r:)title(周期矩形波的形成基波)axis(-4,4.5,-0.5,1.3)su

8、bplot(2,2,2)plot(t,f6),hold ony=time_fun_e(t)plot(t,y,r:)title(周期矩形波的形成基波+2次谐波)axis(-4,4.5,-0.5,1.3)subplot(2,2,3)plot(t,f7),hold ony=time_fun_e(t)plot(t,y,r:)title(基波+2次谐波+3次谐波)axis(-4,4.5,-0.5,1.3)subplot(2,2,4)plot(t,f8),hold ony=time_fun_e(t)plot(t,y,r:)title(基波+2次谐波+3次谐波+4次谐波+6次谐波)axis(-4,4.5,-

9、0.5,1.3)endfunction x=time_fun_x(t)h=1; x1=sym(Heaviside(t+0.5)*h;x=x1-sym(Heaviside(t-0.5)*h;function y=time_fun_e(t)a=0.5;T=5;h=1;tao=0.2*T;t=-8*a:0.01:T-a;e1=1/2+1/2.*sign(t+tao/2);e2=1/2+1/2.*sign(t-tao/2);y=h.*(e1-e2); 2.修改例5-1，将周期信号分解为前10次谐波的叠加，观察运行结果，能得出什么结论？解：function A_sym,B_sym=CTFShchsym1

10、syms t n k x T=5;tao=0.2*T;a=0.5;if nargin4;Nf=9;endif nargin0); （2）要利用sym来创建符号表达式。采用fourier( )及ifourier( )得到的返回函数，仍然是符号表达式。若需对返回函数作图，则应用ezplot( )绘图命令而不能用polt( )命令。如果返回函数中含有诸如狄拉克函数等的项，则用ezplot( )也无法作图。此外，用fourier( )对某些信号求变换时，其返回函数可能会包含一些不能直接表达的式子，甚至可能会出现一些屏幕提示为：“未被定义的函数或变量”的项，此时更不用说对此返回函数作图了。这是fouri

11、er( )的一个局限。另一局限是在很多应用场合，原信号尽管是连续的，但却不能表示成符号表达式，而更多的实际测量现场获得的信号是多给离散的数值量，此时也不能应用fourier( )对离散信号进行处理，而只能采用数值计算的方法求得（参阅参考文献2）。二.傅里叶变换性质时移性质若，则傅里叶变换的时移特性为： （三、实验内容1、将实验原理中的例子在计算机上全部练习一遍；例6-1. 求的傅里叶变换。解：例6-2. 求的傅里叶逆变换。例6-3.设，试画出其幅度频谱。解：例4：已知信号，试用MATLAB画出该信号及其频谱（幅度谱及相位谱）解：clear r=0.02;t=-5:r:5;N=200;W=2*p

12、i*1;k=-N:N;w=k*W/N;f1=1/2*exp(-2*t).*Heaviside(t);F=r*f1*exp(-j*t*w);F1=abs(F);PI=angle(F);subplot(3,1,1);plot(t,f1);grid;xlabel(t);ylabel(f(t);title(f(t);subplot(3,1,2); plot(w,F1);xlabel(w);grid;ylabel(F(jw);subplot(3,1,3);plot(w,PI*180/pi);grid;xlabel(w);ylabel(P(度);2、已知信号f(t)=e-2(t+0.3)(t+0.3)，试

13、画出该信号波形及其频谱；解：r=0.02;t=-5:r:5;N=200;W=2*pi*1;k=-N:N;w=k*W/N;f1=1/2*exp(-2*(t+0.3).*Heaviside(t+0.3);F=r*f1*exp(-j*t*w);F1=abs(F);PI=angle(F);subplot(3,1,1);plot(t,f1);grid;xlabel(t);ylabel(f(t);title(f(t);subplot(3,1,2); plot(w,F1);xlabel(w);grid;ylabel(F(jw);subplot(3,1,3);plot(w,PI*180/pi);grid;xl

14、abel(w);ylabel(P(度);3、已知信号f(t)=e-2(t-0.3)(t-0.3)，试画出该信号波形及其频谱。解：r=0.02;t=-5:r:5;N=200;W=2*pi*1;k=-N:N;w=k*W/N;f1=1/2*exp(-2*(t-0.3).*Heaviside(t-0.3);F=r*f1*exp(-j*t*w);F1=abs(F);PI=angle(F);subplot(3,1,1);plot(t,f1);grid;xlabel(t);ylabel(f(t);title(f(t);subplot(3,1,2); plot(w,F1);xlabel(w);grid;yla

15、bel(F(jw);subplot(3,1,3);plot(w,PI*180/pi);grid;xlabel(w);ylabel(P(度);四、实验思考题观察练习2、3的结果与例4的结果做比较，可得出什么结论？解：时幅度谱保持不变，相位谱随f（t）的变化而变化。信号上机实验七一、实验目的1.了解信号幅度调制的原理；2.能够用MATLAB实现信号的幅度调制。二、实验原理 设信号的频谱为，现将乘以载波信号，得到高频的已调信号，即 成为调制信号。实现信号调制的原理图如图7.1所示。 图7.1 幅度调制原理图从频域上看，已调信号的频谱为原调制信号的频谱搬移到处，幅度降为原的1/2，即： 注意：这里采用

16、的调制方法为抑制载波方式，即的频谱中不含有的频率分量。MATLAB提供了专门的函数modulate( )用于是吸纳信号的调制。调用格式为：y=modulate(x,Fc,Fs,method)y,t=modulate(x,Fc,Fs)其中，x为被调信号，Fc为载波频率，Fs为信号x的采样频率，method为所采用的调制方式，若采用幅度调制、双边带调制、抑制载波调制，则method为am或amdsd-sc。其执行算法为：y=x*cos(2*pi*Fc*t)其中y为已调制信号，t为函数计算时间间隔向量三、实验内容1、将实验原理中的例子在计算机上练习一遍例7-1：设信号，载波为频率为400Hz的余弦信

17、号。试用MATLAB实现调幅制信号，并观察的频谱和的频谱，以及两者在频域上的关系。解：Fs=1000;Fc=400;N=1000;n=0:N-2;t=n/Fs;x=sin(2*pi*50*t);subplot(221)plot(t,x);xlabel(t(s);ylabel(x);title(被调信号);axis(0 0.1 -1 1)Nfft=1024;window=hamming(512);noverlap=256;dflag=none;Pxx,f=psd(x,Nfft,Fs,window,noverlap,dflag);subplot(222)plot(f,Pxx)xlabel(频率(H

18、z);ylabel(功率谱(X);title(被调信号的功率谱)gridy=modulate(x,Fc,Fs,am);subplot(223)plot(t,y)xlabel(t(s);ylabel(y);axis(0 0.1 -1 1)title(已调信号)Pxx,f=psd(y,1024,Fs,window,noverlap,dflag);subplot(224)plot(f,Pxx)xlabel(频率(Hz);ylabel(功率谱(Y);title(已调信号的功率谱);grid例7-2：设，试用MATLAB画出、的时域波形及其频谱，并观察傅立叶变换的频移特性。解：R=0.005;t=-1.

19、2:R:1.2;f=Heaviside(t+1)-Heaviside(t-1);f1=f.*cos(10*pi*t);subplot(221)plot(t,f)xlabel(t);ylabel(f(t);subplot(222);plot(t,f1);xlabel(t);ylabel(f1(t)=f(t)*cos(10*pi*t);W1=40;N=1000;k=-N:N;W=k*W1/N;F=f*exp(-j*t*W)*R; F=real(F);F1=f1*exp(-j*t*W)*R; F1=real(F1);subplot(223);plot(W,F);xlabel(w);ylabel(F(

20、jw);subplot(224);plot(W,F1);xlabel(w);ylabel(F1(jw);2.已知信号载波为频率为400Hz的余弦信号。试用MATLAB实现调幅制信号，并观察的频谱和的频谱，以及两者在频域上的关系。解：Fs=1000;Fc=400;N=1000;n=0:N-2;t=n/Fs;x=sin(2*pi*t)/(2*pi*t);subplot(2,2,1)plot(t,x);xlabel(t(s);ylabel(x);title(被调信号);axis(0 0.1 -1 1)Nfft=1024;window=hamming(512);noverlap=256;dflag=none;Pxx,f=psd(x,Nfft,Fs,window,noverlap,dflag);subplot(222)plot(f,Pxx)