应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 用fortran编写的GaAs MESFET求电流不收敛
61/1返回列表
查看: 703  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

xingfendou

新虫 (初入文坛)

[求助] 用fortran编写的GaAs MESFET求电流不收敛

用fortran语言编写的GaAs MESFET求电流不收敛，但是我认为电荷分配，求电场，散射率，散射机制的选择等都对，就是求的电流不收敛，还有正有负，我已经找了一个月的错误都没找到，求高手指点！
if (iter.gt.itermax-itervel+1) then
         curtot=0.
         do 1130 i=1,iycur
            curtot=curtot+(0.9*charge1(i)+0.1*charge2(i))*
     #                    (0.9*velx1(i)+0.1*velx2(i))
           curtot=curtot+charge2(i)*velx2(i)
1130     continue
         curtot=curtot*constcur
         write(41,*)curtot
      endif
回复此楼

» 猜你喜欢
1楼 2013-05-02 22:48:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

夕阳西下

金虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
do 1130 i=1,iycur
没有判断收敛的条件啊
赞一下(1人) 回复此楼
2楼2013-05-03 11:22:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zeppe

金虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
从选择条件 if (iter.gt.itermax-itervel+1) then 来看，只判断了迭代次数，可能还需要加上某个值小于容许误差。
赞一下 回复此楼
3楼2013-05-03 11:57:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xingfendou

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by 夕阳西下 at 2013-05-03 11:22:35
do 1130 i=1,iycur
没有判断收敛的条件啊

求收敛的条件在求电势（迭代求泊松方程）中，之后求电场，再求漂移速度，然后才求电流，电势和电场分布还行，漂移速度分布大致还行，但是很有毛刺，这是老师给的一程序，有点老了，真的很想做出来
赞一下 回复此楼
4楼2013-05-03 12:16:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xingfendou

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by zeppe at 2013-05-03 11:57:00
从选择条件 if (iter.gt.itermax-itervel+1) then 来看，只判断了迭代次数，可能还需要加上某个值小于容许误差。

电流的收敛应该是由电势，电场，漂移速度的结果引起的，整个程序太大了，但是我感觉应该是小问题，就是不知道到底错哪了
赞一下 回复此楼
5楼2013-05-03 12:23:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zeppe

金虫 (小有名气)

【答案】应助回帖

★ ★ ★
xingfendou: 金币+3, ★有帮助 2013-05-05 11:43:33
引用回帖:
5楼: Originally posted by xingfendou at 2013-05-03 12:23:15
电流的收敛应该是由电势，电场，漂移速度的结果引起的，整个程序太大了，但是我感觉应该是小问题，就是不知道到底错哪了...

可以试试用比较简单的能知道每一步该出什么结果的例子测试，在过程中加几个断点调试一下，看在哪个部分出错了，然后再对那个部分的程序进行调试检查。
赞一下 回复此楼
6楼2013-05-03 13:10:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xingfendou 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭