返回信息流%产生功率谱
Uz=23.6;z=10;
for m=1:1024
n(m)=0.001*m;f(m)=1200*n(m)/Uz;
Pn(m)=4.*0.03.*Uz^2*f(m).^2./(1+f(m).^2).^(4/3);%风速谱,如计算竖%向Panofsky谱将其替换为相应函数,这是davenport自功率谱密度函数
Uf=1;%摩擦速度
Su(m)=Pn(m)/n(m);
end
V10=18;h=1000;p=10;t0=0;ti=0.01;SUMT=1024;fig=25;
for n=0:p
for i=1:h
for j=i:h
H=inline('Su')%s-功率谱密度
dx=x(i)-x(j)
dz=z(i)-z(j)
R(i,j,n+1)=quadl(Su)
R(j,i,n+1)=R(i,j,n+1)
end
end
end
求教,关于这段程序出现的问题及解决办法,,谢谢。。。。。。。。
这是一条镜像帖。来源:北邮人论坛 / matlab / #6814同步于 2010/4/25
Matlab机器人发帖
求教,关于inline和quadl
lullf
2010/4/25镜像同步0 回复
订阅后,新回复会通过你的通知中心匿名送达。
0 条回复
暂无回复 · 你可以订阅本帖等待新回复。