TM11

僅供參考嗷！！

clc;

clear;

count = input(‘迭代次數:’);%迭代次數

w = 2/(1+sin(pi/(count-1)));%收斂因子

h = pi/7;a = pi/11;%矩形波導尺寸

xa = input(‘單邊節點數:’);%單邊節點

tol = 10e-3;

%初始化

u=zeros(xa,xa);

%內點賦值

for i = 2:(xa - 1)

for j = 2:(xa - 1)

u(i , j) = 1;

end

end

%邊界賦值

u(1,:) = 0; u(xa,:)= 0;

u(:,1) = 0; u(:,xa)= 0;

%kc初始化

kc = 0.1;

%迭代初始化

k = 0;

sum1 = 0; sum2 = 0; sum = 0;

%差分(超松弛)

while k < count

err = 0;

temp = u(i,j);

for i = 2 : xa - 1

for j = 2: xa - 1

u(i,j)=u(i,j)+w*((((u(i+1,j)+u(i,j+1)+u(i-1,j)+u(i,j-1)))/(4-(kc*h)^2))-u(i,j));

err = err + abs(temp - u(i,j));

end

end

k = k + 1;

%對kc進行迭代

for i = 2: xa-1

for j = 2 : xa-1

sum1 = sum1 + u(i,j)(u(i+1,j)+u(i,j+1)+u(i-1,j)+u(i,j-1)-4u(i,j));

sum2 = sum2 + u(i,j)^2;

end

end

kc = (sqrt(-sum1/sum2))/h;

%保證精確度

if(err

break

end

end

surf(u);

title(‘TM11’);

總結

以上是生活随笔為你收集整理的怎么用matlab画TM11,矩形波导TM11模matlab仿真的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。