xiuhuwang 发表于 2009-5-18 10:17

如何画出一个函数的时间历程啊？

这个东西算模拟还是？
我都不知道该怎么问

ChaChing 发表于 2009-5-18 11:14

回复 楼主 xiuhuwang 的帖子

help ezplot
help fplot

xiuhuwang 发表于 2009-5-18 11:34

回复 沙发 ChaChing 的帖子

老师让模拟随机共振，采用 Box-Mueller发放，方法会写，但是不知道怎么画出来

format long
omega=0.1;
A=0.1;
gamma=0.1;
Nfft=2^10;

PT=2.0*pi/omega;
M=Nfft;
dt=PT/M;
N=20;
NM=50;
D=0.2;
x0=1.0;
for i=1:1:N*M
    R=rand(2,1);
    gw=(-4*D*dt*dt*log(R(1,1)))^0.5*cos(2.0*pi*R(2,1));
    t=(i-1)*dt;
    x=x0+(x0*(1+A*cos(omega*t))-x0^3.0)*dt+gw;
    x0=x;
end

friendchj 发表于 2009-5-18 11:41

在x0=x;后加上plot(i,x,'.'); hold on，看看画出来的对吗

xiuhuwang 发表于 2009-5-18 11:47

回复 地板 friendchj 的帖子

画出来了图，我再检测一下，先感谢
ps：好像一直不能终止，我用ctrl + c终止才出来的图
这个是怎么回事？

friendchj 发表于 2009-5-18 12:48

回复 5楼 xiuhuwang 的帖子

能终止，时间稍长一点而已，要做N*M次循环。

xiuhuwang 发表于 2009-5-18 13:27

回复 6楼 friendchj 的帖子

哦，好，我继续试试
页: [1]
查看完整版本: 如何画出一个函数的时间历程啊？