下载本文档
版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领
文档简介
1、 ，设一维的凸函数在区间a,b上存在一极值点处取得。黄金分割法在每次迭代中，计算在区间a,b的黄金分割点(1)和(2)的值 y1 和 y2 y1 b-a输出y1 b-a输出a = b = (1)=a+0.382*(b-(2)=a+0.618*(b-0.000000。 det(H) = h11 + h22 + i=1,2,NHijhij的代数 n(n2)Pn，那么有递推P+1 =(+n=2 2 阶行列式 P2=2，他需要做两次Pn=n!。由此可见，计算量跟 N 的阶 H = Hdet(H) = det(L)LUL1，所以他们的行列式可 分解的运算量是 O(N3)，所以总的运算量是O(N3) 输入
2、阶数N 和矩阵LU 输入阶数N 和矩阵LU输入阶数N HLU素的代21 3 M 1.000000 3.000000 21 3 M 1.000000 3.000000 det(M)-5M 0.8147000.0975000.1576000.1419000.9058000.2785000.9706000.4218000.1270000.5469000.9572000.9157000.9134000.9575000.4854000.7922000.6324000.9649000.8003000.959500det(M)- x = w x = w Nhy = N得到所有点的坐标(xij,yij)输入阶
3、数#defineEPS1e-6/doublefun(double double y#defineEPS1e-6/doublefun(double double y= returndoublefindMin(double(*f)(double),double double x1= x2= double z= findMin(fun,- returnx1=a+0.382*(b-x2=a+0.618*(b-y1=y2=a= b= returnMatrix,pMatrix N) pMatrixp =(pMatrix)Matrix,pMatrix N) pMatrixp =(pMatrix) p-N =
4、 p-data=(double *) returndoublegetMatrixItem(pMatrix i,j) if(i=p-N|j=p-f(stderr,Matrixindexisreturn returnp-datai*p-voidsetMatrixItem(pMatrix i,j,doublee) if(i=p-N|j=p-f(stderr,Matrixindexisreturn p-datai*p-N+j=voidMatrix(FILE*f,pMatrix f(f,%f,p- voiddeleteMatrix(pMatrix i,j,N=p-typedefstructdouble*
5、 double returnp- double returnp- return pMatrixsubMatrix =createMatrix(p-N)- 式 det+=det-= returndoublefastDetMatrix(pMatrixp) pMatrixL=createMatrix(p-N),U=L;/L,U共用同一空 doubledouble detMatrix(pMatrix p)if(i=j) /对应U矩阵tmp=getMatrixItem(p,i,j); if(i=j) /对应U矩阵tmp=getMatrixItem(p,i,j); tmp-if(i = j)det*=el
6、se /对应Ltmp=getMatrixItem(p,i,j); tmp-tmp/=getMatrixItem(U,j,j); returnvoid double pMatrixp = f(M= f(ndet(M)= * node, *net,a,b; /, RAMLPARAM LRESULTCALLBACK * node, *net,a,b; /, RAMLPARAM LRESULTCALLBACKWndProc(HWND, HINSTANCEhInstance = WNDCLASS HWND MSG wndclass.styleCS_HREDRAW|CS_VREDRAW|CS_NOCLOS
7、E改变窗 wndclass.lpfnWndProcWndProc窗口函数为 wndclass.cbClsExtra 0 wndclass.cbWndExtra 0 wndclass.hInstance = hInstance; wndclass.hIconLoadIcon(NULL,IDI_APPLICATION用箭头光 wndclass.hCursor= wndclass.hbrBackground = (HBRUSH)GetStockObject(WHITE_BRUSH背景为白色Name = NULL; ame 网格划分window /* 窗口类 if(!RegisterClass(&wndclass)return hwndCreateWindow(网格划分,网格划分 voidshowPic(HDC char dx= dy= dx= dy= x= 10 + y= 50 + LRESULTCALLBACKWndProc(HWND message, LPARAM switchcaseHDChdcSTRUCTps; hdchwnd,&ps获取要重画的窗口的设备hwnd,&ps结束要重画的窗口return * w,h) pi= *)malloc(sizeof * w,
温馨提示
- 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
评论
0/150
提交评论