早教吧作业答案频道 -->数学-->
谁能用Matlab帮我解如下多元非线性方程组,解出结果来的另有重谢!这里面的8个方程,8个字母的全是未知数。
题目详情
谁能用Matlab帮我解如下多元非线性方程组,解出结果来的另有重谢!
这里面的8个方程,8个字母的全是未知数。
这里面的8个方程,8个字母的全是未知数。
▼优质解答
答案和解析
myfun.m
function F = myfun(x)
F=500*800-x(4)*sin(x(7))*160;
500*180-x(3)*sin(x(8))*160;
500-x(4)*cos(x(7)-17.7*pi/180)-x(3)*cos(x(8)+17.7*pi/180);
x(4)*sin(x(7)-17.7*pi/180)-x(3)*sin(x(8)+17.7*pi/180);
500*280-(x(5)- x(1))*150;
x(1)+ x(2)-x(5)-x(6);
0.15* x(1)-x(3);
0.15* x(2)-x(4);
------------------------------
findsolution.m
clear all;
warning off;
A=[];
for i=-pi:0.01:pi
options = optimoptions('fsolve','Display','iter');
x0=[15000/pi,15000/pi,2250/pi,2250/pi,15000/pi + 1125/2,4.2121e+03,i,i];
[x,q]=fsolve(@myfun,x0,options);
L0=[x,i,q];
A=[A;x,i,q];
end
Q=A(:,end);
H=find(Q==min(Q));
reslut=A(H,:)
找不到好的解.初值比较难决定.试了很多还是不尽人意.
以上得到的解为:
reslut =
1.0e+05 *
Columns 1 through 7
0.0477 0.0477 0.0072 -0.0171 0.0534 0.0421 -0.0138
Columns 8 through 10
-0.0000 -0.0000 1.2633
1到8列是解,最后一列是误差.误差太大.
继续努力.
刚改进了算法,得到一个误差比较小的
chuzhi =
-3.5730 3.5699 10.7127 17.8556 24.9984 32.1413 39.2842 46.4270
-----------solution------------
ans =
1.0e+03 *
-0.0036 0.0036 0.0107 2.3407 0.0250 0.0321 -0.5513 0.0464
又得到一个更好的解,这个比较靠谱了.哈哈
-----------------initial value------------
originaldata =
24.1804 31.0785 22.2690 0.0626 34.4995 27.2732 3.5825 3.4617
-----------solution------------
reslut =
1.0e+03 *
0.0242 0.0311 0.0223 -5.6378 0.0345 0.0273 -0.7441 0.0035
得到一个更好的.哈哈
riginaldata =
26.0918 29.0874 37.6150 0.0110 30.5449 12.9429 33.4451 2.3592
-----------solution------------
reslut =
1.0e+03 *
0.0261 0.0291 0.0376 2.6586 0.0305 0.0129 0.0201 0.0024
function F = myfun(x)
F=500*800-x(4)*sin(x(7))*160;
500*180-x(3)*sin(x(8))*160;
500-x(4)*cos(x(7)-17.7*pi/180)-x(3)*cos(x(8)+17.7*pi/180);
x(4)*sin(x(7)-17.7*pi/180)-x(3)*sin(x(8)+17.7*pi/180);
500*280-(x(5)- x(1))*150;
x(1)+ x(2)-x(5)-x(6);
0.15* x(1)-x(3);
0.15* x(2)-x(4);
------------------------------
findsolution.m
clear all;
warning off;
A=[];
for i=-pi:0.01:pi
options = optimoptions('fsolve','Display','iter');
x0=[15000/pi,15000/pi,2250/pi,2250/pi,15000/pi + 1125/2,4.2121e+03,i,i];
[x,q]=fsolve(@myfun,x0,options);
L0=[x,i,q];
A=[A;x,i,q];
end
Q=A(:,end);
H=find(Q==min(Q));
reslut=A(H,:)
找不到好的解.初值比较难决定.试了很多还是不尽人意.
以上得到的解为:
reslut =
1.0e+05 *
Columns 1 through 7
0.0477 0.0477 0.0072 -0.0171 0.0534 0.0421 -0.0138
Columns 8 through 10
-0.0000 -0.0000 1.2633
1到8列是解,最后一列是误差.误差太大.
继续努力.
刚改进了算法,得到一个误差比较小的
chuzhi =
-3.5730 3.5699 10.7127 17.8556 24.9984 32.1413 39.2842 46.4270
-----------solution------------
ans =
1.0e+03 *
-0.0036 0.0036 0.0107 2.3407 0.0250 0.0321 -0.5513 0.0464
又得到一个更好的解,这个比较靠谱了.哈哈
-----------------initial value------------
originaldata =
24.1804 31.0785 22.2690 0.0626 34.4995 27.2732 3.5825 3.4617
-----------solution------------
reslut =
1.0e+03 *
0.0242 0.0311 0.0223 -5.6378 0.0345 0.0273 -0.7441 0.0035
得到一个更好的.哈哈
riginaldata =
26.0918 29.0874 37.6150 0.0110 30.5449 12.9429 33.4451 2.3592
-----------solution------------
reslut =
1.0e+03 *
0.0261 0.0291 0.0376 2.6586 0.0305 0.0129 0.0201 0.0024
看了 谁能用Matlab帮我解如下...的网友还看了以下:
已知集合M具有性质:若A属于M,则2A属于M,现已知-1属于M,则下列元素一定在M中的A-0.5B0 2020-03-30 …
某人买了50元的乘车月票卡,如果此人乘车的次数用m表示,则记录他每次乘车后的余额n元,如下表:乘车 2020-05-13 …
某人买了50元的乘车月票卡,如果此人乘车的次数用m表示,则记录他每次乘车后的余额n元,如下表:乘车 2020-05-13 …
在国家政策的宏观调控下,太仓市某高档小区商品房成交价由今年3月底的8600元/m^2下降到5月底的 2020-05-20 …
已经价格的单位是元/m2,如何求出价格等于元/m?已知,条件1:某圆管表面制作价格0.17元/cm 2020-07-06 …
已知a为IIA族元素,b为ⅢA族元素,它们的原子序数分别为m和n,且A.b为同一周期元素,下列关系 2020-07-29 …
如图所示,将细绳绕过两定滑轮A和B,绳的两端各挂一个质量为m的砝码,AB的中点C处挂一质量为M的物体 2020-11-01 …
设数集M同时满足条件:(1)M中不含元素-1,0,1,(2)若a∈M,则∈M.有下列结论:①集合M中 2020-11-08 …
(2014•中山模拟)四种短周期元素在周期表中的位置如图,其中只有M为金属元素.下列说法不正确的是( 2020-11-12 …
(2012•株洲)M元素的一种微粒的结构示意图为,已知:(1)M元素的单质是空气的重要组成部分,约占 2020-11-12 …