24小时热门版块排行榜

返回列表

lolo123

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 2670.9
散金: 88
帖子: 221
在线: 112.5小时
虫号: 1454002
注册: 2011-10-21
性别: GG
专业: 机械结构强度学

[求助] matlab程序的一个问题求助

一个计算程序前面的一段
syms cell A C D lmt a0 h1 k1 l1 h2 k2 l2 h3 k3 l3 h4 k4 l4 h5 k5 l5 h6 k6 l6 j101 j102 j103 th1 th2 th3 bgk101 bgk102 bgk103 bgk1 bgk2 bgk3 aef1 aef2 aef3 c0 gama1 gama2 gama3 delt1 delt2 delt3 d1 d2 d3 IN;
lmt= 0.15406;
t=96;
IN=[ 101 43.246 0.082;
      102 54.356 0.087;
      103 70.119 0.119;
      002 36.325 0.091;
      100 39.016 0.08;
      110 70.659 0.107];

h1=fix(IN(1,1)/100);
k1=fix(IN(1,1)/10)-fix(IN(1,1)/100)*10;
l1=fix(IN(1,1))-fix(IN(1,1)/10)*10;
j101=IN(1,2);
bgk101=IN(1,3);
h2=fix(IN(2,1)/100);
k2=fix(IN(2,1)/10)-fix(IN(2,1)/100)*10;
l2=fix(IN(2,1))-fix(IN(2,1)/10)*10;
j102=IN(2,2);
bgk102=IN(2,3);
h3=fix(IN(3,1)/100);
k3=fix(IN(3,1)/10)-fix(IN(3,1)/100)*10;
l3=fix(IN(3,1))-fix(IN(3,1)/10)*10;
j102=IN(3,2);
bgk102=IN(3,3);

th1=IN(4,2);bgk1=IN(4,3);
th2=IN(5,2);bgk2=IN(5,3);
th3=IN(6,2);bgk3=IN(6,3);

d1=lmt/(2*sin(j101*pi/360));
d2=lmt/(2*sin(j102*pi/360));
d3=lmt/(2*sin(j103*pi/360));

aef1=h1^2+h1*k1+l1^2;
aef2=h2^2+h2*k2+l2^2;
aef3=h3^2+h3*k3+l3^2;

gama1=l1^2;
gama2=l2^2;
gama3=l3^2;

delt1=sin(j101*pi/180)^2*(1/sin(j101*pi/360)+1/(j101*pi/360));
delt2=sin(j102*pi/180)^2*(1/sin(j102*pi/360)+1/(j102*pi/360));
delt3=sin(j103*pi/180)^2*(1/sin(j103*pi/360)+1/(j103*pi/360));

j101=sin(j101*pi/360)^2;
j102=sin(j102*pi/360)^2;
j103=sin(j103*pi/360)^2;

M1=aef1*j101+aef2*j102+aef3*j103;
M2=gama1*j101+gama2*j102+gama3*j103;
M3=delt1*j101+delt2*j102+delt3*j103;

M11=aef1^2+aef2^2+aef3^2;
M12=gama1*aef1+gama2*aef2+gama3*aef3;
M13=delt1*aef1+delt2*aef2+delt3*aef3;

M21=aef1*gama1+aef2*gama2+aef3*gama3;
M22=gama1^2+gama2^2+gama3^2;
M23=gama1*delt1+gama2*delt2+gama3*delt3;

M31=aef1*delt1+aef2*delt2+aef3*delt3;
M32=gama1*delt1+gama2*delt2+gama3*delt3;
M33=delt1^2+delt2^2+delt3^2;

MP=[M11 M12 M13;
M21 M22 M23;
M31 M32 M33];
MA=[M1  M12 M13;
M2  M22 M23;
M3  M32 M33];
MC=[M11 M1  M13;
   M21 M2  M23;
   M31 M3  M33];
MD=[M11 M12 M1 ;
M21 M22 M2 ;
M31 M32 M3 ];
MP=det(MP);
MA=det(MA);
MC=det(MC);
MD=det(MD);

A=MA/MP;
C=MC/MP;
D=MD/MP;

a0=lmt/sqrt(3*A);
c0=lmt/sqrt(4*C);

if C<0
rh121=l1^2;
rh122=l2^2;
Q121=3*a0^2/(4*d1^2);
Q122=3*a0^2/(4*d2^2);
P1211=2;
P1212=rh121+rh122;
P1221=P1212;
P1222=rh121^2+rh122^2;
P121=Q121+Q122;
P122=Q121*rh121+Q122*rh122;
PP12=[P1211 P1212;
P1221 P1222];
PQ12=[P121 P1212;
P122 P1222];
PR12=[P1211 P121;
P1221 P122];
PP12=det(PP12);
PQ12=det(PQ12);
PR12=det(PR12);
R12=PR12/PP12;
c012=sqrt(0.75*a0^2/R12);

rh131=l1^2;
rh132=l3^2;
Q131=3*a0^2/(4*d1^2);Q132=3*a0^2/(4*d3^2);
P1311=2;
P1312=rh131+rh132;
P1321=P1312;
P1322=rh131^2+rh132^2;
P131=Q131+Q132;
P132=Q131*rh131+Q132*rh132;
PP13=[P1311 P1312;
P1321 P1322];
PQ13=[P131 P1312;
P132 P1322];
PR13=[P1311 P131;
P1321 P132];
PP13=det(PP13);
PQ13=det(PQ13);
PR13=det(PR13);
R13=PR13/PP13;
c013=sqrt(0.75*a0^2/R13);

rh231=l2^2;rh232=l3^2;
Q231=3*a0^2/(4*d2^2);
Q232=3*a0^2/(4*d3^2);
P2311=2;
P2312=rh231+rh232;
P2321=P2312;P2322=rh231^2+rh232^2;
P231=Q231+Q232;
P232=Q231*rh231+Q232*rh232;
PP23=[P2311 P2312;
P2321 P2322];
PQ23=[P231 P2312;
P232 P2322];
PR23=[P2311 P231;
P2321 P232];
PP23=det(PP23);
PQ23=det(PQ23);
PR23=det(PR23);
R23=PR23/PP23;
c023=sqrt(0.75*a0^2/R23);
c0=(c012+c013+c023)/3;
end

a0=a0*10;
c0=c0*10;

format long
cell=[a0 c0];
这时出现 Undefined function or method 'lt' for input arguments of type 'sym'. 我的程序里没有lt这个参数，为什么会出现这种情况？我的金币不是很多，先谢谢了。

回复此楼