牛顿迭代法

Author: computeone

BinaryDivide.cpp

@@ -34,8 +34,9 @@ double binaryDivide(double a,double b,double tol){
 		}
 
 	}
-	return result;
 	cout << "二分法求根结束：" << endl;
+	return result;
+	
 }
 
 double binaryDivide(double a, double b, double tol, double(*function)(double)){
@@ -66,6 +67,7 @@ double binaryDivide(double a, double b, double tol, double(*function)(double)){
 		}
 	}
 
-	return result;
 	cout << "二分法求根结束:" << endl;
+	return result;
+	
 }

CNumericalCalculation.vcxproj

@@ -77,9 +77,11 @@
   </ItemDefinitionGroup>
   <ItemGroup>
     <ClInclude Include="BinaryDivide.h" />
+    <ClInclude Include="FixedIterator.h" />
   </ItemGroup>
   <ItemGroup>
     <ClCompile Include="BinaryDivide.cpp" />
+    <ClCompile Include="FixedIterator.cpp" />
     <ClCompile Include="main.cpp" />
   </ItemGroup>
   <Import Project="$(VCTargetsPath)\Microsoft.Cpp.targets" />

CNumericalCalculation.vcxproj.filters

@@ -18,6 +18,9 @@
     <ClInclude Include="BinaryDivide.h">
       <Filter>头文件</Filter>
     </ClInclude>
+    <ClInclude Include="FixedIterator.h">
+      <Filter>头文件</Filter>
+    </ClInclude>
   </ItemGroup>
   <ItemGroup>
     <ClCompile Include="BinaryDivide.cpp">
@@ -26,5 +29,8 @@
     <ClCompile Include="main.cpp">
       <Filter>源文件</Filter>
     </ClCompile>
+    <ClCompile Include="FixedIterator.cpp">
+      <Filter>源文件</Filter>
+    </ClCompile>
   </ItemGroup>
 </Project>

Debug/BinaryDivide.obj

@@ -1,15 +1,19 @@
-生成启动时间为 4/4/2014 11:27:11 PM。
+生成启动时间为 4/5/2014 6:09:47 PM。
      1>项目“D:\C++\WorkSpace\SunSystem\CNumericalCalculation\CNumericalCalculation.vcxproj”在节点 2 上(Build 个目标)。
      1>ClCompile:
-         D:\Program Files (x86)\Microsoft Visual Studio 12.0\VC\bin\CL.exe /c /ZI /nologo /W3 /WX- /Od /Oy- /D WIN32 /D _DEBUG /D _CONSOLE /D _LIB /D _UNICODE /D UNICODE /Gm /EHsc /RTC1 /MDd /GS /fp:precise /Zc:wchar_t /Zc:forScope /Fo"Debug\\" /Fd"Debug\vc120.pdb" /Gd /TP /analyze- /errorReport:prompt main.cpp
+         D:\Program Files (x86)\Microsoft Visual Studio 12.0\VC\bin\CL.exe /c /ZI /nologo /W3 /WX- /Od /Oy- /D WIN32 /D _DEBUG /D _CONSOLE /D _LIB /D _UNICODE /D UNICODE /Gm /EHsc /RTC1 /MDd /GS /fp:precise /Zc:wchar_t /Zc:forScope /Fo"Debug\\" /Fd"Debug\vc120.pdb" /Gd /TP /analyze- /errorReport:prompt FixedIterator.cpp main.cpp
          main.cpp
-     1>d:\c++\workspace\sunsystem\cnumericalcalculation\main.cpp(7): warning C4244: “参数”: 从“double”转换到“float”，可能丢失数据
+     1>d:\c++\workspace\sunsystem\cnumericalcalculation\main.cpp(8): warning C4244: “参数”: 从“double”转换到“float”，可能丢失数据
+         FixedIterator.cpp
+     1>d:\c++\workspace\sunsystem\cnumericalcalculation\fixediterator.cpp(14): warning C4244: “参数”: 从“double”转换到“float”，可能丢失数据
+         正在生成代码...
        Link:
          D:\Program Files (x86)\Microsoft Visual Studio 12.0\VC\bin\link.exe /ERRORREPORT:PROMPT /OUT:"D:\C++\WorkSpace\SunSystem\Debug\CNumericalCalculation.exe" /INCREMENTAL /NOLOGO kernel32.lib user32.lib gdi32.lib winspool.lib comdlg32.lib advapi32.lib shell32.lib ole32.lib oleaut32.lib uuid.lib odbc32.lib odbccp32.lib /MANIFEST /MANIFESTUAC:"level='asInvoker' uiAccess='false'" /manifest:embed /DEBUG /PDB:"D:\C++\WorkSpace\SunSystem\Debug\CNumericalCalculation.pdb" /SUBSYSTEM:CONSOLE /TLBID:1 /DYNAMICBASE /NXCOMPAT /IMPLIB:"D:\C++\WorkSpace\SunSystem\Debug\CNumericalCalculation.lib" /MACHINE:X86 Debug\BinaryDivide.obj
+         Debug\FixedIterator.obj
          Debug\main.obj
          CNumericalCalculation.vcxproj -> D:\C++\WorkSpace\SunSystem\Debug\CNumericalCalculation.exe
      1>已完成生成项目“D:\C++\WorkSpace\SunSystem\CNumericalCalculation\CNumericalCalculation.vcxproj”(Build 个目标)的操作。
 
 生成成功。
 
-已用时间 00:00:01.23
+已用时间 00:00:02.89

Debug/CNumeric.7D3823D4.tlog/CL.read.1.tlog

@@ -0,0 +1,40 @@
+#include "FixedIterator.h"
+#include <iostream>
+#include <iomanip>
+
+using namespace std;
+
+
+double fixedIterator(double p0, double tol,int N){
+
+	int i = 0;
+	double p;
+	cout << "不动点迭代开始:" << endl;
+	while (i < (N+1)){
+		p = p0 - (powf(p0, 3) + 4 * powf(p0, 2) - 10.0) / (3 * powf(p0, 2) + 8 * p0);
+		cout << setprecision(16)<<"p:" << p << endl << endl;
+		if (fabs(p - p0) < tol){
+			cout << "不动点迭代结束：" << endl;
+			return p;
+		}
+		p0 = p;
+	}
+	cout << "不动点迭代结束:迭代失败!" << endl;
+	return NULL;
+}
+
+double fixedIterator(double p0, double tol,int N, double(*function)(double p)){
+	int i = 0;
+	double p;
+	cout << "不动点迭代开始:" << endl;
+	while (i < (N + 1)){
+		p = function(p0);
+		cout << setprecision(16)<<"p:" << endl << endl;
+		if (fabs(p - p0) < tol){
+			cout << "不动点迭代结束：" << endl;
+		}
+		p0 = p;
+	}
+	cout << "不动点迭代结束：迭代失败！" << endl;
+	return NULL;
+}

Debug/CNumeric.7D3823D4.tlog/CL.write.1.tlog

@@ -0,0 +1,8 @@
+#include <iostream>
+
+using namespace std;
+
+
+//˛ťśŻľăľü´úˇ¨p=g(p)
+double fixedIterator(double p0, double tol,int N);
+double fixedIterator(double p0, double tol,int N, double(*function)(double p));

Debug/CNumeric.7D3823D4.tlog/cl.command.1.tlog

@@ -1,4 +1,5 @@
 #include "BinaryDivide.h"
+#include "FixedIterator.h"
 #include <iostream>
 
 using namespace std;
@@ -7,6 +8,12 @@ double test(double p){
 	return powf(p, 3) + 4 * powf(p, 2) - 10.0;
 };
 
+void testBinaryDivide(){
+	cout << "¸ùΪ£º"<<binaryDivide(1.0, 2.0,0.000002,&test) << endl;
+}
+void testFixedIterator(){
+	cout << "¸ùΪ:" << fixedIterator(1.0, 0.000001, 10) << endl;
+}
 int main(){
-	cout << "¸ùΪ£º"<<binaryDivide(1.0, 2.0,0.000001,&test) << endl;
+	testFixedIterator();
 }