- 1、本文档共69页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
MATLAB周期信号的频谱分析
分析 例题 奇谐函数的频谱图 代码 % exe_zheng_3_7_c.m clear all;close all;clc; syms t x T n; integrand1=(-2*(t+1)/T)*exp(-i*n*2*pi/T*t); integrand2=(2*t/T)*exp(-i*n*2*pi/T*t); Xn=int(integrand1,t,-T/2,0)/T+int(integrand2,t,0,T/2)/T; Xn_N=subs(Xn,T,5); Xn_N=[subs(Xn_N,n,-10:-1) 0 subs(Xn_N,n,1:10)]; n=-10:10; h=stem(n,abs(Xn_N)); set(h,linewidth,2); 频谱图 奇谐函数含有基波(n=1)及奇次(n=3,5,7…)谐波分量 step=0.01; x=0:step:10*pi; y=sinc(x/pi); subplot(2,1,1); h=plot(x,y);set(h,LineWidth,2); Si_wave=cumsum(y)*step; subplot(2,1,2); h=plot(x,Si_wave);set(h,LineWidth,2) 正弦函数的泰勒级数展开 x=0:0.01:2*pi; y=x; K=100 for ii=2:1:K y_add=power(-1,ii+1)*power(x,ii*2-1)/prod(1:(ii*2-1)); y=y+y_add; plot(x,y,k,x,sin(x),y:); end 画Si函数波形 syms x y si_func=int(sin(x)/x,0,y); 7.6 周期信号的功率 (b) 解: (c) n=-9:9; Fn=sinc(n/2); Fn(10)=0; P=sum(power(abs(Fn),2)) 解: (d) n=-99:99; Fn=sinc(n/2); Fn(100)=0; P_n=2*power(abs(Fn(100:199)),2); P_n(1)=0; P_evaluation=cumsum(P_n); 例题 单位冲激序列的频谱 分析:狄氏条件是傅里叶级数存在的充分条件。根据冲激信号的定义和特性,其积分有确定值,傅里叶级数存在。即 满足离散性,谐波性,不满足收敛性,频带无限宽。 余弦形式的频谱图 MATLAB实现 % impulse_cos.m % this program is used to illustrate impuse series can be expanded by cos T=1; w1=2*pi/T; t=-1.2*T:0.01:1.2*T; f=ones(1,length(t))/T; clf; plot(t,f,’r’); hold on color=[r,g,b,c,y]; for n=1:5 fn=cos(n*w1*t)/T*2; f=f+fn; plot(t,fn,color(mod(n,5)+1)); end h=plot(t,f,k); set(h,linewidth,2); 结果显示(N=5) 结果显示(N=21) 总结 pause mesh waterfall cumsum sum clf GUIDE 回调函数编写:get, set, axes num2str strcat mod 7.7 用傅里叶级数综合连续时间信号 例题 % exe3_11_a.m clear all syms t T X1 X_1 X3 X_3 x1=X_1*exp(-i*2*pi/T*t)+X1*exp(i*2*pi/T*t)+X_3*exp(-3*i*2*pi/T*t)+X3*exp(3*i*2*pi/T*t); T=1;X1=5;X_1=5;X3=2;X_3=2; x2=subs(x1);x3=simple(x2); ezplot(x3,[0,2*T]); 更新参数 T=2; X3=1; X_3=1; x3=subs(x1); figure; ezplot(x3,[0,2*T]); 符号函数 x1 = X_1*exp(-2*i*pi/T*t) +X1*exp(2*i*pi/T*t) +X_3*exp(-6*i*pi/T*t) +X3*exp(6*i*pi/T*t) x2 = 5*exp(-2*i*pi*t)+5*exp(2*i*pi*t) +2*exp(-6*i*pi*t)+2*exp(6*i*pi*t) x3 = 5*exp(-i*pi*t)+5*exp(i*pi*t) +exp(-3*i
文档评论(0)