[互助]请帮忙看一下我编的螺旋天线的程序

我自己用Matlab编的矩量法求螺旋天线电流分布的程序，没有错误提示，但是一直在运行，就分了6段，但运行很久都出不来结果，请高手帮忙看一下是哪 .. Pt@%4 :&-h  
'"<6.,Ae  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

程序包括两个部分：luoxuan.m和fun.m   谢谢！

c程序可以运行，只是时间久了些 :a3LS|W  
主要是积分子程序所消耗的时间较长 Di=9m
HC  
-NM0LTF  
同时，楼主画图中的x，y没有注意其包含虚部的影响 noNm^hFL  
!0
Q8iW:  
程序运行结果： 59B&286
1  

图片:2.JPG

我试着把分段增加为72就等了一晚上都没出来，我想是不是死循环了啊？楼上有做过这方面的吗？

“我试着把分段增加为72就等了一晚上都没出来” R"Y?iZed3  
b3^R,6]x&  
不明白楼主说的啥意思？！ =dn1}  
所以没做这方面的研究

无法下载附件啊。

哦 我是把螺旋天线的分段数目增加了，为了使仿真更接近实际，可是就一直有一个警告“maxnium >cBGw'S  
function count exceeded ; singularity likely" ”，不出结果就一直是这个警告，我不知道是怎么回事了，就不知道怎么改了?能告诉我是什么原因吗？楼上的能留个QQ号码？这样交流不方便啊？我的QQ是251190900，邮箱是：poyouxiang@163.com.希望能跟我联系！谢谢！

可以下载啊 是不是你那里网的问题  我刚刚试过   可以的啊

下面这一段有很多地方可以优化！ .BYKdxa  
function y=fun(Sm,n) q 2_N90u  
lambda=1;%波长 %-1BA*J`|  
d=0.000025*lambda;%螺旋线的半径 ~RZJ/%6F  
a=0.15*lambda;%螺旋一圈投影的半径 0i5T] )r  
h=0.05*lambda;%螺距 6oTbn{=UUq  
L=0.1*lambda;%天线的轴向长度 g~ppPAH  
k=2*pi/lambda;%波数 *MnG-\{j  
E=sqrt(h^2+4*pi^2*a^2);%螺旋一圈的长度 \GbHS*\+  
K=L*E/h;%螺旋的总弧长 Cu ['&_@  
C=1/(240*pi*j);%公式中的常系数 lgv-)5|O+H  
F=E/(2*pi);%公式中的常系数 sT[av  
f=@(x)sin(n*pi*(K-x)./K).*(h.^2.*(exp(-j*k*(sqrt(2*a.^2+2*a*d+d.^2+h.^2.*(Sm-x).^2./(4*pi.^2*F.^2)-2*a*(a+d).*cos((Sm-x)./F))))./(sqrt(2*a.^2+2*a*d+d.^2+h.^2.*(Sm-x).^2./(4*pi.^2*F.^2)-2*a*(a+d).*cos((Sm-x)./F))))./(4*pi.^2.*a.^2+h.^2)+(a.^2./F.^2).*cos((Sm-x)./F).*(exp(-j*k*(sqrt(2*a.^2+2*a*d+d.^2+h.^2.*(Sm-x).^2./(4*pi.^2*F.^2)-2*a*(a+d).*cos((Sm-x)./F))))./(sqrt(2*a.^2+2*a*d+d.^2+h.^2.*(Sm-x).^2./(4*pi.^2*F.^2)-2*a*(a+d).*cos((Sm-x)./F))))); 6h0U  
g=quad(f,0,L); =,>TpE  
w=dblquad(@(z,x)((a.^2./F.^2).*sin(n*pi*(K-x)./K).*(exp(-j*k*(sqrt(2*a.^2+2*a*d+d.^2+h.^2.*(Sm-x).^2./(4*pi.^2*F.^2)-2*a*(a+d).*cos((Sm-x)./F)))./(sqrt(2*a.^2+2*a*d+d.^2+h.^2.*(Sm-x).^2./(4*pi.^2*F.^2)-2*a*(a+d).*cos((Sm-x)./F))))).*cos(k*(Sm-z)).*((j*k./(sqrt(2*a.^2+2*a*d+d.^2+h.^2.*(Sm-x).^2./(4*pi.^2*F.^2)-2*a*(a+d).*cos((Sm-x)./F)))+1/(sqrt(2*a.^2+2*a*d+d.^2+h.^2.*(Sm-x).^2./(4*pi.^2*F.^2)-2*a*(a+d).*cos((Sm-x)./F))).^2).*(1-cos((z-x)./F)).*(2*a.^2.*sin((z-x)./F)+h.^2.*(z-x)./(2*pi.^2.*F.^2))-sin((z-x)./F))),0,Sm,0,L); 
')~Y  
y=g-w; [0tfY0  
.%J?T5D  
相同的式子计算一次就可以了，你是计算好几次。增加几个临时变量代替重复的式子再试一下看看。

       我做过代换的，因为是被积函数的表达式，不是一点处的值，所以代换后的结果还是一样，有一个警告，而且一直运行6，7个小时都出不来结果。是不是精度上的问题？不知道哪位有这方面的经验啊？谢谢！