本人新手，在做有关电网稳态信号的检测算法实现，下面程序是电压闪变的内容，但是不太会写，请大家帮忙改正。谢谢！！！
算法：信号模型：u(t)=U*(1+a(t))*cos(w*t),v(t)=m*cos(w*t),a就是想要得到的调幅波。
流程：先将u平方，然后经过小波解调出a的平方，直流的和高频信号（2w及以上）都不要，然后对得到的a方信号给线性变换和开方

所想要的波形：



N=10000;
u=zeros(1,N);%需不需要？
m=0.1,w1=20*pi,w=100*pi,U=1;%这些参数随便改，只要波形对就好
t=1:0.00015625:5;%,采样信号6400，采多久也不太清楚，波形对就好
v=zeros(1,N);%需不需要？
v=m*cos(20*pi*t);
u(1:size(t,2))=U.*[1+v].*cos(w.*t);%这是版主给的写法
figure(1);
plot(u(1:size(t,2)));
x=u.*u;
figure(2);
plot(x);
l=length(x);
[c,l]=wavedec(x,7,'db4');
a7=wrcoef('v',c,l,'db4',7);%不知是不是这样提取a方的，是不是要滤波
figure(3);
plot(v7);
tf=sqrt(v7);%线性变换先没写
figure(4);
plot(tf);


请大家帮帮忙啊！！！谢谢la!!!

[此贴子已经被作者于2007-5-12 12:13:56编辑过]