用matlab的简单计算clcclearsyms E sr2 ze1=0.7sr1=0.8Es=8e6J=0.02v=0.35a0=4.5e6a=1.5r0=1n=e1/(1+e1)z0=0.1z=10sr2=solve((1-2*v)*(1+v)*(2*a*a0*r0^2/z^3*[1+r0^2/z^2]^(-a-1)*z0)/n/(2*Es*(1+(e1*sr1/sr2))^2/(1+e1)/(2+(e1*sr1/sr2)+e1))/(1-v)/(1-sr1+sr1*
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/15 09:30:42
用matlab的简单计算clcclearsyms E sr2 ze1=0.7sr1=0.8Es=8e6J=0.02v=0.35a0=4.5e6a=1.5r0=1n=e1/(1+e1)z0=0.1z=10sr2=solve((1-2*v)*(1+v)*(2*a*a0*r0^2/z^3*[1+r0^2/z^2]^(-a-1)*z0)/n/(2*Es*(1+(e1*sr1/sr2))^2/(1+e1)/(2+(e1*sr1/sr2)+e1))/(1-v)/(1-sr1+sr1*
用matlab的简单计算
clc
clear
syms E sr2 z
e1=0.7
sr1=0.8
Es=8e6
J=0.02
v=0.35
a0=4.5e6
a=1.5
r0=1
n=e1/(1+e1)
z0=0.1
z=10
sr2=solve((1-2*v)*(1+v)*(2*a*a0*r0^2/z^3*[1+r0^2/z^2]^(-a-1)*z0)/n/(2*Es*(1+(e1*sr1/sr2))^2/(1+e1)/(2+(e1*sr1/sr2)+e1))/(1-v)/(1-sr1+sr1*J)-(sr2*(1-sr1+sr1*J)/sr1/(1-sr2+sr2*J)-1))
vpa(sr2,5)
我想使z=1:1:10
搞成一个循环,
用matlab的简单计算clcclearsyms E sr2 ze1=0.7sr1=0.8Es=8e6J=0.02v=0.35a0=4.5e6a=1.5r0=1n=e1/(1+e1)z0=0.1z=10sr2=solve((1-2*v)*(1+v)*(2*a*a0*r0^2/z^3*[1+r0^2/z^2]^(-a-1)*z0)/n/(2*Es*(1+(e1*sr1/sr2))^2/(1+e1)/(2+(e1*sr1/sr2)+e1))/(1-v)/(1-sr1+sr1*
clc
clear
syms E sr2 z
e1=0.7
sr1=0.8
Es=8e6
J=0.02
v=0.35
a0=4.5e6
a=1.5
r0=1
n=e1/(1+e1)
z0=0.1
for z=1:10
sr2=solve((1-2*v)*(1+v)*(2*a*a0*r0^2/z^3*[1+r0^2/z^2]^(-a-1)*z0)/n/(2*Es*(1+(e1*sr1/sr2))^2/(1+e1)/(2+(e1*sr1/sr2)+e1))/(1-v)/(1-sr1+sr1*J)-(sr2*(1-sr1+sr1*J)/sr1/(1-sr2+sr2*J)-1))
vpa(sr2,5)
end