 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 572 | 回復(fù): 1 | ||
wjdwin新蟲 (初入文壇)
|
[求助]
二維TM波，PML邊界，matlab程序。但是仿真結(jié)果有問題，求高手解決
|
|
clc; clear; ft=1.5; fs=20; NX=100;%X方向網(wǎng)格數(shù) NY=150;%Y方向網(wǎng)格數(shù) NA=20;%吸收層厚度 bxrhx=1:100; bxirhx=2:150; byrhy=2:100; byirhy=1:150; bxrez=1:100; byrez=1:150; hx=zeros(100,151); hy=zeros(101,150); ezx=zeros(100,150); ezy=zeros(100,150); sigmamy=hx; sigmamx=hy; sigmax=ezx; sigmay=ezy; sigmamax=0.2; sigma=sigmamax*(NA:-1:1)/NA; sigmam=sigmamax*(2*NA-1:-2:1)/2/NA; sigmamx([1:NA,end:-1:end-NA+1],  =repmat([sigmam,sigmam]',1,150);sigmax([1:NA,end:-1:end-NA+1], =repmat([sigma,sigma]',1,150);sigmamy(:,[1:NA,end:-1:end-NA+1])=repmat([sigmam,sigmam],100,1); sigmay(:,[1:NA,end:-1:end-NA+1])=repmat([sigma,sigma],100,1); f1hx=exp(-sigmamy(bxrhx,bxirhx)); f1hy=exp(-sigmamx(byrhy,byirhy)); f2hx=zeros(size(hx)); f2hy=zeros(size(hy)); f2hx(sigmamy == 0)=1/ft; f2hy(sigmamx==0)=1/ft; f2hx(sigmamy~=0)=(1-exp(-sigmamy(sigmamy~=0)))./sigmamy(sigmamy~=0)/ft; f2hy(sigmamx~=0)=(1-exp(-sigmamx(sigmamx~=0)))./sigmamx(sigmamx~=0)/ft; f2hx=f2hx(bxrhx,bxirhx); f2hy=f2hy(byrhy,byirhy); f1ezx=exp(-sigmax); f1ezy=exp(-sigmay); f2ezx=zeros(size(ezx)); f2ezy=zeros(size(ezy)); f2ezx(sigmax==0)=1/ft; f2ezy(sigmay==0)=1/ft; f2ezx(sigmax~=0)=(1-exp(-sigmax(sigmax~=0)))./sigmax(sigmax~=0)/ft; f2ezy(sigmay~=0)=(1-exp(-sigmay(sigmay~=0)))./sigmay(sigmay~=0)/ft; centerx=floor((1+NX)/2); centery=floor((1+NY)/2); for nn=1:150 ezx(centerx,centery)=ezx(centerx,centery)+sin(2*pi*nn/ft/fs); ezy(centerx,centery)=ezy(centerx,centery)+sin(2*pi*nn/ft/fs); hx(bxrhx,bxirhx)=f1hx.*hx(bxrhx,bxirhx)-f2hx.*(ezx(bxrhx,bxirhx)-ezx(bxrhx,bxirhx-1)+ezy(bxrhx,bxirhx)-ezy(bxrhx,bxirhx-1)); hy(byrhy,byirhy)=f1hy.*hx(byrhy,byirhy)+f2hy.*(ezx(byrhy,byirhy)-ezx(byrhy-1,byirhy)+ezy(byrhy,byirhy)-ezy(byrhy-1,byirhy)); ezx=f1ezx.*ezx+f2ezx.*(hy(bxrez+1, -hy(bxrez,);ezy=f1ezy.*ezy-f2ezy.*(hx(:,byrez+1)-hx(:,byrez)); end 程序可以運(yùn)行。但是自己畫出的等位線圖不對(duì)，contour（ezx+ezy）,contour(hx),contour(hy),都不對(duì)，求解決。 |
回復(fù)此樓» 猜你喜歡
285求調(diào)劑
已經(jīng)有6人回復(fù)
-
085701環(huán)境工程求調(diào)劑
已經(jīng)有6人回復(fù)
-
291求調(diào)劑
已經(jīng)有6人回復(fù)
-
315分求調(diào)劑
已經(jīng)有6人回復(fù)
-
277跪求調(diào)劑
已經(jīng)有8人回復(fù)
-
295求調(diào)劑
已經(jīng)有5人回復(fù)
-
295材料工程專碩求調(diào)劑
已經(jīng)有3人回復(fù)
-
083000學(xué)碩274求調(diào)劑
已經(jīng)有4人回復(fù)
-
275求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩，求調(diào)劑
已經(jīng)有8人回復(fù)
wjdwin
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 41.5
- 帖子: 7
- 在線: 9.6小時(shí)
- 蟲號(hào): 1138733
- 注冊(cè): 2010-11-04
- 專業(yè): 建筑物理
2樓2015-11-12 22:12:34
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 wjdwin 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿北京化工大學(xué)材料與化工 264分各科過A區(qū)國家線 +10 | 哈哈157349 2026-03-21 | 10/500 |
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
|
[考研] 0703化學(xué)338求調(diào)劑！ +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
|
[考研] 343求調(diào)劑 +4 | 贈(zèng)我一本書 2026-03-23 | 4/200 |
|
|
|
[考研] 329求調(diào)劑 +5 | 1() 2026-03-22 | 5/250 |
|
|
|
[考研] 材料科學(xué)與工程 317求調(diào)劑 +4 | JKSOIID 2026-03-26 | 4/200 |
|
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap－ 2026-03-22 | 10/500 |
|
|
|
[考研] 299求調(diào)劑 +7 | 某某某某位 2026-03-21 | 8/400 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
|
[考研] 0854電子信息求調(diào)劑 324 +4 | Promise-jyl 2026-03-23 | 4/200 |
|
|
|
[考研] 食品專碩 一志愿雙一流 328 +3 | xiaom99 2026-03-21 | 4/200 |
|
|
|
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +4 | 17501029541 2026-03-23 | 6/300 |
|
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
|
[考研] 306求調(diào)劑 +5 | 來好運(yùn)來來來 2026-03-22 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
|
[考研] 求助 +5 | 夢(mèng)里的無言 2026-03-21 | 6/300 |
|