找回密码
 注册
Simdroid-非首页
查看: 39|回复: 0

[3. Fortran] (0<u<=1) (advice u>0.5)

[复制链接]
发表于 2007-4-1 12:11:17 | 显示全部楼层 |阅读模式 来自 黑龙江哈尔滨
以下是一个matlab代码,但其中几句(如标题上的)是fortran语言,请教大家是什么意思
function sa=art(a,t,s1,u,ds0)
% t=a*s m条射线 n个网格
% t is travel time           M*1
% a is travel route          M*N
% s is slowness              N*1
% s1 is first model          N*1
% u is loosen factor (0<u<=1)  (advice u>0.5)
% ds0 is precision of remain
[M,N]=size(a);
s=ones(M,N);s(1,=s1';
k=(sum((a.*a)'))';
ds=2*ds0;
    while abs(max(ds))>ds0
        for i=1:M-1
            t1(i,1)=sum(a(i,:).*s(i,:));
            dt(i,1)=t(i,1)-t1(i,1);
            ds=dt(i,1)*a(i,:)/k(i,1);
            s(i+1,:)=s(i,:)+u*ds;
        end
    end
sa=(s(M,:))';
您需要登录后才可以回帖 登录 | 注册

本版积分规则

Simapps系列直播

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-11-1 17:37 , Processed in 0.026219 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表