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

我自己用Matlab编的矩量法求螺旋天线电流分布的程序，没有错误提示，但是一直在运行，就分了6段，但运行很久都出不来结果，请高手帮忙看一下是哪 .. qQ[b VD\*  
-NGK@Yk22  

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

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

c程序可以运行，只是时间久了些 ~Aad9yyi  
主要是积分子程序所消耗的时间较长 a\uie$"cr]  
aFiCZHohw  
同时，楼主画图中的x，y没有注意其包含虚部的影响 Hk_y/97OO  
C) QKPT  
程序运行结果： .N4  

图片:2.JPG

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

“我试着把分段增加为72就等了一晚上都没出来” -\ZcOXpMx=  
j hYToMq  
不明白楼主说的啥意思？！ N- e$^pST  
所以没做这方面的研究

无法下载附件啊。

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

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

下面这一段有很多地方可以优化！ dh $bfAb  
function y=fun(Sm,n) D:K4H+ch  
lambda=1;%波长 qcT'nZ:  
d=0.000025*lambda;%螺旋线的半径 )l#E}Uz  
a=0.15*lambda;%螺旋一圈投影的半径 .c$316  
h=0.05*lambda;%螺距 %!1Q P[}K  
L=0.1*lambda;%天线的轴向长度 GDuMY\1  
k=2*pi/lambda;%波数 b3VS\[p  
E=sqrt(h^2+4*pi^2*a^2);%螺旋一圈的长度 95
V@X ^Ee  
K=L*E/h;%螺旋的总弧长 ^53r/V}%  
C=1/(240*pi*j);%公式中的常系数 0c#/hFn  
F=E/(2*pi);%公式中的常系数 aT`%;i^  
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))))); &:7ZQ1  
g=quad(f,0,L); p"U,G -_  
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); >&U,co$>  
y=g-w; 2j&AiD  
]H1I,`=@  
相同的式子计算一次就可以了，你是计算好几次。增加几个临时变量代替重复的式子再试一下看看。

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