kindly help in the conversion of fortran code to matlab

1 次查看（过去 30 天）

显示更早的评论

Oluwaseyi Aliu 2020-4-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/516813-kindly-help-in-the-conversion-of-fortran-code-to-matlab

评论： Oluwaseyi Aliu 2020-4-10

采纳的回答： Walter Roberson

5 个评论
显示 3更早的评论隐藏 3更早的评论

Walter Roberson 2020-4-10

在 MATLAB Online 中打开

snul = 0.0;
snur = 0.0;
for j = 0:m
    rnul = (th(0,j) - th(1,j)) * n;
    rnur = (th(n-1,j) - th(n,j)) * n;
    snul = snul + rnul;
    snur = snur + rnur;
    fprintf('%g %g %g\n' j/m, rnul, rnur);
end
avnl = snul / m;
avnr = snur / m;
fprintf('%g %g %g\n', ra, avnl, avnr);

I do not see any 100 x 100 matrices, not unless you are saying that th is a matrix. In fortran it is possible but not all that common to declare arrays to start from index 0, or any other index. I would want to see the statements that defined th, such as

REAL*8 th(100,100)

but it could not be exactly like that if th is a matrix being indexed at 0.

Oluwaseyi Aliu 2020-4-10

sir, here is the initialisation;

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Walter Roberson 2020-4-10

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/516813-kindly-help-in-the-conversion-of-fortran-code-to-matlab#answer_425173

在 MATLAB Online 中打开

n = 100; m = 100;
c8 = length(0:8);
f = zeros(c8, n+1, m+1);
feq = zeros(c8, n+1, m+1);
rho = zeros(c8, n+1, m+1);
w = zeros(1, c8);
cx = zeros(1, c8);
cy = zeros(1, c8);
u = zeros(n+1, m+1);
v = zeros(n+1, m+1);
g = zeros(c8, n+1, m+1);
geq = zeros(c8, n+1, m+1);
th = zeros(n+1, m+1);
i = 0;
%then some code to read the data. 
%then
snul = 0.0;
snur = 0.0;
for j = 0:m
    rnul = (th(1+0,1+j) - th(1+1,1+j)) * n;
    rnur = (th(1+n-1,1+j) - th(1+n,1+j)) * n;
    snul = snul + rnul;
    snur = snur + rnur;
    fprintf('%g %g %g\n' j/m, rnul, rnur);
end
avnl = snul / m;
avnr = snur / m;
fprintf('%g %g %g\n', ra, avnl, avnr);