内容发布更新时间 : 2024/5/9 12:19:59星期一 下面是文章的全部内容请认真阅读。
实验程序： 实验结果： a=[1 4 4]; b=[1 3]; t=0:0.01:10; x=exp(-t); lsim(b,a,x,t); 三、实验总结： （1）用matlab求解连续时间信号时域卷积，有了更直观的认识 （2）知道了线性时不变系统时域分析方法，对连续时间信号的冲激响应，阶跃响应，零状态相应有更直观的认识
实验题目 一、实验目的： 实验三 周期信号频域分析 （1）掌握周期信号傅立叶级数分解与合成的计算公式 （2）掌握利用 MATLAB 实现周期信号傅立叶级数分解与综合方法 （3）理解并掌握周期信号频谱特点 二、实验内容： function [A_sym,B_sym]=CTFShchsym % 采用符号计算求一个周期内连续时间函数 f的三角级数展开系数,再用这些 % 展开系数合成连续时间函数 f.傅立叶级数 实验结果： % 函数的输入输出都是数值量 % Nf=6 谐波的阶数 % Nn 输出数据的准确位数 % A_sym 第 1 元素是直流项，其后元素依次是 1,2,3...次谐波 cos 项展开系数 % B_sym 第 2,3,4,...元素依次是 1,2,3...次谐波 sin 项展开系数 % tao=1 tao/T=0.2 syms t n k x T=5;tao=0.2*T;a=0.5; if nargin<4; Nf=10; end if nargin<5;Nn=32; end x=time_fun_x(t); A0=int(x,t,-a,T-a)/T; %求出三角函数展开 系数 A0 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)); %获取串数组 A0 所对应的 ASC2 码数值数组 for k=1:Nf A_sym(k+1)=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