求三角级数法模拟轨道谱程序（matlab）

zp3740233

求三角级数法模拟轨道谱程序（matlab） 谢谢各位高手

ch_j1985

LZ能求到程序的几率很小哈……
与其求现成的，还不如自己动手写，有问题大家一起解决！

ChaChing

什麽是轨道谱?

zp3740233

谢谢各位 我写了这方面的程序但得不到想要的结果 我把我的程序发上来

zp3740233

轨道谱就是轨道不平顺功率谱 是轨检车在待测轨道上测量得到的数据 经过傅里叶分析得到的

zp3740233

这个程序还没有考虑误差的

fmin=0; fmax=0.3; t=10; ff=1/t; n=301; ww=(fmax-fmin)/n;
fai=unifrnd(0,2*pi); m1=1;

sww=zeros(1,301);
for k=0:0.0:0.3, sww(k)=sqrt(sw(k)*ww); end
m=zeros(1,100); w=zeros(1,100);
for k=1:100, m(k)=m1+(k-1/2)*ww; end
for x=1:100, w(x)=w(x)+sww(k)*cos(m(k)*x+fai); end

subplot(211); plot(linspace(0,10,length(w)),w,'r-');
y=fft(w,2048); p=40*(y.*conj(y))/2048; ff=0.5*(0:1023)/2048;
subplot(212); loglog(ff,p(1:1024),'k-')

[ 本帖最后由 ChaChing 于 2009-2-23 21:51 编辑 ]

ChaChing

LZ程序错误真的太多! LZ一定初学? 多看看Help及试例题!
就 sww=zeros(1,301); for k=0:0.0:0.3, sww(k)=sqrt(sw(k)*ww); end来说吧!
sww有定义, 那sw试什麽? 还好loop根本没跑, 不然就报错!
试试两者差异
for k=0:0.1:0.3, k, end
for k=0:0:0.3, k, end