matlab部分程序.doc

1. 用MATLAB命令产生如下信号，并绘出波形图 (1) t=-2*pi:pi/50:2*pi; a=t=0; f(a)=2*exp(-1.5*t(a)); f(~a)=0; plot(t,f) （2） t=-2*pi:pi/50:2*pi; y=2*sin(2*pi*t+pi/4); plot(t,y) (3) t=-5:0.01:5; y=2*rectpuls(t-0.5,1); plot(t,y) axis equal 2. 已知信号的波形如图所示，试用MATLAB命令画出、、、的波形图。 function y=p(t) y1=(t0t=-2); y2=(-t+1).*(t=1t=0); y=y1+y2 t=-5:0.01:5; y=p(t-2); plot(t,y) axis equal t=-5:0.01:5; y=p(3*t); plot(t,y) axis equal t=-5:0.01:5; y=p(-t); plot(t,y) axis equal t=-5:0.01:5; y=p(-3*t-2); plot(t,y) axis equal 3.已知描述系统的微分方程和激励信号如下 ， 要求：（1）从理论上求解系统的冲激响应和零状态响应，并根据求解结果用MATLAB绘制其时域波形； （2）分别用MATLAB的impulse()函数和lsim()函数绘制系统的冲激响应和零状态响应，验证（1）中的结果。 clf; a=[1 4 4]; b=[0 1 3]; subplot(2,1,1); impulse(b,a) t=0:0.01:3.5; h=(1+t).*exp(-2*t); subplot(2,1,2); plot(t,h) axis([0,3.5,0,1]) figure(2) t=0:0.01:3.5; x=exp(-t); subplot(2,1,1); lsim(b,a,x,t) y=2*exp(-t)-(2+t).*exp(-2*t); subplot(2,1,2); plot(t,y) axis([0,3.5,0,1]) 4.已知周期矩形脉冲信号的波形如图所示，设脉冲幅度为1，宽度为，周期为，用MATLAB绘制信号的频谱图，并研究脉冲宽度变化时（分别取、；、；、），对其频谱的影响。 t=4; k=-10:0.01:10; T=8; w=2*pi/T; y=1/2*sin(k*w*t)./(k*w*t) plot(k*w,y) t=2; k=-10:0.01:10; T=8; w=2*pi/T; y=1/2*sin(k*w*t)./(k*w*t) plot(k*w,y) t=1; k=-10:0.01:10; T=8; w=2*pi/T; y=1/2*sin(k*w*t)./(k*w*t) plot(k*w,y) 5.设系统的频率响应为，若外加激励信号为，用MATLAB求其响应。 a=[1 3 2]; b=[0 0 1]; t=-10:0.01:10; x=5*cos(t)+2*cos(10*t); y=lsim(b,a,x,t); plot(t,y,t,x) 6.下图是用RLC元件构成的二阶低通滤波器。设，，，试用MATLAB的freqs()函数绘出该系统的幅频特性曲线和相频特性曲线。（求模：abs()函数，求相角：angle()函数） a=[0.08 0.4 1]; b=[0 0 1]; w = logspace(-1,1); freqs(b,a,w) 7.假设基带信号为，被调制成，在接收端又被解调为，并通过低通滤波器 得到信号。试用MATLAB绘制上述各个信号的时域波形和频谱。 function [t,omg,FT,IFT]=prefourier(Trg,N,OMGrg,K) T=Trg(2)-Trg(1); t=linspace(Trg(1),Trg(2)-T/N,N); OMG=OMGrg(2)-OMGrg(1); omg=linspace(OMGrg(1),OMGrg(2)-OMG/K,K); FT=T/N*exp(-j*kron(omg,t)); IFT=OMG/2/pi/K*exp(j*kron(t,omg)); clc,clear,clf [t,omg,FT,IFT]=prefourier([-5 5],length(-5:0.01:5),[-250 250],2000); g=3*cos(10*t)+2*cos(20*t); G=FT*g; figure(1); subplot(2,1,1) plot(t,g) title(g(t)=3cos(10t)