实验报告信号与系统

《实验报告信号与系统》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、信号与系统实验院系：水利电力学院班级：光伏二班小组成员：赵宝慧1300309042周洪煜1300309046张云才1300309043周生丽1300309045使用multisim完成实验实验一：二阶电路的瞬态响应电路图观察电感上的电流i的阶跃响应波形电路图电容上的阶跃响应波形图阶跃信号的设置仿真结果分析t=0:0.01:2*pi;f1=4/pi*sin(t);%基波f3=4/pi*(sin(3*t)/3);%三次谐波f5=4/pi*(sin(5*t)/5);f7=4/pi*(sin(7*t)/7);f9=4/pi*(sin(9*t)/9);y1=f1+f3;y2=f1+f3+f5;

2、y3=f1+f3+f5+f7+f9;subplot(2,2,1);plot(t,f1),holdony=1*sign(pi-t);plot(t,y,'c:');title('周期矩形波的形成-基波')subplot(2,2,2);plot(t,y1);holdon;y=1*sign(pi-t);plot(t,y,'c:');title('周期矩形波的形成-基波+3次谐波')subplot(2,2,3);plot(t,y2)holdon;y=1*sign(pi-t);plot(t,y,'c:');title('基波+3次谐波+5次谐波');subplot(2,2,4);plot(t,y

3、3);holdon;y=1*sign(pi-t);plot(t,y,'c:')title('-基波+3次谐波+5次谐波+7次谐波+9次谐波')实验二：方波信号的分解本实验通过MATLAB软件完成程序如下：t=0:0.01:2*pi;f1=4/pi*sin(t);%基波f3=4/pi*(sin(3*t)/3);%三次谐波f5=4/pi*(sin(5*t)/5);f7=4/pi*(sin(7*t)/7);f9=4/pi*(sin(9*t)/9);f11=4/pi*(sin(11*t)/11);f13=4/pi*(sin(13*t)/13);f15=4/pi*(sin(15*t)/15);

4、y1=f1+f3;y2=f1+f3+f5;y3=f1+f3+f5+f7+f9;y4=f1+f3+f5+f7+f9+f11;y5=f1+f3+f5+f7+f9+f11+f13;y6=f1+f3+f5+f7+f9+f11+f13+f15;subplot(2,2,1);plot(t,f1),holdony=1*sign(pi-t);plot(t,y,'c:');title('周期矩形波的形成-基波')subplot(2,2,2);plot(t,y4);holdon;y=1*sign(pi-t);plot(t,y,'c:');title('周期矩形波的形成-基波+11次谐波')subplot

5、(2,2,3);plot(t,y5)holdon;y=1*sign(pi-t);plot(t,y,'c:');title('基波+11次谐波+13次谐波');subplot(2,2,4);plot(t,y6);holdon;y=1*sign(pi-t);plot(t,y,'c:')title('-基波+11次谐波+13次谐波+15次谐波')要求对得到的图形结果记录并进行分析。试编写11次、13次、15次谐波的叠加程序。实验三：滤波器MATLAB实现：clear;b=[-1.000e-1000];%生成向量ba=[-1.000e-5-2.000e-5-1];%生成向量a[h,w]=fr

6、eqs(b,a,100)%求系统响应函数H(jw),设定100个频率点figure(1)subplot(3,1,1);plot(w,abs(h));%绘制幅频特性gridxlabel('w');ylabel('abs(h)');subplot(3,1,2);semilogx(w,20*log10(abs(h)));%绘制对数频率特性grid;xlabel('w');ylabel('分贝');subplot(3,1,3);plot(w,angle(h));%绘制相频特性grid;xlabel('w');ylabel('angle(h)');记录所得波形并分析。Multisim实现电路图

7、如下：观察所得波形，记录波形并进行分析。实验四：采样定理1、MATLAB实现连续信号的采样与重构原理wm=1;%信号带宽wc=wm;%滤波器截止频率,Ts=pi/wm;%采样间隔ws=2*pi/Ts;%采样角频率n=-100:100;%时域采样点数nTs=n*Ts;f=sinc(nTs/pi);Dt=0.005;t=-15:Dt:15;fa=f*Ts*wc/pi*sinc((wc/pi)*(ones(length(nTs),1)*t-nTs'*ones(1