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

看看看看 ,jn?s^X6Dj  
是不是真地?????

看看是什么东西！

我要看 L9T u
>4  
还应厂个什么啊

HALLEN积分方程

楼主好的，狂顶

怎么能看到代码呢?

我编了一个三角形基函数的可是有问题  请大家帮我看看哪的问题,谢谢！ _1<'"u#6w  
function y=trfun(k1,m) ceZ8}Sh  
lambda=1;%波长 UVc<C 1q  
a=0.0000275*lambda;%半径 vo ;F;  
L=0.25*lambda;%天线半长度 oL'1Gm@X?  
k=2*pi/lambda;%波数 GrR0RwnH)?  
f=@(x)(1-abs(x)/(N-1)).*exp(-j*k*sqrt(a.^2+(k1-x).^2))./sqrt(a.^2+(k1-x).^2);%海伦积分公式 d-;9L56{P  
y=quad(f,(m-1)/(N-1),(m+1)/(N-1));%在取值点k1的积分 hNnX-^J<o  
(\t_Hs::a  
xDJ+BQ<1A  
clc; ej%;%`C-  
clear all NOr <,  
clf; !?o$-+a|  
lambda=1;%波长 qmA2bw]  
a=0.0000275*lambda;%半径 eQFb$C]R}y  
N=51;%取点数 ZQ~myqx,+L  
L=0.25*lambda;%天线的半长度 @f{_=~+  
k=2*pi/lambda;%波数 R pI<]1
  
C=1/(60*j);%公式中的一个常系数 BwJ^_:(p~  
A=zeros(N,N);%生成N*N的零矩阵 {Mr~%y4  
for m=1:N c~UAr k S  
    k1=lambda*(m-1)/24;%取点位置 knOnUU  
    for n=1:(N-1) RO{@RhnV  
        q(m,n)=trfun(k1,m);%调用trfun函数计算出在k1点处电流方程第m个系数时的积分 +Y)#yGUn  
        A(m,n)=q(m,n);%生成的N*(N-1)的矩阵赋值给A ydND$@; Z  
    end HNy/ -  
end =6L*!JP<  
for i=1:N Y/,$Y]%g  
    A(i,N)=cos((i-1)*pi/12);%给A矩阵最后一列赋值 b"M`@';+  
end ml.;wB|  
c=zeros(N,1); pEVgJ/>  
for i=1:N D!}K)T1~R  
    c(i)=C*sin((i-1)*pi/12); ]hE+$sKd  
end .S!>9X,  
B=A\c; Vp$wHB&  
x=-L:1/10000:L; xN3 [Kp  
y=(sum(B))*(1-abs(x)/(N-1)); a3UPbl3^  
plot(x,y)%绘图出电流 g[s\~MF@s  
xlabel('L/\lambda'),ylabel('电流(A)'); bugFl>  
title('电流分布图(V=1伏)') qX{X4b$  
legend('N=51') 3`k[!!   
grid on 3XwU6M$5g  
hold on

xiexiexiexie

不错的  我也想学一下

回一个，让我看看吗