声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1028|回复: 2

[编程技巧] 请问拟合曲线比标准曲线短一截怎么处理?谢谢

[复制链接]
发表于 2009-3-2 20:33 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
以免出现乱码,直接把程序贴在下面了.(matlab)
n=100;wu=0.3;wl=0
ww=(wu-wl)/n
fai=unifrnd(0,2*pi)    %产生(0,2*pi的均匀分布函数)
for i=1:n
    wk(i)=wl+(i-0.5)*ww
end
w=zeros(1,n);sv=zeros(1,n)
w=0:0.0001:0.3
   for i=1:500
sv(i)=9.404e-3*[(w(i)^2+9.701e-2)/(w(i)^4+w(i)^2*3.768e-2+2.666e-5)];   %标准值函数表达式
   end
for i=1:n
    sigma(i)=4*sv(i)*ww   
end
for i=1:n
   ak(i)=normrnd(0,sigma(i))   %产生高斯分布随机数
end
x=zeros(1,n)
for i=1:n
    x(i)=x(i)+ak(i)*sin(wk(i)+fai)   %拟合表达式
end
figure(1)
[p,f]=psd(x,length(x),25)   %求功率谱
loglog(0.01*f,p,'k-')
    hold on
w=zeros(1,100)   %以下产生红线的标准值功率谱
sw=zeros(1,100)
w=0:0.001:0.3
   for i=1:301
sw(i)=9.404e-3*[(w(i)^2+9.701e-2)/(w(i)^4+w(i)^2*3.768e-2+2.666e-5)];
   end
loglog(w,0.01*sw,'r-')

附图中黑色模拟值,红色为标准值,请问有什么方法可以调整好模拟值让两者误差在规定范围内?谢谢大家!
回复
分享到:

使用道具 举报

发表于 2009-3-2 22:11 | 显示全部楼层
 楼主| 发表于 2009-3-3 09:23 | 显示全部楼层

解释一下

就是说,想让黑线短一些,使之和红线在统一的范围内。最好能够让图中黑线和红线尽量重合。谢谢大家的热心帮助。
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-9-22 19:37 , Processed in 0.053972 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表