梯形波的傅里叶级数分解

已于 2023-09-06 07:53:50 修改
阅读量9.5k 收藏 33
点赞数 26
分类专栏： 信号 文章标签： 硬件工程
于 2022-09-04 01:18:04 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_36160506/article/details/126684862
版权

文章目录

1.前言

想必大家搜索看到这篇文章的时候，大概已经是踏入信号与系统的大门，不满足于书上的简单的矩形波三角波，想探探梯形波，今天它来了~

前几天在CSDN 硬件工程师炼成之路 大佬 的一篇文章《不同的反射系数，信号的是怎样的呢？》写的很好，里面大概是讲了源端发射到终端，用信号反射的观点，去计算源端和终端的电压波形。

在这里插入图片描述

用的是梯形波：

在这里插入图片描述
这种边沿可变+脉宽可控的梯形波，是非常有意义的。

因为实际的TX端，边沿不可能无限陡峭，脉宽也要做成可控。这个波形比单纯的矩形波更常用。

傅里叶级数真的优雅，一般书本上的习题也都是三角波、矩形波，积分起来相对容易。脑子里突然一热，想把边沿可变+脉宽可控的梯形波，用傅里叶级数分解掉。
在这里插入图片描述

有点飘了，数学底子不行，不会用对称性，积分也不熟练了，没什么技巧，就是死做。

2.动手开搞

周五晚上就酝酿了，徘徊退缩了很久想一想还是积一下，网上仅有的资源也几乎只有一个答案，而且边沿不灵活，数学表达式也不是我想要的。

列表达式：

在这里插入图片描述
然后开始漫长的积分，很长很长，很细很细，像是高中的笨小孩，死做椭圆的题，设未知数，代入方程越解越长：

在这里插入图片描述

过程非常的痛苦，实不相瞒持续了3个小时，输完所有公式，从18.30持续纠错到近了22.30，洋洋洒洒7-8页的重复劳动。

*所有详细的推导放在网盘里了

在这里插入图片描述

3.导入matlab 进行plot

matlab是个高级计算器了，可以帮你验证下积分的结果正确与否。
把word里的公式全部打入 matlab文件，看看波图

3.1 A=1,T=100,m=10,h=30,p=20

（1）振幅A=1，周期T=100，上升时间m=10，正脉宽=30，下降时间p=20；
在这里插入图片描述

3.2 A=1,T=100,m=5,h=40,p=5

（2）振幅A=1，周期T=100，上升时间m=5，正脉宽=40，下降时间p=5；
在这里插入图片描述

3.3 A=2,T=150,m=20,h=50,p=1

（3）振幅A=2，周期T=150，上升时间m=20，正脉宽=50，下降时间p=1；

在这里插入图片描述

3.4 奇偶对称性

在这里插入图片描述
可以看到An是偶函数，Bn是奇函数，A0是f(t)的直流分量。
波形非常的完美，不过无法想象，为什么要在那些点跳变，把这三个值相加就正好是一个梯形波。

在这里插入图片描述

4.代码


clear all;
close all;

%parameter set


n=1:1:10;

A=2;T=150;m=20;h=50;p=20;

t=0:0.01:2*T;






str=strcat('A=',num2str(A),',T=',num2str(T),',m=', ...
num2str(m),',h=',num2str(h),',p=',num2str(p) )

%(2./T)
%n.*omg
omg=2.*pi./T;
%(n.*omg)


%	A./(n.*omg)

g1 =(1./(n.*omg));


X1= (A./(n.*omg)).*sin(n.*omg.*m) + ...
(A./m).*(	1./	(  (n.*omg).*(n.*omg)	)	).*cos(n.*omg.*m)- ...
(A./m).*(	1./	(  (n.*omg).*(n.*omg)	)	);


X2=(A./(n.*omg)).*sin(n.*omg.*m+n.*omg.*h) - ...
(A./(n.*omg)).*sin(n.*omg.*m);



% v1 (A./p).*(1./(n.*omg))
v1 = (A./p).*(1./(n.*omg));

X3_1= -v1.*(m+h+p) .* sin(n.*omg.*m+n.*omg.*h+n.*omg.*p) + ...
v1 .*(m+h).*sin(n.*omg.*m+n.*omg.*h) + ...
-v1.*g1.*cos(n.*omg.*m+n.*omg.*h+n.*omg.*p) + ...
v1.*g1.*cos(n.*omg.*m+n.*omg.*h);





v2 =(A.*m)./(p.*n.*omg);
X3_2 = v2 .* sin(n.*omg.*m+n.*omg.*h+n.*omg.*p) - ...
v2 .* sin(n.*omg.*m+n.*omg.*h);


v3 = (A)./ (n.*omg);
X3_3=v3 .*  sin(n.*omg.*m+n.*omg.*h+n.*omg.*p) - ...
v3 .* sin(n.*omg.*m+n.*omg.*h);

v4 =(A.*h)./(p.*n.*omg);
X3_4 = v4 .* sin(n.*omg.*m+n.*omg.*h+n.*omg.*p) - ...
v4 .* sin(n.*omg.*m+n.*omg.*h);

X3=X3_1 + X3_2 + X3_3 + X3_4;

an=(2./T).*X1+(2./T).*X2+(2./T).*X3;





% 



K1=(-A.*g1) .*cos(n.*omg.*m)+ (A./m).*g1.*g1.*sin(n.*omg.*m);
K2= -A.*g1 .* cos(n.*omg.*m + n.*omg.*h) + A.*g1.* cos(n.*omg.*m);
g2=(A./p);
K3_1=g2.*g1.*(m+h+p).*cos(n.*omg.*m + n.*omg.*h+n.*omg.*p) - ...
g2.*g1.*(m+h).*cos(n.*omg.*m + n.*omg.*h) - ...
g2.*g1.*g1.*sin(n.*omg.*m + n.*omg.*h+n.*omg.*p) + ...
g2.*g1.*g1.*sin(n.*omg.*m + n.*omg.*h);
g3=(A.*m)./(p.*n.*omg);

K3_2 = -g3.*cos(n.*omg.*m + n.*omg.*h+n.*omg.*p) + ...
g3.*cos(n.*omg.*m + n.*omg.*h) ;

K3_3 = -A.*g1.*cos(n.*omg.*m + n.*omg.*h+n.*omg.*p) + ...
A.*g1.*cos(n.*omg.*m + n.*omg.*h);


g4 = (A.*h)./(p.*n.*omg);

K3_4 = -g4.* cos(n.*omg.*m + n.*omg.*h+n.*omg.*p) + ...
+g4.* cos(n.*omg.*m + n.*omg.*h);

K3=K3_1 + K3_2 + K3_3 + K3_4;

bn=(2./T).*K1+(2./T).*K2+(2./T).*K3;




ft=zeros(1,length(t));


A0= A.*m./(2.*T) + ...
A.*h./T + ...
-(A./(2.*T.*p))       .*       (     (m+h+p).* (m+h+p)  -(m+h).*(m+h)      ) +...
(1./T) .* (A.*m./p + A +A.*h./p ) .*p ;


for j =1:1:n(end)
	ft=ft + an(j) .* cos(j.*omg.*t)  + bn(j) .* sin(j.*omg.*t);
end

figure(1)
plot(t, A0 + ft)

m_max= max(A0 + ft);
m_min= min(A0 + ft);

axis([0,t(end) m_min-1.1 m_max+1.1])

legend(str) 


xlabel('time axis')
ylabel('amplitude axis')
set(gca,'Fontname','Times New Roman','FontSize',15);


if 0
    figure(2)
    AN=zeros(1,length(t));
    BN=zeros(1,length(t));
    for j =1:1:n(end)
        AN=AN + an(j) .* cos(j.*omg.*t);
    end

    for j =1:1:n(end)
        BN=BN + bn(j) .* sin(j.*omg.*t);
    end

    plot(t, AN,':r','Linewidth',1);
    hold on;
    plot(t, BN,':b','Linewidth',1);
    hold on;
    plot(t, A0.*ones(1,length(t)),':k','Linewidth',1);
    hold on;
    plot(t, A0+BN+AN,'-k','Linewidth',1);
    hold on;

    legend('An','Bn','A0','f(t)');
    xlabel('time axis')
    ylabel('amplitude axis')
    set(gca,'Fontname','Times New Roman','FontSize',15);
end



% + bn(j) .* sin(j.*omg.*t)

完整的文档放在网盘自取：
链接：https://pan.baidu.com/s/1kHA5rh2I8Fx54K3No3O79Q?pwd=AAAA
提取码：AAAA



欢迎关注我的公众号：沈土豪的书屋

傅里叶变换
qflook的博客
06-03 5449
傅里叶变换（Fourier transform）是一种线性积分变换，用于信号在时域（或空域）和频域之间的变换，在物理学和工程学中有许多应用。一般情况下，若“傅里叶变换”一词不加任何限定语，则指的是“连续傅里叶变换”（连续函数的傅里叶变换）。傅里叶变换的前提条件是符合狄里赫利条件，即要绝对可积，周期内有限个不连续点等。
矩形梯形波的傅里叶变换及截断影响
04-22
文档及matlab 程序，，不同谐波次数逼近原信号，
几种常见波形的傅里叶级数展开式
热门推荐
吃蛋糕的小学生的博客
01-03 5万+
公式 给定一个周期为的函数，那么它可以表示为无穷级数： 其中傅里叶系数为： 几种常见波形的傅里叶级数展开式 梯形波（奇函数） 傅里叶展开为： 脉冲波（偶函数） 傅里叶展开为： 方波（奇函数） 傅里叶展开为： 三角波（奇函数） 傅里叶展开为： 锯齿波（非奇非偶函数） 傅里叶级数展开式为： ...
方波等各种信号的傅里叶分解
08-03
傅里叶分解 方波信号，单边指数信号等展开为傅里叶级数 ……
傅里叶变换分类
u011555996的博客
02-12 8431
傅里叶变换 傅里叶变换（Fourier transform）是一种线性的积分变换，从时间转换为频率的变化 1 1. 连续傅里叶变换 这是将频率域的函数F(ω)表示为时间域的函数f（t）的积分形式 连续傅里叶变换的逆变换 (inverse Fourier transform)为： 一般可称函数f（t）为原函数，而称函数F(ω)为傅里叶变换的像函数，原函数和像函数构成一个傅里叶变换对（tran...
【信号分解的艺术】：傅立叶级数与变换的比较分析
本文综述了信号分解领域中傅立叶分析的理论基础及其实践应用，阐述了傅立叶级数、连续时间傅立叶变换和离散时间傅立叶变换的概念与数学定义，探讨了快速傅立叶变换（FFT）算法的原理和优化。文章进一步对比了傅立叶...
傅立叶变换..pdf
最新发布
10-04
傅立叶算法的基本思想源于傅里叶级数。该算法假设输入信号为一周期性信号，即输入信号中除基频分量外，只包含恒定的直流分量和各 种整 次 谐波分量。傅立叶算法可以将输入信号分解为多个谐波分量，每个谐波分量对应...
MATLAB傅里叶级数计算器的GUI实现
傅里叶级数是一个数学工具，它允许我们将周期函数或信号分解为一系列正弦和余弦...通过傅里叶级数，可以将复杂的信号分解为简单的正弦波和余弦波的组合，从而更容易地理解信号的特性、进行滤波处理或信号合成等操作。
【引言与基础知识】傅里叶变换的历史背景：傅里叶级数与傅里叶积分
它的历史可追溯至18世纪，当时的科学家们致力于研究如何将复杂的波形分解为简单的正弦波组合。这一理论基础是由法国数学家让-巴蒂斯特·约瑟夫·傅里叶提出，他的工作奠定了傅里叶级数和傅里叶变换
信号与系统中关于傅里叶的几种变换（学习笔记）
米安r的博客
05-01 4704
信号与系统中关于傅里叶的几种变换简介时域连续连续周期连续非周期时域离散周期非周期证明 简介 这篇文章是参考奥本海姆的Signals and Systems和B站上西电的信号与系统这门课，字母的选择可能会与视频里的有些许的差别。这里的讨论主要是对几个基本公式的简易推导。先有一个简单的结论吧，就是 时域离散，频域周期 时域周期，频域离散 然后这个结论会在接下来的几个式子中得到证明。 时域连续 傅里叶发...
MATLAB中方波、三角波、梯形波的生成：方波、三角波、梯形波的仿真-matlab开发
05-29
Qadeer Sir 方波、三角波和梯形波的生成/模拟
matlab绘制梯形并标注参数
04-18
本文的目的: 绘制梯形并标注参数 参数为：斜率，上边沿的宽度以及各个转折点的标注 可以自行调参 Amp = 2; Fre = 2; Ste = 56; Breite = (1/(4*Fre))*((Ste/4/Fre-1)/(Ste/4/Fre))*2; x1 = (0+Ste/4/Fre)/Ste; x2 = (1+Ste/(4*Fre))/Ste; x3 = x2 + Breite; x4 = (0+(3*Ste)/(4*Fre))/Ste; figure(1) x = 0:0.00001:x4; y = trapmf(x,[x1 x2 x3 x4]); plot(x,y,'k',[x2 x2],[0 1.1],'--k',[0 x4],[0 0],'--k'); text(x2+Breite/3,0.95,['\uparrow ' 'Breite=' (num2str(Breite,2))]); text(x2+Breite/40,0.5,['\leftarrowSteigung=' (num2str(Ste,2))]); xlabel(['P = [',(num2str(x1,3)),', ',(num2str(x2,3)),', ',(num2str(x3,3)),', ',(num2str(x4,3)),']']) xlim([0 x4]) ylim([-0.05 1.05])
trapezoid:产生梯形波-matlab开发
06-01
这会创建一个梯形波，类似于内置的平方函数，但没有占空比选项
方波、锯齿波、三角波、梯形波、正弦波（使用电脑的两款软件联机实现）
06-04
org 0000h loot: jnb p1.0,aa jnb p1.1,bb jnb p1.2,cc jnb p1.3,dd jnb p1.4,ee sjmp loot aa: mov r7,#10 net: mov p2,#0 acall delay mov p2,#0ffh acall delay djnz r7,net sjmp loot bb: mov r6,#10 lll: mov A,#0 ll: mov p2,A inc a mov r2,#10 djnz r2,$ cjne a,#0ffh,ll djnz r6,lll sjmp loot cc: mov r5,#10 looop:mov a,#0 loop: mov p2,a inc a mov r2,#10 djnz r2,$ cjne a,#0ffh,loop gg: mov p2,a dec a mov r2,#10 djnz r2,$ jnz gg djnz r5,looop sjmp loot dd: start:mov r1,#10 start1:mov a,#0 xx: mov p2,a mov r2,#10 djnz r2,$ inc a cjne a,#0ffh,xx acall delay2 yy: mov p2,a dec a mov r2,#10 djnz r2,$ jnz yy mov p2,a acall delay2 djnz r1,start1 sjmp loot ee: MOV R6,#10 LOOP2: MOV DPTR,#table LOOP3:MOV A,#0 MOVC A,@A+DPTR MOV P2,A ACALL DELAY4 INC DPTR mov r4,dph CJNE r4,#01h,loop3 mov r4,dpl cjne r4,#0f7h,loop3 DJNZ R6,LOOP2 ljmp loot DELAY4:mov r2,#10 djnz r2,$ ret TABLE: DB 128, 130, 132, 135, 137, 139, 141, 144, 146, 148 DB 150, 152, 155, 157, 159, 161, 163, 165, 168, 170 DB 172, 174, 176, 178, 180, 182, 184, 186, 188, 190 DB 192, 194, 196, 198, 200, 201, 203, 205, 207, 209 DB 210, 212, 214, 215, 217, 219, 220, 222, 223, 225 DB 226, 227, 229, 230, 232, 233, 234, 235, 237, 238 DB 239, 240, 241, 242, 243, 244, 245, 246, 247, 247 DB 248, 249, 250, 250, 251, 252, 252, 253, 253, 254 DB 254, 254, 255, 255, 255, 255, 255, 255, 255, 255 DB 255, 255, 255, 255, 255, 255, 255, 255, 255, 254 DB 254, 254, 253, 253, 252, 252, 251, 250, 250, 249 DB 248, 247, 247, 246, 245, 244, 243, 242, 241, 240 DB 239, 238, 237, 235, 234, 233, 232, 230, 229, 227 DB 226, 225, 223, 222, 220, 219, 217, 215, 214, 212 DB 210, 209, 207, 205, 203, 201, 200, 198, 196, 194 DB 192, 190, 188, 186, 184, 182, 180, 178, 176, 174 DB 172, 170, 168, 165, 163, 161, 159, 157, 155, 152 DB 150, 148, 146, 144, 141, 139, 137, 135, 132, 130 DB 128, 126, 124, 121, 119, 117, 115, 112, 110, 108 DB 106, 104, 101, 99, 97, 95, 93, 91, 88, 86 DB 84, 82, 80, 78, 76, 74, 72, 70, 68, 66 DB 64, 62, 60, 58, 56, 55, 53, 51, 49, 47 DB 46, 44, 42, 41, 39, 37, 36, 34, 33, 31 DB 30, 29, 27, 26, 24, 23, 22, 21, 19, 18 DB 17, 16, 15, 14, 13, 12, 11, 10, 9, 9 DB 8, 7, 6, 6, 5, 4, 4, 3, 3, 2 DB 2, 2, 1, 1, 1, 0, 0, 0, 0, 0 DB 0, 0, 0, 0, 0, 0, 1, 1, 1, 2 DB 2, 2, 3, 3, 4, 4, 5, 6, 6, 7 DB 8, 9, 9, 10, 11, 12, 13, 14, 15, 16 DB 17, 18, 19, 21, 22, 23, 24, 26, 27, 29 DB 30, 31, 33, 34, 36, 37, 39, 41, 42, 44 DB 46, 47, 49, 51, 53, 55, 56, 58, 60, 62 DB 64, 66, 68, 70, 72, 74, 76, 78, 80, 82 DB 84, 86, 88, 91, 93, 95, 97, 99, 101, 104 DB 106, 108, 110, 112, 115, 117, 119, 121, 124, 126 DB 128 delay:mov r2,#200 djnz r2,$ ret delay2: mov r3,#60 zz: mov r2,#0 djnz r2,$ djnz r3,zz ret end
信号处理、傅里叶变换 、小波变换的入门资料
bestlinjialyin的专栏
09-28 3855
概念 http://en.wikipedia.org/wiki/Signal_processing Signal processing is an area of systems engineering, electrical engineering and applied mathematics that deals with operations on or analysis of an
傅里叶变换性质和常见信号的傅里叶变换
集成电路设计那些事儿
02-13 9172
微信公众号：xiaoshi_IC 开年第一篇，单位里没有网，只能用12年前的老电脑发点文章，老电脑还行，虽然装个软件cpu会跑到100%，但是不卡，哈哈。 老电脑见证了摩尔定律，我也没想到最终会做IC，FPGA相关的。 废话不多说了，直接上图，熟练的掌握这些性质，还是很有用的。 ...
五、卷积与傅立叶变换
qq_27531383的博客
05-19 2万+
一、卷积1、一维的卷积连续：在泛函分析中，卷积是通过两个函数f(x)f(x)和g(x)g(x)生成第三个函数的一种算子，它代表的意义是：两个函数中的一个(我取g(x)g(x)，可以任意取)函数，把g(x)g(x)经过翻转平移,然后与f(x)f(x)的相乘，得到的一个新的函数，对这个函数积分，也就是对这个新的函数求它所围成的曲边梯形的面积。 def： 设f(t),g(t)f(t),g(t)是Rn\
MATLAB之傅里叶变换，快速傅里叶变换FFT
Heart_Sea的博客
04-14 4350
文章目录傅里叶变换及傅里叶逆变换定义窗函数/矩形脉冲信号的傅里叶变换基于MATLAB的快速傅里叶变换FFT 傅里叶变换及傅里叶逆变换定义 能从时域的非周期连续信号转化到频域非周期连续信号。 窗函数/矩形脉冲信号的傅里叶变换 结论： 随着脉冲宽度τ的减小，主叶变得更宽，而且更多的能量被移到更高的频率。则反之。 当信号脉冲在时间上扩展时，它的变换在频率上压缩。则反之。 基于MATLAB的快速...
三相电压型方波逆变电路，逆变的波形电压计算。1/3U，2/3U来源于两个电阻并联再与第三电阻串联分压。两个电容串联接在直流电源两端也分压，分压一半。6阶波，4阶波概念以及进行傅里叶变换。基波谐波
二进制模----细微行动改变影响世界
06-01 1万+
https://zhuanlan.zhihu.com/p/274965441?utm_id=0
hive left semi join not in
09-07
在Hive中，left semi join和not in是两种不同的操作。left semi join是一种连接操作，它只传递左表的join key给Map阶段，并返回只包含左表的结果。而not in是一种条件操作，用于筛选不在指定列表中的值。 在left semi join中，只有左表中的匹配项会被返回，而不会返回右表中的任何记录。这意味着在left semi join中，如果右表中有重复记录，左表中的匹配项也只会返回一次。这种机制使得left semi join在处理有重复值的情况下更加高效。 对于not in操作，它使用一个子查询来排除在指定列表中的值。如果某个值在子查询的结果中找不到匹配项，则该值会被返回。not in可以用于在查询中筛选出不满足某些条件的数据。 综上所述，left semi join和not in在功能上有一些差异。left semi join用于连接两个表并返回左表的匹配项，而not in用于筛选不在指定列表中的值。<span class="em">1</span><span class="em">2</span><span class="em">3</span><span class="em">4</span>
评论 7
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

沈土豪的书屋

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值