MATLAB求解符号线性方程组时,结果误差很大?
MATLAB求解符号线性方程组时,结果误差很大?
syms F1 F2 F3 F4 F5 C1 C2 C3 C4 E I D2 D3 D4 W0 L1 L2 L3 L4;
f1=C1+(F4*L1^2)/(2*E*I)-C2;
f6=(F4-F5)*(L1+L2)^2/(2*E*I)-F4*L1*(L1+L2)/(E*I)+C2-((F4-F3-F5)*(L1+L2)^2)/(2*E*I)-((F3-F4)*L1+F3*L2)*(L1+L2)/(E*I)-C3;
f7=(F4-F3-F5)*(L1+L2+L3)^2/(2*E*I)+((F3-F4)*L1+F3*L2)*(L1+L2+L3)/(E*I)+C3-F1*(L1+L2+L3)^2/(2*E*I)+(L1+L2+L3+L4)*F1*(L1+L2+L3)/(E*I)-C4;
f2=-(F5*L1^3)/(6*E*I)+C1*L1;
f3=(F4-F5)*L1^3/(6*E*I)-(F4*L1^3)/(2*E*I)+C2*L1+D2;
f4=(F4-F5)*(L1+L2)^3/(6*E*I)-(F4*L1*(L1+L2)^2)/(2*E*I)+C2*(L1+L2)+D2+W0;
f5=(F4-F3-F5)*(L1+L2)^3/(6*E*I)+((F3-F4)*L1+F3*L2)*(L1+L2)^2/(2*E*I)+C3*(L1+L2)+D3+W0; %x=L1+L2;W=W0
f8=(F4-F3-F5)*(L1+L2+L3)^3/(6*E*I)+((F3-F4)*L1+F3*L2)*(L1+L2+L3)^2/(2*E*I)+C3*(L1+L2+L3)+D3; %x=L1+L2+L3;W=0
f9=F1*(L1+L2+L3)^3/(6*E*I)-(L1+L2+L3+L4)*F1*(L1+L2+L3)^2/(2*E*I)+C4*(L1+L2+L3)+D4; %x=L1+L2+L3;W=0
f10=-F1*(L1+L2+L3+L4)^3/(3*E*I)+C4*(L1+L2+L3+L4)+D4; %x=L1+L2+L3+L4;W=0
f11=-F1+F2-F3+F4-F5;
f12=-F3*(L3+L4)+F4*(L2+L3+L4)-F5*(L1+L2+L3+L4)+F2*L4;
[F1,F2,F3,F4,F5,C1,C2,C3,C4,D2,D3,D4]=solve(f1,f2,f3,f4,f5,f6,f7,f8,f9,f10,f11,f12,'F1','F2','F3','F4','F5','C1','C2','C3','C4','D2','D3','D4')
上面是我的代码,求出的结果带入f1~f12,发现都不为0,这是为什么?
syms F1 F2 F3 F4 F5 C1 C2 C3 C4 E I D2 D3 D4 W0 L1 L2 L3 L4;
f1=C1+(F4*L1^2)/(2*E*I)-C2;
f6=(F4-F5)*(L1+L2)^2/(2*E*I)-F4*L1*(L1+L2)/(E*I)+C2-((F4-F3-F5)*(L1+L2)^2)/(2*E*I)-((F3-F4)*L1+F3*L2)*(L1+L2)/(E*I)-C3;
f7=(F4-F3-F5)*(L1+L2+L3)^2/(2*E*I)+((F3-F4)*L1+F3*L2)*(L1+L2+L3)/(E*I)+C3-F1*(L1+L2+L3)^2/(2*E*I)+(L1+L2+L3+L4)*F1*(L1+L2+L3)/(E*I)-C4;
f2=-(F5*L1^3)/(6*E*I)+C1*L1;
f3=(F4-F5)*L1^3/(6*E*I)-(F4*L1^3)/(2*E*I)+C2*L1+D2;
f4=(F4-F5)*(L1+L2)^3/(6*E*I)-(F4*L1*(L1+L2)^2)/(2*E*I)+C2*(L1+L2)+D2+W0;
f5=(F4-F3-F5)*(L1+L2)^3/(6*E*I)+((F3-F4)*L1+F3*L2)*(L1+L2)^2/(2*E*I)+C3*(L1+L2)+D3+W0; %x=L1+L2;W=W0
f8=(F4-F3-F5)*(L1+L2+L3)^3/(6*E*I)+((F3-F4)*L1+F3*L2)*(L1+L2+L3)^2/(2*E*I)+C3*(L1+L2+L3)+D3; %x=L1+L2+L3;W=0
f9=F1*(L1+L2+L3)^3/(6*E*I)-(L1+L2+L3+L4)*F1*(L1+L2+L3)^2/(2*E*I)+C4*(L1+L2+L3)+D4; %x=L1+L2+L3;W=0
f10=-F1*(L1+L2+L3+L4)^3/(3*E*I)+C4*(L1+L2+L3+L4)+D4; %x=L1+L2+L3+L4;W=0
f11=-F1+F2-F3+F4-F5;
f12=-F3*(L3+L4)+F4*(L2+L3+L4)-F5*(L1+L2+L3+L4)+F2*L4;
[F1,F2,F3,F4,F5,C1,C2,C3,C4,D2,D3,D4]=solve(f1,f2,f3,f4,f5,f6,f7,f8,f9,f10,f11,f12,'F1','F2','F3','F4','F5','C1','C2','C3','C4','D2','D3','D4')
上面是我的代码,求出的结果带入f1~f12,发现都不为0,这是为什么?
赞 creak 幼苗
共回答了13个问题采纳率：84.6% 举报
原因使用solve函数求解方程组,如果返回多个参数,则参数是按照字母顺序排列,而与你在输入参数表中指定的顺序无关. 可参考关于solve函数的文档： Three different types of output are possible. For one equation and one
output, the resulting solution is returned, with multiple solutions to
a nonlinear equation in a symbolic vector. For several equations and
an equal number of outputs, the results are sorted in lexicographic
order and assigned to the outputs. For several equations and a single
output, a structure containing the solutions is returned.
建议一种方式是,修改solve函数的输出参数顺序：[C1,C2,C3,C4,D2,D3,D4,F1,F2,F3,F4,F5]=solve(f1,f2,f3,f4,f5,f6,f7,f8,f9,f10,f11,f12,'F1','F2','F3','F4','F5','C1','C2','C3','C4','D2','D3','D4')
或者,按照下面的方式返回参数：s=solve(f1,f2,f3,f4,f5,f6,f7,f8,f9,f10,f11,f12,'F1','F2','F3','F4','F5','C1','C2','C3','C4','D2','D3','D4');
fns = fieldnames(s);
for i=1:length(fns)
eval([fns{i} '=s.' fns{i}]);
end
output, the resulting solution is returned, with multiple solutions to
a nonlinear equation in a symbolic vector. For several equations and
an equal number of outputs, the results are sorted in lexicographic
order and assigned to the outputs. For several equations and a single
output, a structure containing the solutions is returned.
建议一种方式是,修改solve函数的输出参数顺序：[C1,C2,C3,C4,D2,D3,D4,F1,F2,F3,F4,F5]=solve(f1,f2,f3,f4,f5,f6,f7,f8,f9,f10,f11,f12,'F1','F2','F3','F4','F5','C1','C2','C3','C4','D2','D3','D4')
或者,按照下面的方式返回参数：s=solve(f1,f2,f3,f4,f5,f6,f7,f8,f9,f10,f11,f12,'F1','F2','F3','F4','F5','C1','C2','C3','C4','D2','D3','D4');
fns = fieldnames(s);
for i=1:length(fns)
eval([fns{i} '=s.' fns{i}]);
end
1年前
2
可能相似的问题
-
1年前1个回答
-
1年前1个回答
-
1年前2个回答
-
1年前2个回答
-
1年前5个回答
-
闭合导线坐标计算.方位角OK了,为什么纵坐标标增值误差很大?
1年前1个回答
你能帮帮他们吗