| 查看: 1103 | 回复: 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 |
回复此楼» 猜你喜欢
085600，材料与化工321分，求调剂
已经有13人回复
-
311求调剂
已经有8人回复
-
348求调剂
已经有3人回复
-
0805求调剂
已经有8人回复
-
266分，一志愿电气工程，本科材料，求材料专业调剂
已经有3人回复
-
385分 生物学（071000）求调剂
已经有6人回复
-
材料调剂
已经有3人回复
-
一志愿北交大材料工程总分358
已经有9人回复
-
311求调剂
已经有12人回复
-
319求调剂
已经有11人回复
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
