《信号与系统实验程序.docVIP

实验一 熟悉Matlab软件环境 1.1 ones(3,2) zeros(3,2) rand(3,2) randn(3,2) eye(3) 1.2 求1～10的平均值、总和 （1）mean.m function[m,s]=mean(a) l=length(a); s=sum(a); m=s/l; （2） a=1:9; [s,m]=mean(a) 1.3 y=sint t=0:1/100:2*pi; y=sin(t); plot(t,y) line([0,2*pi],[0,0]) axis([0 2*pi -2 2]) 1.4 y=exp(-0.1*t)*sin(t+1) t=0:1/100:2*pi; x=sin(t+1); z=exp(-0.1*t); y=x.*z; plot(t,y) line([0,2*pi],[0,0]) axis([0 2*pi -2 -2]) 实验二 连续时间信号的描述及运算 2.1 s2chongji.m t0=0; tf=5; dt=0.01; t1=1; t=[t0:dt:tf]; st=length(t); n1=floor((t1-t0)/dt); %求t1对应的样本序号 x1=zeros(1,st); %把全部信号先初始化为零 x1(n1)=1/dt; %给出t1处的脉冲信号 stairs(t,x1) %绘图 axis([0,5,0,22]) 2.2 s2jieyue.m f=heaviside(t) %heaviside单位阶跃函数 plot(t,f) axis([0,3,-0.2,1.2]) 2.3 s2zhishu.m t=0:0.05:3; f=exp(-1.5*t); plot(t,f) axis([0 3 0 1.2]) title(单边指数信号) text(3.1,0.05,t) 2.4 s24.m syms t f=sym((t/2+1)*(heaviside(t+2)-heaviside(t-2))) subplot(2,3,1),ezplot(f,[-3,3]) y1=subs(f,t,t+2) subplot(2,3,2),ezplot(y1,[-5,1]) y2=subs(f,t,t-2) subplot(2,3,3),ezplot(y2,[-1,5]) y3=subs(f,t,-t) subplot(2,3,4),ezplot(y3,[-3,3]) y4=subs(f,t,2*t) subplot(2,3,5),ezplot(y4,[-2,2]) y5=-f subplot(2,3,6),ezplot(y5,[-3,3]) 实验三 连续时间信号的描述及运算 3.1 s31.m b=[1] a=[2 1 8] subplot(1,2,1) impulse(b,a) %impluse绘出由向量a和b表示的连续系统在指定时间范围内的冲激响应h(t)的时域波形图 subplot(1,2,2) step(b,a) %step绘出由向量a和b表示的连续系统在指定时间范围内的阶跃响应g(t)的时域波形图 3.2 s32.m a=[1 2 1]; b=[1 2]; p=0.01;%定义采样间隔 t=0:p:5;%定义时间间隔 x=exp(-2*t);%定义输入信号 lsim(b,a,x,t);%对系统输出信号进行仿真 实验四 连续时间信号的频域分析 4.1 s41.m syms t f w ;%定义符号变量 f=heaviside(t+1)-heaviside(t-1); F=fourier(f);%符号函数f的傅里叶变换 subplot(2,1,1); ezplot(f); subplot(2,1,2); ezplot(F);%ezplot绘制二维符号函数曲线 4.2 s42.m t=0:.01:2*pi;%设定一个时间数组，有101个点 y=sin(t); plot(t,y),figure(gcf),pause %频率为w=1(f=1/2π)基波 y=sin(t)+sin(3*t)/3; plot(t,y),pause %叠加三次谐波 y=sin(t)+sin(3*t)/3+sin(5*t)/5+sin(7*t)/7+sin(9*t)/9;%用1，3，5，7，9次谐波叠加 plot(t,y) %为了绘制三维曲面，要把各次波形数据存为一个三维数组，因此必须重新定义y y=zeros(10,max(size(t))); x=zeros(size(t)); for k=1:2:19 x=x+sin(k*t)/k;