| 查看: 1101 | 回复: 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 |
回复此楼» 猜你喜欢
805600专硕材料与化工348分求调剂
已经有6人回复
-
重庆大学材料与化工085600，初试370+，求求调剂建议
已经有5人回复
-
309求调剂
已经有3人回复
-
一志愿北交大材料工程，总分358
已经有4人回复
-
302求调剂一志愿北航070300，本科郑大化学
已经有11人回复
-
化学工程专硕324分，一志愿中国矿业大学求调剂
已经有7人回复
-
材料求调剂
已经有11人回复
-
【求调剂】新能源材料本科，一志愿211，初试321
已经有4人回复
-
336求调剂
已经有5人回复
-
材料求调剂
已经有8人回复
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
