风速时程曲线仿真(请大家关注)

bainhome

function fengsu
clear all;
m=0.01;
v10=20;                                              %10米高处风速为20m/s
k=0.005;
t=0;
i=1;
s=38.0818;
V=s;
thta=rand(100,1);                                    %随机数
for t=1:100
    for i=2:100
        x=1200*m*i/v10;
        s11=4*k*v10^2*x.^2./(m*i)./(1+x.^2).^(4/3);  %Davenport谱
        s12=s11.*exp(-(m*i)./v10);
        S=[s11 s12;s12 s11];
        H=chol(S);                                   %Cholesky分解
        a=det(H(1,1));
        b=cos(((m*i)*t)+2*pi*thta(i,1));
        c=a.*b;
        v(t)=2*(m).^(1/2)*c;                         %风荷载模拟
        V=V+v(t);
    end
end
plot(t,V,'ro')                                            %显示风荷载
xlabel('t(s)');
ylabel('v(t)');
axis([95 105 -500 1000])
稍稍改动一下，你的程序没有保存前面的结果，所以只有最后t＝100一个点

semeng

怎么就一个点?

zhouyinhua

同问，同问，同问