Matlab编程问题:编写牛顿法程序来求解下面的方程组:(x-2)^2+(y-3+2x)^2=5 2(x-3)^2+(y
来源:学生作业帮 编辑:神马作文网作业帮 分类:综合作业 时间:2024/11/16 04:00:37
Matlab编程问题:编写牛顿法程序来求解下面的方程组:(x-2)^2+(y-3+2x)^2=5 2(x-3)^2+(y/3)^2=4
% 用牛顿迭代法解非线性方程组
% F1=(x-2)^2+(y-3+2x)^2-5=0
% F2=2(x-3)^2+(y/3)^2-4=0
% F=[F1;F2]=[5*x^2 + 4*x*y - 16*x + y^2 - 6*y + 8
% 2*x^2 - 12*x + y^2/9 +14 ]
%迭代初值设为:X0=[x,y]=[0,0]'
% X(k+1)=X(k)+delt_X
% Jaco_F(X(k))*delt_X(k)=-F(X(k))
% Jaco_F为Jacobian矩阵
% Jaco_F=[10*x+4*y-16 2*y+4*x-6
% 4*x-12 2*y/9]
%error 为 ||X(k+1)-X(k)||,当error小于0.0001时,迭代停止
x=0;
y=0;
error=10;
while(error>0.0001)
F = [ 5*x^2+4*x*y-16*x+y^2-6*y+8
2*x^2-12*x+y^2/9+14 ];
Jaco_F = [10*x+4*y-16 2*y+4*x-6
4*x-12 2*y/9];
X = [x y]' - (inv(Jaco_F))*F; %迭代后的X值
error=norm(X-[x,y]'); % error=||X(k+1)-X(k)||
x=X(1);
y=X(2);
end
x
y
%%%%%%%%%%%%%到此结束%%%%%%%%%%%%%
程序结果:
x =
1.7362
y =
-2.6929
% F1=(x-2)^2+(y-3+2x)^2-5=0
% F2=2(x-3)^2+(y/3)^2-4=0
% F=[F1;F2]=[5*x^2 + 4*x*y - 16*x + y^2 - 6*y + 8
% 2*x^2 - 12*x + y^2/9 +14 ]
%迭代初值设为:X0=[x,y]=[0,0]'
% X(k+1)=X(k)+delt_X
% Jaco_F(X(k))*delt_X(k)=-F(X(k))
% Jaco_F为Jacobian矩阵
% Jaco_F=[10*x+4*y-16 2*y+4*x-6
% 4*x-12 2*y/9]
%error 为 ||X(k+1)-X(k)||,当error小于0.0001时,迭代停止
x=0;
y=0;
error=10;
while(error>0.0001)
F = [ 5*x^2+4*x*y-16*x+y^2-6*y+8
2*x^2-12*x+y^2/9+14 ];
Jaco_F = [10*x+4*y-16 2*y+4*x-6
4*x-12 2*y/9];
X = [x y]' - (inv(Jaco_F))*F; %迭代后的X值
error=norm(X-[x,y]'); % error=||X(k+1)-X(k)||
x=X(1);
y=X(2);
end
x
y
%%%%%%%%%%%%%到此结束%%%%%%%%%%%%%
程序结果:
x =
1.7362
y =
-2.6929
matlab编程求解微分方程组:x=x*y+t,y=x^2*sin(y)+t^2
c语言编程 编写程序求下列表达式的值 x+a%3*(x+y)%2/4.设x=2.5,a=7,y=4
请matlab高手教教 matlab编写一段程序,画出函数y=x^3*e^(-x^2),-1
基于MATLAB的遗传算法程序 对函数y=(x)1^2-(x)2^2+3 求解全局最大值.
怎么编写绘制 z=x^2+y^2图形的Matlab程序.
求方程组x^2+y^2=1','x*y=2的解的matlab程序(使用solve)
C语言编程题:输入X和Y的值,编写程序计算f(x,n)=x^2+x^3-x^4+...(-1)^n-1x^n,其中用函数
5x+y=2 x-3y=4 二元一次方程组求解
MATLAB求方程y’= y - 2x/y,y(0) = 1 (0≤x≤1,h = 0.1) 的数值解,要求编写程序
用c++编写高斯消元法求解线性方程组的计算程序,并计算下例:2X+2Y+3Z=3、4X+7Y+7Z=1、-2X+4Y+5
方程组{上面的4x+3y=17下面2x+5y=19
编写一个程序,输入x的值,输出函数y=5x^5+2x^4+3x^3+2x^2+x