求助帖，怎么把电流波形里的基波波形分离出来呢？

月野

我用simulink做仿真，有这么一个波形，我现在想要得到的是它的基波波形，该怎么得到呢？ : N0 G6 G8 e" k! B9 e' z$ A) h 我知道powergui里面可以对它进行频谱分析，可以得到基波幅值和相位，前面0.1S基波幅值是54.96,相位是-13.3,后面0.1-0.2S基波幅值是105.7，相位是-19.5，如果我按照前面0.1S基波波形54.96sin(wt-13.3),后面基波幅值是105.7sin(wt-19.5)的话中间连接不上。因为0.1S开始对它FFT分离，选取不同周期的时候相位是有些差别的。 . H. _; }- ~6 N$ P& ~5 i: r MATLAB有没有什么办法可以得到准确得到它的基波波形呢？

redplum

这位仁兄，你的波形好像是晶闸管调压电路，能否分享一下simulink文件，基波需要自己计算

tommy_jiang

前面和后面基波波形本来就不同，分别计算就行了。

月野

回复 2# redplum
# h) B* G+ w/ d

1 ?# w% Q/ r! x" J0 n1 D9 Z2 R' G   是个有源滤波器的模型，波形是整流桥负载变化的情况，我是想看一下负载电流突变的时候，我检验出来的基波跟实际基波的比较，看看它动态响应 yanzheng.rar (16.61 KB, 下载次数: 5)

2013-6-9 10:49 上传

点击文件名下载附件
simulink文件

月野

回复 3# tommy_jiang

  ]/ h+ ?9 n- H* ~2 r! Y9 V  |9 P
  Z  s' a' h) s, W7 O# c   可是这样子画出来中间突变的地方衔接不起来啊，实际情况不应该是连续的么？

图

redplum

我原来曾经做过晶闸管调压的仿真，楼主的波形和3/2模式一致，此种情况基波为为50Hz
% D+ L. p+ P  N+ R3 r0 N，谐波为6K±1次。具体计算公式：
" C" U/ O* Y& M# JD1=fft(VA);
: ~$ M  g* F, f* ?D11=abs(D1);                            % N个点的离散傅立叶变换(DFT)
" n' o5 l% }4 j% q2 qD2=(D11)/N;                             % D2中的第一个数是基波分量的幅值，D2中的第2至第（N/2-1）个数为第1至（（N-1）/2）次谐波分量的幅值的一半                                       
% h6 p# }# Y$ b7 q3 W$ ?- F2 _" yD3=D2;
: s- }2 b& E) S! SD3(1,2:end)=2*D2(1,2:end);
D4=D3(1,1:N/2);                         % D4中的数依次为基波，1，2，3...各次谐波分量的幅值
" s( b) u" ~& n8 ^9 v' b# s8 n%********************************************
D1_i=fft(IA);
D11_i=abs(D1_i);                    % N个点的离散傅立叶变换(DFT)
D2_i=(D11_i)/N;                      % D2中的第一个数是基波分量的幅值，D2中的第2至第（N/2-1）个数为第1至（（N-1）/2）次谐波分量的幅值的一半                                       
# ~1 D4 `6 k# q; E/ X2 CD3_i=D2_i;
- F0 D8 `& ^2 Y. A4 \. _7 j9 ?D3_i(1,2:end)=2*D2_i(1,2:end);
D4_i=D3_i(1,1:N/2);               % D4中的数依次为基波，1，2，3...各次谐波分量的幅值

redplum

这其中VA为Simulink里面波形存成array格式的第二列

wlm_28

路过看看，没有用过simulink