MATLAB实习报告.docx

1、MATLAB实习报告 西安邮电学院信号与系统实验 报 告院系名称：电子工程学院院 专业班级：集成电路设计与集成系统0801学生姓名：赵启学号：04086017时间：2010年 6 月 20 日 实验一、熟悉MATLAB环境 A=1 2 3;B=4 5 6;C=A+B;D=A-B;E=A.*B;F=A./B;G=A.B;subplot(2,4,1),stem(A)subplot(2,4,2),stem(B)subplot(2,4,3),stem(C)subplot(2,4,4),stem(D)subplot(2,4,5),stem(E)subplot(2,4,6),stem(F)subplot(

2、2,4,7),stem(G)t=0:0.001:10;f=4+exp(-t);plot(t,f)plot(t,f)t=0:0.001:10;f=2*exp(-2*t)+3*exp(-3*t);plot(t,f) t=-4:0.001:4;f=cos(2*pi*t+pi/3);plot(t,f)t=0:0.01:3;f=exp(-t).*sin(2*pi*t);plot(t,f) t=-6:0.01:6;f=heaviside(t+3)-heaviside(t-3);plot(t,f)t=-10:0.001:10;f=et(t+3)-et(t-3);plot(t,f)t=-4*pi:0.001:

3、4*pi;f=sinc(t/pi);plot(t,f)grid on 实验二、离散时间信号的基本运算1、A: 单位采样序列function dwxulie(k1,k2,k0)k=k1:k2;n=length(k);f=zeros(1,n);f(1,-k0-k1+1)=1;stem(k,f,filled)axis(k1,k2,0,1.5)title(单位序列d(k)B:单位阶跃序列function jyxulie(k1,k2,k0)k=k1:-k0-1;kk=-k0:k2;n=length(k);nn=length(kk)u=zeros(1,n); %k0前信号赋值为零uu=ones(1,nn)

4、;%k0后信号赋值为一stem(kk,uu,filled)hold onstem(k,u,filled)hold offtitle(单位阶跃序列)axis(k1 k2 0 1.5)C:实指数虚列function dszsu(c,a,k1,k2)k=k1:k2;x=c*(a.k);stem(k,x,filled)hold onplot(k1,k2,0,0)hold offD:复数指数虚列function dfzsu(n1,n2,r,w)k=n1:n2;f=(r*exp(i*w).k;Xr=real(f);Xi=imag(f);Xa=abs(f);Xn=angle(f);subplot(2,2,1

5、), stem(k,Xr,filled),title(实部);subplot(2,2,3), stem(k,Xi,filled),title(虚部);subplot(2,2,2), stem(k,Xa,filled),title(模);subplot(2,2,4), stem(k,Xn,filled),title(相角);E：正余弦序列function yuxuan(w,o)n=0:40;f=cos(w*n+o);stem(n,f,filled)F：周期序列function f=zhouqi(f,T)n=0:40;f=f(mod(n,T)+1);stem(n,f);2、实现IMPSEQ(N0,

6、N1,N2)function jyxulie(k1,k2,k0)k=k1:-k0-1;kk=-k0:k2;n=length(k);nn=length(kk)u=zeros(1,n); uu=ones(1,nn);stem(kk,uu,filled)hold onstem(k,u,filled)hold offtitle()axis(k1 k2 0 1.5)4、两信号相加functionf,k=lsxj(f1,f2,k1,k2)k=min(min(k1),min(k2):max(max(k1),max(k2); s1=zeros(1,length(k);s2=s1;s1(find(k=min(k

7、1)&(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)；11、t0=0:0.5:1;f=sin(2*pi.*t);subplot(1,3,1);stem(t0,f);t1=0:0.25:1;f=sin(2*pi.*t1);subplot(1,3,2);stem(t1,f);t2=0:0.1:1;f=sin(2*pi.*t2);subplot(1,3,3);stem(t2,f);实验三

8、、系统零状态响应的时域求解1、a=1 2 100;b=0 0 1;t=-10:0.1:10;subplot(1,3,1);impulse(b,a);subplot(1,3,2);step(b,a);f=10*sin(2*pi.*t);subplot(1,3,3);lsim(b,a,f,t);3、a=1 -1 0.9 -0.5;b=1;subplot(1,2,1);impz(b,a);subplot(1,2,2);dstep(b,a);4、（1）、x=3 11 7 0 -1 4 2;h=2 3 0 -5 2 1;conv(x,h)（2）、function y,nyconv_m(x,nx,h,nh

9、);y=conv(x,h);ny0=nx(1)+nh(1); ny3=length(x)+length(h)-2;ny=ny0:nyo+ny3;end实验四、周期信号频谱分析与综合1、function A,B,C,fn,t,x=CTFSfftT=5;if (nargin2 | isempty(M);M=8;endif nargin3;Nf=input(plear Input :Nf=);endN=2M; f=1/T; w0=2*pi*f;dt=T/N; n=0:1:(N-1); t=n*dt; x=time_fun(t,T); X=fft(x); cn=X/N; z_cn=find(abs(c

10、n)N/2;error( Nf int2str(N/2-1);endA(Nf+2:end)=;B(Nf+2:end)=;n1=-N/2*2*pi/T:2*pi/T:N/2*2*pi/T;fn=n1*f; y=time_fun_e subplot(2,1,2),stem(n1,abs(C);title()axis(-150,150,0,0.12)hold off%-function x=time_fun(t,T)A=1;tao=1;T=5;x=zeros(size(t);ii=find(t=-tao /2& t=tao/2);x(ii)=ones(size(ii).*A;x(t=0)=1;%-f

11、unction y=time_fun_eT=5;t=-2*T:0.01:2*T;tao=T/5;y=rectpuls(t,1); subplot(2,2,2)plot(t,y)hold ony=rectpuls(t-5,1); plot(t,y)hold ony=rectpuls(t+5,1); plot(t,y)title(T=5tao=1 )axis(-10,10,0,1.2)4、functionA_sym,B_sym=CTFShchsymsyms t n k xT=5;tao=0.2*T;a=0.5;if nargin4;Nf=6;endif nargin5;Nn=32;endx=tim

12、e_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); 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);endif nargout=0c=A_sym;disp(c) d=B_sym;disp(d) t=-8*a:0.01:T-a;f1=0.2

13、/2+0.1871.*cos(2*pi*1*t/5)+0.*sin(2*pi*1*t/5); f2=0.1514.*cos(2*pi*2*t/5)+0.*sin(2*pi*2*t/5); f3=0.1009.*cos(2*pi*3*t/5)+0.*sin(2*pi*3*t/5); f4=0.0468.*cos(2*pi*4*t/5)+0.*sin(2*pi*4*t/5); f5=-0.0312.*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 o

14、ny=time_fun_e % plot(t,y)title()subplot(2,2,2)plot(t,f6),hold ony=time_fun_eplot(t,y)title(+2 )subplot(2,2,3)plot(t,f7),hold ony=time_fun_eplot(t,y)title(+2 +3 )subplot(2,2,4)plot(t,f8),hold ony=time_fun_eplot(t,y)title(+2 +3 +4 +6 )end%-function x=time_fun_x(t)h=1;x1=sym(Heaviside(t+0.5)*h;x=x1-sym

15、(Heaviside(t-0.5)*h;%-function y=time_fun_ea=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); 实验五、零极点分析及系统的频响特性1、num=3 -9 6;den=1 3 2;z,p=tf2zp(num,den);subplot(1,3,1);zplane(z,p);c=1;d=1 0;z,p=tf2zp(c,d);subplot(1,3,2);zplane(z,p);e=1 1;f=1 2 5

16、;z,p=tf2zp(e,f);subplot(1,3,3);zplane(z,p);2、（1）调用function ljdt(A,B)p=roots(A); q=roots(B); p=p;q=q; x=max(abs(p q 1); x=x+0.1;y=x;clfhold onaxis(-x x -y y) w=0:pi/300:2*pi;t=exp(i*w);plot(t) axis(square)plot(-x x,0 0)plot(0 0,-y y) text(0.1,x,jImz)text(y,1/10,Rez)plot(real(p),imag(p),x) plot(real(q

17、),imag(q),o) title(pole-zero diagram for discrete system) hold offnum=1 0 0;denn=1 1 0.5;ljdt(num,denn);（2）、num=1 0.5 0 0;den=1 -1.25 0.75 -0.125;ljdt(num,den);3、（1）clear allclose alla=2 0;b=1 sqrt(2) 1;z,p=tf2zp(a,b);w=logspace(-1,1);m=0;n=;f1=0;f2=10;k=0.1;subplot(1,3,1);zplane(z,p);subplot(1,3,2)

18、;impulse(a,b);subplot(1,3,3);splxy(f1,f2,k,m,n);（2）clear allclose alla=3 -9 2;b=1 3 2;z,p=tf2zp(a,b);w=logspace(-1,1);m=a;n=;f1=0;f2=10;k=0.1;subplot(1,3,1);zplane(z,p);subplot(1,3,2);impulse(a,b);subplot(1,3,3);splxy(f1,f2,k,m,n);实验六、信号的调制与解调1、2、n=0:256;s=sin(2*pi*n)/256);sm=demod(s,100000,1000000,amssb);stem(n,sm)