[资料共享]基于hallen方程的矩量法分析半波振子

看看看看 ,>rvl P  
是不是真地?????

看看是什么东西！

我要看 >":xnX#  
还应厂个什么啊

HALLEN积分方程

楼主好的，狂顶

怎么能看到代码呢?

我编了一个三角形基函数的可是有问题  请大家帮我看看哪的问题,谢谢！ l:V R8g[  
function y=trfun(k1,m) $
~G,T g  
lambda=1;%波长 !RmVb}m  
a=0.0000275*lambda;%半径 D9oNYF-V  
L=0.25*lambda;%天线半长度 SraZxuPg>  
k=2*pi/lambda;%波数 ':wf%_Iw  
f=@(x)(1-abs(x)/(N-1)).*exp(-j*k*sqrt(a.^2+(k1-x).^2))./sqrt(a.^2+(k1-x).^2);%海伦积分公式 C:J;'[,S  
y=quad(f,(m-1)/(N-1),(m+1)/(N-1));%在取值点k1的积分 /YvXyi>^"%  
J;0;oXwJ<  
Y7}>yC/GY  
clc; Pqomi!1  
clear all [b3!H{b#  
clf; (G"'Fb6d  
lambda=1;%波长 '9V/w[mI  
a=0.0000275*lambda;%半径 V^2_]VFj  
N=51;%取点数 ]^~}/@  
L=0.25*lambda;%天线的半长度 crM5&L9zF  
k=2*pi/lambda;%波数 ^ F]hW  
C=1/(60*j);%公式中的一个常系数 | ;tH?E  
A=zeros(N,N);%生成N*N的零矩阵 m;OvOc,  
for m=1:N -&8( MT*  
    k1=lambda*(m-1)/24;%取点位置 l'+3 6  
    for n=1:(N-1) [MTd<@  
        q(m,n)=trfun(k1,m);%调用trfun函数计算出在k1点处电流方程第m个系数时的积分 5.~Je6K U  
        A(m,n)=q(m,n);%生成的N*(N-1)的矩阵赋值给A QO'Hyf t  
    end 6ZQ |L=Ytp  
end G68KoM  
for i=1:N 8v"rM >[  
    A(i,N)=cos((i-1)*pi/12);%给A矩阵最后一列赋值 M5`v^>  
end 7=6:ZSI  
c=zeros(N,1); t/iI!}  
for i=1:N 2a$.S" ?  
    c(i)=C*sin((i-1)*pi/12); ff#7}9_mh  
end [A-_?#cZ  
B=A\c; _ >OP  
x=-L:1/10000:L; a$My6Qa#  
y=(sum(B))*(1-abs(x)/(N-1)); ")|/\ w,  
plot(x,y)%绘图出电流 ]nTeTW  
xlabel('L/\lambda'),ylabel('电流(A)'); D|LO!,=b  
title('电流分布图(V=1伏)') d/7fJ8y8  
legend('N=51') 9jkz83/+<  
grid on xo"GNFh!  
hold on

xiexiexiexie

不错的  我也想学一下

回一个，让我看看吗