[求助]这个积分需要加点什么条件能够得到合理的结果

[风雨浪涛]

有哪位朋友能帮忙看看这个积分，需要加点什么条件能够得到合理的结果。
clear
clc

syms k v10 n t
xn=1200*n/v10;
s1=4*k*v10^2*xn^2/n/(1+xn^2)^(4/3)*cos(2*pi*n*t);
F=int(s1,n,0,inf);
先谢谢了！

[laosam280]

至少要知道t是怎么变化的，即使是求参变量的积分。
另外，n的变化是不是0，1，2，。。。，这样的话，就不可以直接用积分命令，而且把n当成积分符号了

[laosam280]

另外，我这里有个能求解简单奇异积分的Gauss-Lagarre公式：
function I = IntGaussLager(f,n,AK,XK)
if(n<6 && nargin == 2)
AK = 0;
XK = 0;
else
I=sum(AK.*subs(sym(f),findsym(sym(f)),XK));
end

switch n
case 2,
I=0.853553*subs(sym(f),findsym(sym(f)),-0.585786)+...
0.146447*subs(sym(f),findsym(sym(f)),3.414214);

case 3,
I=0.711093*subs(sym(f),findsym(sym(f)),0.415575)+...
0.278518*subs(sym(f),findsym(sym(f)),2.294280)+...
0.0103893*subs(sym(f),findsym(sym(f)),6.289945);

case 4,
I=0.603154*subs(sym(f),findsym(sym(f)),0.322548)+...
0.357419*subs(sym(f),findsym(sym(f)),1.745761)+...
0.0388879*subs(sym(f),findsym(sym(f)),4.536620)+...
0.000539295*subs(sym(f),findsym(sym(f)),9.395071);

case 5,
I=0.521756*subs(sym(f),findsym(sym(f)),0.263560)+...
0.398667*subs(sym(f),findsym(sym(f)),1.413403)+...
0.0759424*subs(sym(f),findsym(sym(f)),3.596426)+...
0.00361176*subs(sym(f),findsym(sym(f)),7.085810)+...
0.0000233700*subs(sym(f),findsym(sym(f)),12.640801);
end

[风雨浪涛]

感谢你在百忙之中回复我的帖子。我的源程序中的k,v10都是已知量，只不过用字母来表示。频率n是从0变到无穷，t没有具体的变化范围。我只是想通过积分把F表示成t的函数，然后再利用这个函数画图，此时t我取的是从-200s到200s。请你再给细看一下好吧，谢谢了。

[laosam280]

你说源程序中的k,v10都是已知量，那都是常数吗？如果是常数的话，我可以再计算这个积分，尝试其他的数值积分公式。我的邮箱[email protected]，请你把有关的资料发给我，谢谢