在符号积分中出现的错误,请大虾们指教
%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的警告消息,是那里的问题?
麻烦哪位高人热心解答下，不胜感激！ 原帖由 wenbing1981 于 2007-7-9 20:49 发表 http://www.chinavib.com/forum/images/common/back.gif
%9个空间点 时间间隔0.1，跨度80m，高5.8m
%使用AR模型
clear;
syms f; %频率
k=0.022; %地面粗糙度系数
ds=9; %要模拟的空间点数
t=0.1; %时间间隔
js=4; %模型的阶 ...
没有显式解 出现了这样的警告，是哪里的问题
怎么解决呢 原帖由 wenbing1981 于 2007-7-10 14:53 发表 http://www.chinavib.com/forum/images/common/back.gif
出现了这样的警告，是哪里的问题
怎么解决呢
你看看置顶贴，里面有介绍如何查看出错代码
页:
[1]