求助，電流分布圖只有正半軸，沒有負半軸。
clear all;
clc;
clf;
tic; %計時
lambda=1;
N=31;a=0.0000001;%已知天線和半徑
ii=0;
for h=0.2:0.1:0.9
ii = ii+1;
L=h*lambda;
len=L/N;%將線分成奇數段,注意首末兩端的電流為0
e0=8.854e-012;u0=4*pi*10^(-7);k=2*pi/lambda;
c=3e+008;w=2*pi*c;%光速,角頻率
ata=sqrt(u0/e0);
z(1)=-L/2+len/2;
for n=2:N
z(n)=z(n-1)+len;
end
for m=1:N
for n=1:N
if (m==n)
p(m,n)=log(len/a)/(2*pi)-j*k*len/4/pi;
else
r(m,n)=sqrt((z(m)-z(n))^2+a^2);
p(m,n)=len*exp(-j*k*r(m,n))/(4*pi*r(m,n));
end
end
end
for m=1:N
q(m)=cos(k*z(m));
s(m)=sin(k*z(m));
t(m)=sin(k*abs(z(m)))/(j*2*ata);
end
pp=p(N+1:N^2-N);
pp=reshape(pp,N,N-2);
mat=[pp,q',s'];%構造矩陣
I=mat\t';
II=[0;I(1:N-2);0];%加上兩端零電流
Current=abs(II);
x=linspace(-L/2,L/2,N);
figure(1);
string=['b','g','r','y','c','k','m','r'];
string1=['ko','bo','yo','co','mo','ro','go','bo'];
plot(x,Current,string(ii),'linewidth',1.3);
xlabel('L/\lambda'),ylabel('電流分布');
grid on
hold on
%legend('L=0.1\lambda','L=0.3\lambda','L=0.5\lambda','L=0.7\lambda','L=0.9\lambda','L=1.1\lambda','L=1.3\lambda','L=1.5\lambda')
legend('L=0.1\lambda','L=0.3\lambda','L=0.5\lambda','L=0.7\lambda','L=0.9\lambda','L=1.1\lambda','L=1.3\lambda','L=1.5\lambda')
Zmn=1/I((N+1)/2);%%%%%%V=1v
theta=linspace(0,2*pi,360);
運行結果如圖所示
0.png (38.7 KB, 下載次數: 40)
下載附件
2020-3-1 19:57 上傳
按道理來說應該是這樣子的
1.png (24.53 KB, 下載次數: 42)
下載附件
2020-3-1 19:57 上傳
請問該怎么修改呢？
pp=p(N+1:N^2-N);
pp=reshape(pp,N,N-2);
mat=[pp,q',s'];%構造矩陣
I=mat\t';
II=[0;I(1:N-2);0];%加上兩端零電流
Current=abs(II);
是從這一塊下手嗎？
matlab2019a Win10 1903
|
QQ好友和群
QQ空間
騰訊微博
騰訊朋友
收藏
淘帖
頂
踩
