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

你先测试一下你的函数是否正确，如果正确，代换后的速度应该比原来的快3~4倍

for m=1:N+1 YS9)%F=X  
    Sm=K*(2*m-1)/2*N+1;%取点位置---》这个地方是不是有问题？是否应该是：Sm=K*(2*m-1)/(2*N+1) $ ga,$G  
    for n=1:N 2Sy:wt  
        q(m,n)=fun(Sm,n);%调用fun函数计算出在Sm点处电流方程第n个系数时的积分 Pmuk !V}f  
        A(m,n)=q(m,n);%生成的(N+1)*N的矩阵赋值给A $43CNnf3N  
     end y#ON=8l  
end `ER#S_}  
E
xi#@-  
('_S1?y  
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))))); <+ <o X"I  
g=quad(f,0,L); "H=fWz5z  
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); ]~m=b
`o  
y=g-w; dZ,7q_r,~  
下面直接用数值代替：1、F=E/(2*pi)；2、4*pi.^2*F、等等，你仔细检查一下，我要是你的老师，我得狠狠批评你。

你好！我试着做代换，几种方法都不行，syms, @, regexprep, 以及单引号'   ',    这些都与积分相矛盾了，还有就是int 和quad这两种积分应该用哪一中呢？我把做代换的子函数传上来了，希望你能帮我看一下，给一些建议，我的仿真做的比较少，所以还请您赐教！谢谢！

下面是我想的四种代换方法

我感觉这个地方有问题： yi1V\8DC  
R 9Yk9v  
Sm=K*(2*m-1)/2*N+1;%取点位置，我感觉是不是应该写成： Sm=K*(2*m-1)/(2*N+1);%取点位置,如果是，这个地方也可以优化：就是在添加一个临时变量：K1 = K / ( 2 * N + 1 )，然后：Sm = K1 * ( 2 * m - 1 );

学习一下   

你应该反馈一下问题是否解决，这样大家就都明白了

引用第16楼kerbcurb于2008-11-03 11:04发表的 : 8xs[{?|:  
你应该反馈一下问题是否解决，这样大家就都明白了

3gaijVN  
}a,j1r_Hl&  
说的很对啊，大家问题得到解决最好都反馈一下