| 查看: 1165 | 回复: 5 | ||
破碎东瓜金虫 (初入文坛)
|
[求助]
谁能帮我把一段matlab程序转为fortran啊！！！（很短的一段）
|
|
while err>errlim s=find(P>0); sn=find(P<=0); dd=real(ifft2(fA.*fft2(P,n*2,n*2))); u=dd(1:n,1:n)*2.0*dx/pi^2; rk=h+u; do=mean(rk(s)); rk=rk-do; G=sum(sum(rk(s).*rk(s))); pk(s)=rk(s)+G/Gold*pk(s); pk(sn)=0; Gold=G; dd=real(ifft2(fA.*fft2(pk,2*n,2*n))); qk=dd(1:n,1:n)*2.0*dx/pi^2; rb=mean(qk(s)); qk=qk-rb; dp=sum(sum(rk(s).*pk(s)))/sum(sum(qk(s).*pk(s))); P(s)=P(s)-dp*pk(s); s=find(P<=0);P(s)=0; sol=find((P==0)&(rk<0)); P(sol)=P(sol)-dp*rk(sol); W=sum(sum(P))*dx*dx; P=charge/W*P; err=sqrt(Gold*dx*dx); end |
回复此楼» 猜你喜欢
云南大学材料与能源学院解琳课题组钙钛矿博士招生
已经有4人回复
-
博士申请
已经有6人回复
-
评审有感
已经有11人回复
-
面上本子正文33页，违规吗？会被低分嘛？
已经有10人回复
-
国自然上会要求
已经有9人回复
-
今年审到国自然15份，谈谈感受
已经有16人回复
-
上海大学实验技术岗位非升即走
已经有8人回复
-
考博自荐
已经有6人回复
-
青C资助名额大幅增加！
已经有16人回复
-
重磅！青年科学基金项目（C类）资助增幅预计超过50%
已经有10人回复
wxws.2008
木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 5015.5
- 散金: 256
- 红花: 9
- 帖子: 354
- 在线: 147.9小时
- 虫号: 639701
- 注册: 2008-10-29
- 性别: GG
- 专业: 荒漠化与水土保持
2楼2015-12-12 09:57:31
破碎东瓜
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1520.2
- 散金: 53
- 帖子: 46
- 在线: 32.1小时
- 虫号: 2231948
- 注册: 2013-01-09
- 性别: GG
- 专业: 机械摩擦学与表面技术
3楼2015-12-12 15:36:58
wxws.2008
木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 5015.5
- 散金: 256
- 红花: 9
- 帖子: 354
- 在线: 147.9小时
- 虫号: 639701
- 注册: 2008-10-29
- 性别: GG
- 专业: 荒漠化与水土保持
4楼2015-12-12 18:07:54
lishouyin
铁杆木虫 (著名写手)
- 应助: 25 (小学生)
- 金币: 12773.3
- 散金: 614
- 红花: 10
- 沙发: 1
- 帖子: 2824
- 在线: 793.6小时
- 虫号: 1983857
- 注册: 2012-09-07
- 性别: GG
- 专业: 电路与系统
5楼2015-12-13 10:32:35
破碎东瓜
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1520.2
- 散金: 53
- 帖子: 46
- 在线: 32.1小时
- 虫号: 2231948
- 注册: 2013-01-09
- 性别: GG
- 专业: 机械摩擦学与表面技术
6楼2015-12-13 14:49:55
20