[就酱的新征程]数字处理仿真与应用实验二：仿真观测对于不同数量的谐波，其线性叠加与周期方波信号的误差

本文链接：https://blog.csdn.net/qq_43345190/article/details/106560708

@[就酱的新征程]数字处理仿真与应用实验二：仿真观测对于不同数量的谐波，其线性叠加与周期方波信号的误差

这学期数字处理仿真与应用实验的记录
不多说，进入正题吧~

实验要求

仿真观测对于不同数量的谐波，其线性叠加与周期方波信号的误差。

实验结果

在这里插入图片描述 ## 代码

%---------------------%
%------数字信号处理仿真第一次实验------%
%---------------------%
%程序功能：
%----------1.仿真观测对于不同数量的谐波，其线性叠加与周期方波信号的误差。

%编程人员：就酱
%编程更新时间：2020年5月20日

clc;
clear all;
close all;

syms x;
T=10;
t=0:0.001:16;

f=max(square(pi*0.5*t,50),0);%创建方波最大值是1，最小值是0

figure(1)
subplot(411);plot(t,f);grid on;hold on;title('标准方波');
xlabel('t(s)');ylabel('f(t)');axis([0 4*pi -0.5 1.5]);

A0=1/2;
F=0;
for i=1:10  %谐波的阶数
    As=int(2*cos(2*pi*i*x/T)/T,x,0,T/2);%傅里叶系数an
    Bs=int(2*sin(2*pi*i*x/T)/T,x,0,T/2);%傅里叶系数bn
    F=F+As*cos(2*pi*i*t/T)+Bs*sin(2*pi*i*t/T);%求傅里叶级数展开
end
F=F+A0;
subplot(412);plot(t,F);grid on;hold on;title('n=10方波级数分解');
xlabel('t(s)');ylabel('f(t)');axis([0 4*pi -0.5 1.5]);

A0=1/2;
F_1=0;
for i=1:20  %谐波的阶数
    As=int(2*cos(2*pi*i*x/T)/T,x,0,T/2);%傅里叶系数an
    Bs=int(2*sin(2*pi*i*x/T)/T,x,0,T/2);%傅里叶系数bn
    F_1=F_1+As*cos(2*pi*i*t/T)+Bs*sin(2*pi*i*t/T);%求傅里叶级数展开
end
F_1=F_1+A0;
subplot(413);plot(t,F_1);grid on;hold on;title('n=20方波级数分解');
xlabel('t(s)');ylabel('f(t)');axis([0 4*pi -0.5 1.5]);

A0=1/2;
F_1=0;
for i=1:50  %谐波的阶数
    As=int(2*cos(2*pi*i*x/T)/T,x,0,T/2);%傅里叶系数an
    Bs=int(2*sin(2*pi*i*x/T)/T,x,0,T/2);%傅里叶系数bn
    F_1=F_1+As*cos(2*pi*i*t/T)+Bs*sin(2*pi*i*t/T);%求傅里叶级数展开
end
F_1=F_1+A0;
subplot(414);plot(t,F_1);grid on;hold on;title('n=50方波级数分解');
xlabel('t(s)');ylabel('f(t)');axis([0 4*pi -0.5 1.5]);