|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
% 9个空间点 时间间隔0.1,跨度80m,高5.8m
%使用AR模型
clear;
syms f; %频率
k=0.022; %地面粗糙度系数
ds=9; %要模拟的空间点数
t=0.1; %时间间隔
js=4; %模型的阶数
sd=25; %高度十米处的平均风速
kuadu=80; %结构跨度
g=5.8; %结构高度
for p=1:ds %每个空间点的平均风速
v(p)=25;
end
%生成每个空间点的坐标
x=linspace(0,kuadu,ds);
y=4*g/kuadu.*x.*(1-x/kuadu);
z=linspace(0,0,ds);
for p=1:ds
for q=1:p
h(p,q)=exp(-2*f*(256*(x(p)-x(q))^2+64*(y(p)-y(q))^2+100*(z(p)-z(q))^2)^0.5/(v(p)+v(q))); %风相干函数
end
end
%%
for p=1:ds
for q=1:p
if p==q
s=4*k*sd^2*(1200*f/sd)^2/(f*(1+(1200*f/sd)^2)^(4/3))*cos(2*pi*f*0*t); %风自相关功率谱
else
s=4*k*sd^2*(1200*f/sd)^2/(f*(1+(1200*f/sd)^2)^(4/3))*h(p,q)*cos(2*pi*f*0*t); %风互相关功率谱
end
r0(p,q)=int(s,f,0,inf); %积分区间
end
end
for p=1:ds
for q=1:p
if p==q
s=4*k*sd^2*(1200*f/sd)^2/(f*(1+(1200*f/sd)^2)^(4/3)); %自相关功率谱
else
s=4*k*sd^2*(1200*f/sd)^2/(f*(1+(1200*f/sd)^2)^(4/3))*h(p,q); %互相关功率谱
end
r1(p,q)=int(s*cos(2*pi*f*t),f,0,inf); %积分区间
end
end
出现了Explicit integral could not be found的警告消息,是那里的问题?
麻烦哪位高人热心解答下,不胜感激! |
|