BBYR Achieve
返回信息流
这是一条镜像帖。来源:北邮人论坛 / matlab / #11457同步于 2014/12/29
Matlab机器人发帖

求助循环平稳信号的仿真??

liumi1125
2014/12/29镜像同步0 回复
我想请教一下,我在做循环平稳信号的高阶统计量的仿真时,用了循环体但是没有体现出循环。而且变量的幅值也是没有体现。clc clear close all N = 128; n = [0:N-1]; Nfft =128; fs = 12000; t = n/fs; A=1;%定义输入信号幅值 fn=3000; fa=12000; fr=8000; fc=10000; fA=20; P=1/20; T=1/120; y=randn(size(t)); e1=mean(y);%%%%均值方差 e1=0; Ai=zeros(1,128*3); e2=mean(Ai); e2=P; i=3; for i=1:3 for j=1:length(t) Ai(j+128*(i-1))=cos(2*pi*fA*t(j)+30)+1+y(j); end end a=0.01*P; mu=0; sigma=a; tao=normrnd(mu,sigma,1,128); B=2;%定义衰减震荡函数幅值 w=zeros(1,128*3); h=zeros(1,128*3); i=3; for i=1:3 for j=1:length(t) w(j+128*(i-1))=t(j)-i*P-tao(j); h=B*exp(-30*w).*sin(2*pi*fn*w);%定义衰减震荡函数 end end figure plot(w,h);%输出t2-h图像 xlabel('w','fontsize',13);%x轴为t2,字号13 ylabel('h','fontsize',13);%y轴为h,字号13 title('h=exp(-w/3).*sin(2*pi*fn*w)');%显示标题 s=Ai.*h; figure(2) plot(w,s);%输出t2-h图像 xlabel('w','fontsize',13);%x轴为t2,字号13 ylabel('h','fontsize',13);%y轴为h,字号13 title('h=exp(-w/3).*sin(2*pi*fn*w)');%显示标题
订阅后,新回复会通过你的通知中心匿名送达。
0 条回复
暂无回复 · 你可以订阅本帖等待新回复。