周期信号的频谱
更新时间:2023-12-22 03:13:01 阅读量:2 教育文库 文档下载
实验二 周期信号的频谱
三 实验的参考程序
%傅立叶级数的部分和,最高谐波次数为3,21,41和81的波形比较 sy2_1.m %锯形波形 clear all
n_max=[3 21 41 81]; N=length(n_max); t=-1.1:.001:1.1; omega_0=2*pi; for k=1:N n=[];
n=[1:n_max(k)];
b_n=2./(pi*n).*(-1).^(1+n); x=b_n*sin(omega_0*n'*t);
subplot(N,1,k),plot(t,x,'linewidth',2); axis([-1.1 1.1 -1.5 1.5]);
line([-1.1 1.1],[0 0],'color','r'); line([0 0],[-1.5 1.5],'color','r');
bt=strcat('最高谐波次数=',num2str(n_max(k))); title(bt); end
%三角波形 clear all
n_max=[3 21 41 81]; A=1;
N=length(n_max); t=-1.1:.002:1.1; omega_0=2*pi; figure; for k=1:N n=[];
n=[1:n_max(k)]; b=(pi*n).^2;
a=(sin(pi*n/2)).^2; a_n=4./b.*a;
x=A*a_n*cos(omega_0*n'*t);
subplot(N,1,k),plot(t,x,'linewidth',2); axis([-1.1 1.1 -1.5 1.5]);
line([-1.1 1.1],[0 0],'color','r'); line([0 0],[-1.5 1.5],'color','r');
bt=strcat('最高谐波次数=',num2str(n_max(k))); title(bt); end
首先定义一个抽样函数(Sa.m)放在实验程序相同的目录下, Sa.m
%抽样函数(连续或离散) %高度为1,
%调用 y=Sa(t) 函数幅度为1,第一个过零点为pi Function f=Sa(t) F=sinc(t./pi);
% 画频谱,T变化,tau不变 sy2_2.m
clear all
tau=0.5; t=[2 5 10];
tau_T=tau./t; %占空比 omega_0=2*pi./t n0=0;n1=50; N=length(t); for k=1:N
n=[n0:omega_0(k):n1]; F_n=tau_T(k)*Sa(tau*n/2) Fn_max=max(F_n); Fn_min=min(F_n);
subplot(N,1,k),stem(n,F_n,'.');
axis([n0 n1 Fn_min-0.01 Fn_max+0.01]);
line([n0 n1],[0 0],'color','r'); %画直线,表示横轴,线为红色
title(strcat('幅度频谱:脉冲宽度=0.5,周期=',num2str(t(k))));% 在图上写标题 end
% sy2_3.m clear all
tau=[0.2 0.5 1];t=5;
tau_T=tau./t; % 占空比 omega_0=2*pi/t
n0=0;n1=50;N=length(tau); n=[n0:omega_0:n1]; for k=1:N
F_n=tau_T(k)*Sa(tau(k)*n/2); Fn_max=max(F_n); Fn_min=min(F_n);
subplot(N,1,k),stem(n,F_n,'.'); % 画幅度频谱 axis([n0 n1 Fn_min-0.01 Fn_max+0.01]);
line([n0 n1],[0 0 ],'color','r'); % 画直线,表示横轴,线为红色
title(strcat('幅度频谱:周期=5,脉冲宽度=',num2str(tau(k)))); % 在图上写标题 end
%周期锯形信号的分解 clear all
t=0:0.001:2*pi; w=1;
N=[1;3;5;7;9;11;13;15;17;19]; k=1; while k<=10 n=N(k,:); b=4./(pi*n); x=b*sin(w*n*t);
subplot(5,2,k),plot(t,x); xlabel('t'),ylabel('x(t)'); axis([0,2*pi,-2,2]),
legend(strcat(num2str(N(k)),'次谐波')) k=k+1; end
2x(t) x(t)20-2 02t43次谐波 0-2 022t41次谐波6 x(t)6 7次谐波20-2 02t4x(t)0-2 022t45次谐波6 x(t)6 11次谐波20-2 02t4x(t)0-2 022t49次谐波6 x(t)6 15次谐波20-2 02t4x(t)0-2 022t413次谐波6 x(t)6 19次谐波20-2 02t4x(t)0-2 02t417次谐波66
%周期锯形信号的合成 n_max=[7 21 41 99]; N=length(n_max); t=-2.1:.002:2.1; omega_0=pi; for k=1:N n=[];
n=[1:2:n_max(k)]; a0=1/2;
a_n1=2./(pi*n); a_n2=sin(pi*n/2); a_n=a_n1.*a_n2;
x=a0+a_n*cos(omega_0*n'*t);
subplot(N,1,k),plot(t,x,'linewidth',2); axis([-2.1 2.1 0 1.5]);
line([-2.1 2.1],[0 0],'color','r'); line([0 0],[-1.5 1.5],'color','r');
bt=strcat('最高次谐波数=',num2str(n_max(k))); title(bt); end
最高谐波次数=71.510.501.510.501.510.501.510.50-2-1.5-1-0.500.5最高谐波次数=2111.52-2-1.5-1-0.500.5最高谐波次数=4111.52-2-1.5-1-0.500.5最高谐波次数=9911.52-2-1.5-1-0.500.511.52
正在阅读:
周期信号的频谱12-22
西政安徽校友会通讯录01-19
编译原理课程设计报告05-17
和萨克斯做伴的我作文700字07-02
抖音文案:拍摄短视频脚本文案(建议收藏)03-27
模拟炒股大赛02-17
青贮发酵剂的原理09-29
电子技术复习题(1)05-09
报关单填制规范04-10
心事作文400字04-01
- exercise2
- 铅锌矿详查地质设计 - 图文
- 厨余垃圾、餐厨垃圾堆肥系统设计方案
- 陈明珠开题报告
- 化工原理精选例题
- 政府形象宣传册营销案例
- 小学一至三年级语文阅读专项练习题
- 2014.民诉 期末考试 复习题
- 巅峰智业 - 做好顶层设计对建设城市的重要意义
- (三起)冀教版三年级英语上册Unit4 Lesson24练习题及答案
- 2017年实心轮胎现状及发展趋势分析(目录)
- 基于GIS的农用地定级技术研究定稿
- 2017-2022年中国医疗保健市场调查与市场前景预测报告(目录) - 图文
- 作业
- OFDM技术仿真(MATLAB代码) - 图文
- Android工程师笔试题及答案
- 生命密码联合密码
- 空间地上权若干法律问题探究
- 江苏学业水平测试《机械基础》模拟试题
- 选课走班实施方案
- 频谱
- 周期
- 信号
- 景区运营部工作总结
- 用友ERP-U860应用中常见问题解析 - 图文
- 中国传统文化-古为今用,辩证发展
- 泉州市洛江区对外贸易经济合作局
- 广东高考数学理试题下载A卷word版
- 工程竣工资料第四册(给排水)
- (新课标)2019版高考物理一轮复习 主题五 能量和动量 综合技能演练5-1-6.doc
- 同桌的你观后感800字
- 2018国家公务员考试时政热点:双11别让利兴冲昏了理性
- 二期导向竣工验收监理总结(修改)2014.1.20
- 2016年国家公务员行测选词填空每日一练(2)
- 山东省曲阜市书院街道办事处圣林中学八年级地理上册 3.2 节约与保护水资源教案(新版)商务星球版
- 初中语文知识点《文学常识及鉴赏》《近代文学》课时作业(含答案考点及解析)
- 社区XX第六次全国人口普查工作汇报(1)
- 第一讲 绪论
- 20种塑胶特性
- 上海市桃浦中学16届和谐校园材料 - 图文
- 小吃2
- 浙江省人力资源和社会保障厅、浙江省档案局关于印发浙江省档案
- 圆锥曲线三种弦长问题