MATLAB周期信号的频谱分析.ppt

分析 例题 奇谐函数的频谱图 代码 % 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