1955阅读
0回复

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

离线yyqxb

UID ：8061

注册：2008-01-30
登录：2008-05-27
发帖：7
等级：旁观者

0楼发表于: 2008-05-20 21:06:41

— 本帖被 casey 从矩量法MOM 复制到本区(2008-05-20) —

矩量法分析半波振子天线的采用的是三角形基函数
function y=trfun(k1,m)
lambda=1;%波长
a=0.0000275*lambda;%半径
L=0.25*lambda;%天线半长度
k=2*pi/lambda;%波数
f=@(x)(1-abs(x)/(N-1)).*exp(-j*k*sqrt(a.^2+(k1-x).^2))./sqrt(a.^2+(k1-x).^2);%海伦积分公式
y=quad(f,(m-1)/(N-1),(m+1)/(N-1));%在取值点k1的积分

clc;
clear all
clf;
lambda=1;%波长
a=0.0000275*lambda;%半径
N=51;%取点数
L=0.25*lambda;%天线的半长度
k=2*pi/lambda;%波数
C=1/(60*j);%公式中的一个常系数
A=zeros(N,N);%生成N*N的零矩阵
for m=1:N
k1=lambda*(m-1)/24;%取点位置
for n=1:(N-1)
q(m,n)=trfun(k1,m);%调用trfun函数计算出在k1点处电流方程第m个系数时的积分
A(m,n)=q(m,n);%生成的N*(N-1)的矩阵赋值给A
end
end
for i=1:N
A(i,N)=cos((i-1)*pi/12);%给A矩阵最后一列赋值
end
c=zeros(N,1);
for i=1:N
c(i)=C*sin((i-1)*pi/12);
end
B=A\c;
x=-L:1/10000:L;
y=(sum(B))*(1-abs(x)/(N-1));
plot(x,y)%绘图出电流
xlabel('L/\l ..

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

发帖回复