信号与系统上机实验第二次matlab实现.docxVIP

信号与系统第二次上机 第三次试验 实验一： 仿照例程，实现下述周期信号的傅立叶级数分解与合成：I (a)首先，推导出求解绻，乞，?的公式，计算出前10次系数; (b)利用MATLAB求解％,①,?的值，其中比，?求解前10次系数，并给 出利用这些系数合成的信号波形。 function [A_sym,B_sym]=CTFShchsym % N仁10 谐波的阶数 % Nn 输出数据的准确位数 % A_sym 第1元素是直流项，其后元素依次是123…次谐波cos项展开系数 % B_sym 第2,3,4,...元素依次是1,2,3...次谐波sin项展开系数 % tao=l tao/T=0.25 syms t n k x T=4;tao=0.25*T;a=0.5; 讦 nargin4;Nf= 10;end if nargin5;Nn=32;end x=time_fun_x(t); A0=int(x,t,-a,T-a)rr； %求出三角函数展开系数AO As=int(2*x*cos(2*pi*n*t/T)/T,t,-a,T-a); %求出三角函数展开系数 As Bs=int(2*x*sin(2*pi*n*t/T)/T,t,-a,T-a); %求出三角函数展开系数 Bs A_sym( 1 )=double(vpa(A0,Nn)); %获取串数组AO所对应的ASC2码数值数组 for k=l:Nf A_sym(k+l)=double(vpa(subs(As,n,k),Nn)); %获取串数组 A 所对应的 ASC2 码数值数组 B_sym(k+1 )=double(vpa(subs(Bs,n,k),Nn)); %获取串数组 B 所对应的 ASC2 码数值数组 end if nargout==0 c=A_sym;disp(c) %输出c为三角级数展开系数:第1元素是直流项，其后元素依次 是1,2,3…次谐波cos项展开系数 d=B_sym;disp(d)谐波sin项展开系数%输出d为三角级数展开系数：第2,3,4,... d=B_sym;disp(d) 谐波sin项展开系数 %基波% 2 %基波 % 2次谐波 % 3次谐波 % 4次谐波 fl =c( 1 )+c(2).*cos(2*pi* 1 *t/4)+d(2).*sin(2*pi* 1 *t/4);; f2=c(3).*cos(2*pi*2*t/4)+d(3).*sin(2*pi*2*t/4);; f3=c(4).*cos(2*pi*3*t/4)+d(4).*sin(2*pi*3*t/4); f4=c(5).*cos(2*pi*4*t/4)+d(5).*sin(2*pi*4*t/4);; f5=c(6).*cos(2*pi*5*t/4)+d(6).*sin(2*pi*5*t/4); f6=c(7).*cos(2*pi*6*t/4)+d(7).*sin(2*pi*6*t/4); f7=c(8).*cos(2*pi*7*t/4)+d(8).*sin(2*pi*7*t/4); f8=c(9).*cos(2*pi*8*t/4)+d(9).*sin(2*pi*8*t/4); f9=c( 10).*cos(2*pi*9*t/4)+d( 10).*sin(2*pi*9*t/4); fl 0=c( 11 ).*cos(2*pi* 10*t/4)+d( 11). *sin(2*pi* 10*t/4); fl2=fll+f3+f4+f5+f6; fl3=fl2+f7+f8+f9+flO; subplot(2,2,l) plot(t,fl),hold on y=time_fun_e(t) %调用连续时间函数■周期矩形脉冲 plot(t,y,”) mi*周期矩形波的形成一基波+直流j axis([-4,4.5,-0.5,1.5]) subplot(2,2,2) plot(t,fll),hold on y=time_fun_e(t) plot(t,y；r:*) title。周期矩形波的形成一基波+2次谐波+3次谐波J axis([-4,4.5,-0.5,1.5]) subplot(2,2,3) plot(t,fl2),hold on y=time_fun_e(t) plot(t,y,r:*) titleCS波+2次谐波+3次谐波….+6次谐波) axis(卜 4,4.5,?0.5,1.5]) subplot(2,2,4) plot(t,fl3),hold on y=time_fun_e(t) plot(t,y；r:) titled基波+2次谐波+3次谐波+4次谐波+10次谐波J axis(卜 4,4.5,?0.5,1.5]) end % function x=time_fun_x(t) %该函数是CTFShc