实验报告信号与系统.docx

信号与系统实验 院系：水利电力学院 班级：光伏二班 小组成员：赵宝慧 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; y3=f1+f3+f5+f7+f9;subplot(2,2,1);plot(t,f1),hold ony=1*sign(pi-t);plot(t,y, c:); title(周期矩形波的形成-基波)subplot(2,2,2);plot(t,y1); hold on;y=1*sign(pi-t);plot(t,y, c:); title(周期矩形波的形成-基波+3次谐波)subplot(2,2,3);plot(t,y2)hold on;y=1*sign(pi-t);plot(t,y, c:); title(基波+3次谐波+5次谐波);subplot(2,2,4) ;plot(t,y3);hold on;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);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),hold ony=1*sign(pi-t);plot(t,y, c:); title(周期矩形波的形成-基波)subplot(2,2,2);plot(t,y4); hold on;y=1*sign(pi-t);plot(t,y, c:); title(周期矩形波的形成-基波+11次谐波)subplot(2,2,3);plot(t,y5)hold on;y=1*sign(pi-t);plot(t,y, c:); title(基波+11次谐波+13次谐波);subplot(2,2,4) ;plot(t,y6);hold on;y=1*sign(pi-t);plot(t,y, c:)title(-基波+11次谐波+13次谐波+15次谐波)要求对得到的图形结果记录并进行分析。试编写11次、13次、15次谐波的叠加程序。实验三： 滤波器matlab实现：clear; b=-1.000e-10 0 0; %生成向量ba=-1.000e-5 -2.000e-5 -1;%生成向量ah,w=freqs(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 实现电路图如下：观察所得波形，记录波形并进行分析。实验四： 采样定理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,length(t); %信号重构error=abs(fa-sinc(t/pi); %求重构信号与原信号的误差t1=-15:0.5:15;f1=sinc(t1/pi);subplot(3 ,1 ,1);stem(t1,f1);xlabel(kts);ylabel(f(kts);title(sa(t)=sinc(t/pi)的临界采样信号);subplot(3 ,1 ,2);plot(t,fa);xlabel(t);ylabel(fa(t);grid;title(由sa(t)=sinc(t/