[求助]一个严重待修改的电压闪变程序
本人新手,在做有关电网稳态信号的检测算法实现,下面程序是电压闪变的内容,但是不太会写,请大家帮忙改正。谢谢!!!算法:信号模型: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编辑过]