matlab代码转换成fortran代码，是关于求解方程的！

U=16;
p=1025;s=20;c=8;lmda=2.5/c;g=9.8;
mf=pi*p*s^2*c^2/4/(s^2+c^2)^(1/2)/5;
A33f=2*mf;
xf=83.3-c/2;
A35f=-2*mf*xf;
A55f=2*mf*xf^2;
cl=5.7/(1+2/lmda);
C35f=p*U^2*s*cl;
C55f=-p*U^2*s*cl*xf;
A53f=A35f;
B33f=p*U*s*cl;
B35f=-p*U*s*cl*xf;
B55f=p*U*s*cl*xf^2;
B53f=B35f; M=0.80584*10^4*1025;C33=2.13476*10^7;C53=1.6096*10^8;C55=3.312161*10^10;C35=C53; Isita=1.4363*10^10;
ZZ=zeros(53,7); QQ=zeros(53,7); zz=zeros(53,7); qq=zeros(53,7)
for m=1:1:4;
if m<=1
beta =pi*2/3;
else
beta=((m-1)*15+135)*2*pi/360;
end for
n=1:1:51;
w=0.08+n*0.02;
A33=[];
A35=[];
A55=[];
B33=[];
B35=[];
B53=[];
B55=[];
F3c=[];
F3s=[];
F5c=[];
F5s=[];
we=w-w^2/g*U*cos(beta);k=w^2/g;z=-5.58;
A11=(M+A33(n,m))*(-we^2)+C33-i*we*B33(n,m);
A12=(-we^2*(A35(n,m))+C35-i*we*B35(n,m));
B1=F3c(n,m)+i*F3s(n,m);
A21=C53-we^2*A53( n,m)-i*we*(B53( n,m));
A22=(-we^2*(Isita+A55( n,m))+C55)-i*we*(B55(n,m));
B2=F5c(n,m)+F5s(n,m)*i;
AA=[A11, A12;A21, A22;];
b=[B1;B2;];
xx=AA\b;
ZZ(n,m)=abs(xx(1));
QQ(n,m)=abs(xx(2));

a11=(M+A33(n,m)+A33f)*(-we^2)+C33-i*we*(B33(n,m)+B33f);
a12=-we^2*(A35(n,m)+A35f)+C35+C35f-i*we*(B35(n,m)+B35f);
b1=F3c(n,m)+i*F3s(n,m)+p*U*s*cl*(-w)*exp(k*z+i*k*xf*cos(beta));
a21=C53-we^2*(A53( n,m)+A53f)-i*we*(B53( n,m)+B53f);
a22=-we^2*(Isita+A55( n,m)+A55f)+C55+C55f-i*we*(B55(n,m)+B55f);
b2=F5c(n,m)+F5s(n,m)*i-p*U*s*cl*xf*(-w)*exp(k*z+i*k*xf*cos(beta));
aa=[a11, a12;a21, a22;];
d=[b1;b2;];
yy=aa\d;
zz(n,m)=abs(yy(1));
qq(n,m)=abs(yy(2));
end;
end;
QQ=QQ*360/2/pi;
qq=qq*360/2/pi;
哈尔滨工程大学本科生毕业论文

YY频道7450欢迎您
ansys、flux之类的功能比较正规点， 设置一下就可以做分析， 但是你有什么新奇点的想法大概是很难实现的了
matlab有个FEM函数包， 也许做这方面能容易些吧。

如果你是公司里要做事情， 用ansys flux这种现成的就好了。 如果是搞研究要琢磨算法或者写论文的话也许matlab好一些。 如果研究完了要实用的话， 可以等matlab做通了转fortran， 速度能快一些

写代码比较省心， 要是能设置好分布式计算的话， 算几天什么数据都有了 。 而且用ansys的话论文里对算法、精度、开销、细节之类的也不太好分析吧
我是说ansys， flux之类的东西方法都是固定的， 不是用来研究的， 如果就具体某类问题你选了一个比较适合的方法，做了一些适当的改进的话， 自己写个程序就要比用这些更实用点。
是的

你好！ 你想问的是什么？我的论文也是MATLAB 如有疑问，请追问。