资源描述：
《信号与系统实验程序new》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、实验一熟悉Matlab软件环境1.1ones(3,2)zeros(3,2)rand(3,2)randn(3,2)eye(3)1.2求1～10的平均值、总和（1）mean.mfunction[m,s]=mean(a)l=length(a);s=sum(a);m=s/l;（2）a=1:9;[s,m]=mean(a)1.3y=sintt=0:1/100:2*pi;y=sin(t);plot(t,y)line([0,2*pi],[0,0])axis([02*pi-22])1.4y=exp(-0.1*t)*sin(t+1)t=0:1/100:2*pi;x=sin(t+1)
2、;z=exp(-0.1*t);y=x.*z;plot(t,y)line([0,2*pi],[0,0])axis([02*pi-2-2])实验二连续时间信号的描述及运算2.1s2chongji.mt0=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.2s2jieyue.mf=heaviside(
3、t)%heaviside单位阶跃函数plot(t,f)axis([0,3,-0.2,1.2])2.3s2zhishu.mt=0:0.05:3;f=exp(-1.5*t);plot(t,f)axis([0301.2])title('单边指数信号')text(3.1,0.05,'t')2.4s24.msymstf=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=
4、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=-fsubplot(2,3,6),ezplot(y5,[-3,3])实验三连续时间信号的描述及运算3.1s31.mb=[1]a=[218]subplot(1,2,1)impulse(b,a)%impluse绘出由向量a和b表示的连续系统在指定时间范围内的冲激响应h(t)的时域波形图
5、subplot(1,2,2)step(b,a)%step绘出由向量a和b表示的连续系统在指定时间范围内的阶跃响应g(t)的时域波形图3.2s32.ma=[121];b=[12];p=0.01;%定义采样间隔t=0:p:5;%定义时间间隔x=exp(-2*t);%定义输入信号lsim(b,a,x,t);%对系统输出信号进行仿真实验四连续时间信号的频域分析4.1s41.msymstfw;%定义符号变量f=heaviside(t+1)-heaviside(t-1);F=fourier(f);%符号函数f的傅里叶变换subplot(2,1,1);ezplot(f);su
6、bplot(2,1,2);ezplot(F);%ezplot绘制二维符号函数曲线4.2s42.mt=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)%为了绘制三维曲面，要把各次波形数据存为一个三维数组，因此必须重新定
7、义yy=zeros(10,max(size(t)));x=zeros(size(t));fork=1:2:19x=x+sin(k*t)/k;y((k+1)/2,:)=x;end%将各波形叠合绘出pause,figure(1),plot(t,y(1:9,:)),gridline([0,pi+0.5],[pi/4,pi/4])%加上方波幅度线及标注text(pi+0.5,pi/4,'pi/4')%将各半波形绘成三维网格图，看出增加谐波阶次对方波逼近程度的影响halft=ceil(length(t)/2);pause,figure(2),mesh(t(1:halft)
8、,[1:10],y(:,
