[求助]我编的程序不知道哪的问题

矩量法分析半波振子天线的采用的是三角形基函数 E1W:hGI  
function y=trfun(k1,m) kp`0erJqw  
lambda=1;%波长 'F5)ACA%  
a=0.0000275*lambda;%半径 hTw}X.<4  
L=0.25*lambda;%天线半长度 NG9vml  
k=2*pi/lambda;%波数 ]kS7n@8  
f=@(x)(1-abs(x)/(N-1)).*exp(-j*k*sqrt(a.^2+(k1-x).^2))./sqrt(a.^2+(k1-x).^2);%海伦积分公式 ;$;rD0i|  
y=quad(f,(m-1)/(N-1),(m+1)/(N-1));%在取值点k1的积分 k/#>S*Ne  
u^=@DO
'  
Gvn: c/m;  
clc; 6V;:+"
BkJ  
clear all }m_t$aaUc1  
clf; 5Y-2 #  
lambda=1;%波长 UC|JAZL  
a=0.0000275*lambda;%半径 QymD-A"P  
N=51;%取点数 ORH93`  
L=0.25*lambda;%天线的半长度 uaxB -PZ  
k=2*pi/lambda;%波数 ragSy8M  
C=1/(60*j);%公式中的一个常系数 p >aw  
A=zeros(N,N);%生成N*N的零矩阵 \!wh[qEQ\  
for m=1:N NI\H \#bJ  
    k1=lambda*(m-1)/24;%取点位置 F^rl$#pCS  
    for n=1:(N-1) xOIg|2^8  
        q(m,n)=trfun(k1,m);%调用trfun函数计算出在k1点处电流方程第m个系数时的积分 b1=pO]3u  
        A(m,n)=q(m,n);%生成的N*(N-1)的矩阵赋值给A Wk[)+\WQ?  
    end q\H7&w  
end EUZq$@uWL  
for i=1:N !$r9C/k  
    A(i,N)=cos((i-1)*pi/12);%给A矩阵最后一列赋值 v3^t/[e~:  
end t w!.%_1^  
c=zeros(N,1); k4R4YI"jV  
for i=1:N 0N VI+Z$  
    c(i)=C*sin((i-1)*pi/12); *Sb2w*c>  
end /{we;Ut=g  
B=A\c; Nza; O[  
x=-L:1/10000:L; -+P7:4/  
y=(sum(B))*(1-abs(x)/(N-1)); e\.HWV]I  
plot(x,y)%绘图出电流 }vD;DSz:  
xlabel('L/\l .. @?/\c:cp  
*?N<S$m  

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