利用MATLAB软件分析系统频率响应.docVIP

康学院电子与信息工程系——信号与系统实验报告 书山有路勤为径，学海无涯苦作舟！ 任课教师 余顺园 2012年 第 PAGE 5页 实验名称专业班级学 号姓 名利用MATLAB软件分析系统的频率响应备注：（1）、按照要求独立完成实验项目内容，报告中要有程序代码和程序运行结果和波形图等原始截图。 （2）、实验结束后，把电子版实验报告按要求格式改名（例：09号-张三-实验一）后，上传至指定ftp服务器目录下（homework_upload）的相应文件里，并由实验教师批阅记录后；实验室统一刻盘留档。 ftp:6 账号：microele 密码：ele1507 实验七 利用MATLAB软件分析系统的频率响应 一、实验目的： 1、利用MATLAB求解系统的频率响应。 二、实验原理 MATLAB提供了函数freqs来计算系统的频率响应。 三、实验内容：（包括代码与产生的图形） 6-16 w=linspace(0,5,200); b=[1]; a=[1 2 2 1]; H=freqs(b,a,w); subplot(2,1,1); plot(w,abs(H)); set(gca,xtick,[0 1 2 3 4 5]); set(gca,ytick,[0 0.4 0.7071]);grid; xlabel(\omega) subplot(2,1,2); plot(w,angle(H)); set(gca,xtick,[0 1 2 3 4 5]);grid; xlabel(\omega); 6-17 RC=0.04; t=linspace(-2,2,1024); w1=5;w2=100; H1=j*w1/(j*w1+1/RC); H2=j*w2/(j*w2+1/RC); f=cos(5*t)+cos(100*t); y=abs(H1)*cos(w1*t+angle(H1))+abs(H2)*cos(w2*t+angle(H2)); subplot(2,1,1); plot(t,f); subplot(2,1,2); plot(t,y); 6-1 w=linspace(0,5,200); b=[10]; a=[1 10]; H=freqs(b,a,w); plot(w,abs(H)); RC=0.04; t=linspace(0,0.2,1024); w1=100;w2=3000; H1=(1/RC)/(j*w1+1/RC); H2=(1/RC)/(j*w2+1/RC); f=cos(100*t)+cos(3000*t); y=abs(H1)*cos(w1*t+angle(H1))+abs(H2)*cos(w2*t+angle(H2)); subplot(2,1,1); plot(t,f); subplot(2,1,2); plot(t,y);