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

矩量法分析半波振子天线的采用的是三角形基函数 Z*w(
{k7]  
function y=trfun(k1,m) VtzZ1/JE  
lambda=1;%波长 Pi=FnS  
a=0.0000275*lambda;%半径 |)To 0Z  
L=0.25*lambda;%天线半长度 $/ "+t.ir3  
k=2*pi/lambda;%波数 6R^F^<<  
f=@(x)(1-abs(x)/(N-1)).*exp(-j*k*sqrt(a.^2+(k1-x).^2))./sqrt(a.^2+(k1-x).^2);%海伦积分公式 A@|Z^T:  
y=quad(f,(m-1)/(N-1),(m+1)/(N-1));%在取值点k1的积分 1J1Jp|j.  
ZiY2N*,VO  
{J1rjrPo  
clc; TJRp/BP  
clear all
wuY-f4  
clf; <-N eusx%  
lambda=1;%波长 5P #._
Em  
a=0.0000275*lambda;%半径 uq9mq"  
N=51;%取点数 <<S4l~"o  
L=0.25*lambda;%天线的半长度 !{V`N|0  
k=2*pi/lambda;%波数 7Fpa%N/WL  
C=1/(60*j);%公式中的一个常系数 **__&Xp1  
A=zeros(N,N);%生成N*N的零矩阵 "JmbYb#Z  
for m=1:N XsJ`x  
    k1=lambda*(m-1)/24;%取点位置 [V_
mF  
    for n=1:(N-1) :T$}@& -  
        q(m,n)=trfun(k1,m);%调用trfun函数计算出在k1点处电流方程第m个系数时的积分 Bn8&~  
        A(m,n)=q(m,n);%生成的N*(N-1)的矩阵赋值给A XelFGTE  
    end (9( xJ)  
end 9$v\D3<Z  
for i=1:N {(-923|,  
    A(i,N)=cos((i-1)*pi/12);%给A矩阵最后一列赋值 0y<9JvN$9  
end ApSseBhh  
c=zeros(N,1); n1|]ji[c  
for i=1:N hoBFC1  
    c(i)=C*sin((i-1)*pi/12); 4N=,9  
end An(gHi;1$  
B=A\c; i/ o  
x=-L:1/10000:L; )?&mCI*  
y=(sum(B))*(1-abs(x)/(N-1)); rMWvW(@@D  
plot(x,y)%绘图出电流 *g}==o`  
xlabel('L/\l .. OO/>}? ob  
zx"EAF{  

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