|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
各位大侠,在matlab中用randn函数产生的是一系列白噪声,但是我需要对这一串 白噪声进行求积分运算,请问各位大侠如何解决呢?谢谢谢谢
这是程序:
init=205515866;randn('seed',init);
Sig4=rand(1,1000);
SampFreq=500000;
t=0:1/SampFreq:0.014;
Sig1=(t>=0&t<=0.014).*(1+sin(2*pi*15000*t)).*cos(2*pi*60000*t+sin(2*pi*15000*t));
Sig2=(t>=0&t<=0.028).*(1+sin(2*pi*20000*t)).*cos(2*pi*150000*t+sin(2*pi*20000*t));
Sig3=(t>=0.1128&t<=0.0084).*cos(2*pi*150000*t.*(1+sin(2*pi*20000*t))).*cos(2*pi*150000*t+sin(2*pi*20000*t));
Sig4=2*randn(size(Sig1+Sig2+Sig3));
Sig=Sig1+Sig2+Sig3+Sig4;
是求信号Sig的积分,Sig1、Sig2、Sig3都已经求出来啦,程序如下,Sig4的积分怎么求呢?
f1=inline('((1+sin(2*pi*15000*t)).*cos(2*pi*60000*t+sin(2*pi*15000*t))).^2','t');
f2=inline('((1+sin(2*pi*20000*t)).*cos(2*pi*150000*t+sin(2*pi*20000*t))).^2','t');
f3=inline('(cos(2*pi*150000*t.*(1+sin(2*pi*20000*t))).*cos(2*pi*150000*t+sin(2*pi*20000*t))).^2','t');
[a1,kk1]=quad(f1,0,0.014);
[a2,kk2]=quad(f2,0,0.028);
[a3,kk3]=quad(f3,0.0084,0.1128);
a1,a2,a3返回的值就是积分值
|
|