java编程.docVIP

作业： 已知二元函数 试以（0，0）为初始点，分别使用Newton法和最速下降法求函数f(x)的极小点，要求编写计算程序求解（编程语言不限），并给出计算结果。 注：（1）程序及计算结果需要用A4纸打印。 （2）作业在11月22日上交。 （3）自己独立完成，严禁互相抄袭。若发现抄袭，平时成绩按照零分计算。 用Newton法和最速下降法求函数f(x)的极小点 java实现编程 2010-11-14 11:28 同学让帮忙做数值分析，要求用编码实现Newton法和最速下降法，觉得对不少朋友应该有用，所以代码写完之后贴到这里，供大家参考。。。 题目如下： 已知二元函数 试以（0，0）为初始点，分别使用Newton法和最速下降法求函数f(x)的极小点，要求编写计算程序求解（编程语言不限），并给出计算结果。 算法描述如下： 用java实现的算法代码如下： 1.牛顿法： package Test_Math; public class NewTon_way { /*** @param args* @这是牛顿法求极值* @author 竹下听箫 * @date 2010.11.13*/ /*** Old_fun()是原函数*/public double Old_fun(int coeff[], double x[]) {????? return (double)coeff[0]*Math.pow(x[0],2)+ coeff[1]*Math.pow(x[1],2)+coeff[2]*x[0]*x[1]+ coeff[3]*x[0]+coeff[4]*x[1];} /***grads()是求得的导函数*/public void grads(int coeff[], double x[],double grads_x[]) {?????? ?? grads_x[0] = (double)2*coeff[0]*x[0]+coeff[2]*x[1]+coeff[3];?? grads_x[1] = (double)2*coeff[1]*x[1]+coeff[2]*x[0]+coeff[4];}/***求矩阵*/public void Matrix(int coeff[], double temp_matrix[]){?? temp_matrix[0] = (double)2.0*coeff[0];?? temp_matrix[1] = (double)coeff[2];?? temp_matrix[2] = (double)coeff[2];?? temp_matrix[3] = (double)2.0*coeff[1];}/***对应的逆矩阵*/public void invertH(double temp_matrix[], double invert[]){?? double deter;????? deter = temp_matrix[0]*temp_matrix[3]-temp_matrix[1]*temp_matrix[2];?? invert[0] = temp_matrix[3]/deter;?? invert[1] = -temp_matrix[1]/deter;????? invert[2] = -temp_matrix[2]/deter;?? invert[3] = temp_matrix[0]/deter;} public static void main(String[] args) {?? // TODO Auto-generated method stub?? double eps = 0.00001 ; //精度值?? int coeff[] = new int[5];//这是初始函数的系数值数组，可以自己设定?? coeff[0] =2;?? coeff[1] =1;?? coeff[2] =2;?? coeff[3] =1;?? coeff[4] =-1;???? ?? double x[] = new double[2];//初始的值???? x[0] = 0;???? x[1] = 0;???? double grads_x[] = new double[2];???? double matrix[] = new double[4];???? double invert[] = new double[4];????? double last_val = 0;?? NewTon_way newTon = new NewTon_way(); ?? newT