请教一个怪异现象！！StruveH函数
请教一个怪异现象！！StruveH函数对于这个函数StruveH0自己定义.m文件：function y=StruveH0(z)
m=[];
for k=0:20 %%%%%%%%%%%%%%%% k=0,1,2,3....
q=(((-1)^k)/(prod(1:2:2*k+1))^2)*(z^(2*k+1));
m=;
end
y=(2/pi)*sum(m);
下面是通过plot画这个H0函数的.m文件：
symsa1;
Y=[];
for a1=0:17
w=StruveH0(a1);
Y=;
end
figure(1)
plot(Y)
a1代表StruveH0函数里的自变量z
现在可以顺利运行，可就是出了问题
1，虽说K的取值上限是无穷，那我们可以先任意取个值试试，结果随着K的增大，反而StruveH0的图形发生变形，可越大，StruveH0还会出现NAN（非数值）。
2，假定K取30，（我发现最好的取值），a1最多取到24，要不然往后StruveH0就会发生突变，再也不是衰减的了。
3.理论上StruveH0这个函数应该是0-inf衰减的。可真正用公式表示出来后，发现只能在0-20衰减。
请教！！！谢谢。。我已经快抓狂了,不知如何是好！
页:
[1]