




下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、信号与系统第二次实验(题目二)2、给定一个连续LTI系统,其微分方程为:输入信号: (其中=0.06Hz,=5Hz,=32Hz) u 理论计算:1) 计算系统的幅度响应,判断该系统是哪一类频率选择性滤波器?(低通、高通、带通、带阻)2) 求系统的输出信号。u 仿真分析:程序如下:ts=0; %起始时刻te=50; %终止时刻fs=100; %采样频率t=ts:1/fs:te;N=1000;x=sin(2*pi*0.06*t)+sin(2*pi*5*t)+sin(2*pi*32*t);%信号描述X=fft(x,N)/N; fx=(-(N-1)/2:(N-1)/2)/N*fs; %建立频率轴num
2、=1 1; %分子系数den=1 110 1000; %分母系数sys=tf(num,den); %建立系统syssubplot(334);plot(t,x);xlabel(时间);ylabel(时域);title(x(t)时域波形);mag=abs(fftshift(X);theta=angle(fftshift(X)/pi*180;%修正幅度谱%delta_1=1,2*ones(1,length(fftshift(X)-1);mag = mag.*delta_1; %修正后的幅度谱%修正相位谱%delta_2= (mag0.1); % 将幅值为0的频率分量的相位置为0”theta = th
3、eta.*delta_2;%修正后的相位谱subplot(335);plot(fx,mag);xlabel(频率/Hz);ylabel(幅值);title(x(t)幅频谱);subplot(336);bar(fx,theta,3); xlabel(频率/Hz);ylabel(相位);title(x(t)相频谱);y=lsim(sys,x,t); %将输入信号x加在系统sys上,得到输出信号ysubplot(337);plot(t,y);xlabel(时间);ylabel(时域);title(y(t)时域波形);Y=fft(y,N)/N; mag=abs(fftshift(Y); theta=(
4、angle(fftshift(Y)/pi*180);%修正幅度谱%delta_1=1,2*ones(1,length(fftshift(Y)-1);mag = mag.*delta_1; %修正后的幅度谱%修正相位谱%delta_2=(mag0.0005); % 将幅值为0的频率分量的相位置为0”theta = theta.*delta_2;%修正后的相位谱subplot(338);plot(fx,mag); xlabel(频率/Hz);ylabel(幅值);title(y(t)幅频谱);subplot(339);bar(fx,theta,3); xlabel(频率/Hz);ylabel(相位
5、);title(y(t)相频谱); %画bode图%w=0:0.1:10000; W=logspace(-1,2,100); mag,pha=bode(num,den,w); magdB=20*log10(mag);set(gcf,color,w) %设置背景颜色为白色subplot(321);semilogx(w,magdB);grid on;xlabel(omega);ylabel(20log_1_0|H(jomega)|);title(伯德图(a) 模);subplot(322);semilogx(w,pha);grid on;set(gca,yticklabel,-pi/2,-pi/4,0,pi/4,pi/2)xlabel(omega);ylabel(argH(jomega);title(伯德图(b) 相位);1) 绘制系统的频率响应Bode图;频率响应Bode图如下:2) 绘制输入、输出信号的时域波形和频谱,分析和解释滤波过程;输入信号x(t)时域波形及频谱图如下:输出信号y(t)时域波形及频谱图如下:分析和解释滤波过程:根据计算和bode图可以看出3) 对输入信号的采样时长改变时,频谱怎样变化?为什么?采样频率为50、200时频谱如下:由图可看出,采样时长变大,采样点减少,分散度变大;采样时间变小,采样点增多,分散度变小
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2025中文居间服务合同样本
- 2025华彩债券投资基金合同
- 2025股权转让中介合同
- 原单位未提供2025年的劳动合同证明如何办理
- 2025建筑项目施工主要结构劳务分包合同
- 食品生产工艺与质量控制试题及答案
- 2025保险公司赔偿贸易借款合同
- 2025租赁合同模板大全分享
- 2025《租赁合同》范本
- 滇西应用技术大学《塑料制品检测与标准》2023-2024学年第二学期期末试卷
- 湖北省十一校2024-2025学年高三第二次联考数学试卷(解析版)
- 《手工制作》课件-幼儿园挂饰
- 人武专干考试题型及答案
- 2025届高三化学二轮复习 化学反应原理综合 课件
- 2025年北京五湖四海人力资源有限公司招聘笔试参考题库含答案解析
- 常见的酸和碱第2课时酸的化学性质 2024-2025学年九年级化学人教版(2024)下册
- 欢乐购物街-认识人民币(说课稿)-2024-2025学年人教版数学一年级下册
- 2025年中国南方航空股份有限公司招聘笔试参考题库含答案解析
- 疫苗与接种管理制度
- 电子商务案例分析知到智慧树章节测试课后答案2024年秋西安邮电大学
- 《通信用开关电源的元器件降额准则-》
评论
0/150
提交评论