matlab程序出错,目的是要求自然风压时程,用谐波叠加法,Davenport谱clear,clct0=clock;N=1024; %采样点数wu=5*2*pi; %截断频率dm=wu/N; %频率步长dt=2*pi/(2*wu); %时间步长k=0.003; %地面粗糙度d=0.01;f=d:d:10;v10=45; %10
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/03 02:24:17
![matlab程序出错,目的是要求自然风压时程,用谐波叠加法,Davenport谱clear,clct0=clock;N=1024; %采样点数wu=5*2*pi; %截断频率dm=wu/N; %频率步长dt=2*pi/(2*wu); %时间步长k=0.003; %地面粗糙度d=0.01;f=d:d:10;v10=45; %10](/uploads/image/z/810854-62-4.jpg?t=matlab%E7%A8%8B%E5%BA%8F%E5%87%BA%E9%94%99%2C%E7%9B%AE%E7%9A%84%E6%98%AF%E8%A6%81%E6%B1%82%E8%87%AA%E7%84%B6%E9%A3%8E%E5%8E%8B%E6%97%B6%E7%A8%8B%2C%E7%94%A8%E8%B0%90%E6%B3%A2%E5%8F%A0%E5%8A%A0%E6%B3%95%2CDavenport%E8%B0%B1clear%2Cclct0%3Dclock%3BN%3D1024%3B+%25%E9%87%87%E6%A0%B7%E7%82%B9%E6%95%B0wu%3D5%2A2%2Api%3B+%25%E6%88%AA%E6%96%AD%E9%A2%91%E7%8E%87dm%3Dwu%2FN%3B+%25%E9%A2%91%E7%8E%87%E6%AD%A5%E9%95%BFdt%3D2%2Api%2F%282%2Awu%29%3B+%25%E6%97%B6%E9%97%B4%E6%AD%A5%E9%95%BFk%3D0.003%3B+%25%E5%9C%B0%E9%9D%A2%E7%B2%97%E7%B3%99%E5%BA%A6d%3D0.01%3Bf%3Dd%3Ad%3A10%3Bv10%3D45%3B+%2510)
matlab程序出错,目的是要求自然风压时程,用谐波叠加法,Davenport谱clear,clct0=clock;N=1024; %采样点数wu=5*2*pi; %截断频率dm=wu/N; %频率步长dt=2*pi/(2*wu); %时间步长k=0.003; %地面粗糙度d=0.01;f=d:d:10;v10=45; %10
matlab程序出错,
目的是要求自然风压时程,用谐波叠加法,Davenport谱
clear,clc
t0=clock;
N=1024; %采样点数
wu=5*2*pi; %截断频率
dm=wu/N; %频率步长
dt=2*pi/(2*wu); %时间步长
k=0.003; %地面粗糙度
d=0.01;
f=d:d:10;
v10=45; %10m高处平均风速
x=1200*f/v10;
s=4*k*v10*v10.*x.^2/f./(1+x).^2.^(4/3); %Davenport谱
z1=4; %取第一点为4m
z2=2; %取第二点为2m
r=0.12; %考虑地面粗糙度的无量纲幂指数
v4=v10*(z1/10)^r; %计算n米高处平均风速
v2=v10*(z2/10)^r;
c=10; %指数衰减系数(取平均值)
v1=zeros(2*N,1);
v2=zeros(2*N,1);
thta1=rand(N,1);
thta2=rand(N,1); %随机数
node=1;
for K=1:node
for j=1:2*N
sum1=0;
sum2=0;
for l=1:N
m1=l*dm-0.5*dm;
m2=l*dm;
x1=1200*m1/(2*pi*v10);
s11=2*pi*4*k.*v4.*v4.*x1*x1./m1./(1+x1*x1).^(4/3);
x2=1200*m2/(2*pi*v10);
s22=2*pi*4*k.*v2.*v2.*x2*x2./m2./(1+x2*x2).^(4/3);
s12=sqrt(s11*s22).*exp(-2*m2*c*abs(z1-z2)./(2*pi*(v4+v2)));
s21=sqrt(s11*s22).*exp(-2*m1*c*abs(z1-z2)./(2*pi*(v4+v2)));
s=[s11 s12;s21 s22];
H=chol(s);
a1=abs(H(1,1));
H1=H';
a21=abs(H1(2,1));
a22=abs(H1(2,2));
b1=cos((m1*dt*(j-1))+2*pi*thta1(1,1));
b2=cos((m2*dt*(j-1))+2*pi*thta2(1,1));
c1=a1*b1;
c21=a21*b1;
c22=a22*b2;
d1=(dm).^0.2*c1;
d2=(dm).^0.2*(c21+c22);
sum1=sum1+d1;
sum2=sum2+d2;
end
sum1=2*sum1;
sum2=2*sum2;
v1(j,K)=sum1;
v2(j,K)=sum2;
end
end
time=etime(clock,t0);
%A、B两点风速时程对比
t=(0:2*N-1)*dt;
subplot(2,1,1);
plot(t,v1,'b-');
xlabel('t(s)');
ylabel('v(t)');
hold on
plot(t,v2,'r-');
xlabel('t(s)');
ylabel('v(t)');
axis([0 1000 -50 50]);
matlab程序出错,目的是要求自然风压时程,用谐波叠加法,Davenport谱clear,clct0=clock;N=1024; %采样点数wu=5*2*pi; %截断频率dm=wu/N; %频率步长dt=2*pi/(2*wu); %时间步长k=0.003; %地面粗糙度d=0.01;f=d:d:10;v10=45; %10
s11与另外三个类型不同
所以数组s=[s11 s12;s21 s22];错误
应该是这里有问题:
v4=v10*(z1/10)^r; %计算n米高处平均风速
v2=v10*(z2/10)^r;
c=10; %指数衰减系数(取平均值)
v1=zeros(2*N,1);
v2=zeros(2*N,1);
改为
v3=v10*(z2/10)^r;
c=10; %指数衰减系数(取平均值)
v1=zeros(2*N,1);
v2=zeros(2*N,1);
.
s22=2*pi*4*k.*v3.*v3.*x2*x2./m2./(1+x2*x2).^(4/3);
就可以了
结果是
希望对你有帮助