代碼片段和文件信息
//?Homework_1.cpp?:?定義控制臺(tái)應(yīng)用程序的入口點(diǎn)。
//
#include?“stdafx.h“
#include?
#include?
using?namespace?std;
bool?matrixVectorMultiply(double?*adouble?*xdouble?*yint?N);
double?vectorInnerProduct(double?*xdouble?*yint?N);
bool?powerIterationMethod(double?*adouble?*udouble?*lambdaint?Nint?*time);
bool?inversePowerIterationMethod(double?*adouble?*udouble?*lambdaint?Nint?*time);
bool?theGaussSeidelIterationMethod(double?*adouble?*xdouble?*yint?N);
bool?directDoolittleMethod(double?*aodouble?*udouble?*yoint?N);
bool?detMatrix(double?*aodouble?*detint?N);
bool?directDoolittleSolveMethod(double?*adouble?*budouble?*bldouble?*cudouble?*cldouble?*udouble?*yoint?N);
bool?directDoolittleSplitMethod(double?*adouble?*budouble?*bldouble?*cudouble?*clint?N);
int?_tmain(int?argc?_TCHAR*?argv[])
{
double?a[501]am[501]u[501]lambda[3]lambdamaxlambdamindet;
double?b=0.16;
double?c=-0.064;
int?times=0;
for(int?i=0;i<501;i++)
{
a[i]=1.64-0.024*(i+1);
a[i]*=sin(0.2*(i+1));
a[i]-=0.64*exp(0.1/(i+1));
u[i]=1;
}
powerIterationMethod(au&lambda[0]501×);
for(int?i=0;i<501;i++) u[i]=1;
inversePowerIterationMethod(au&lambda[2]501×);
for(int?i=0;i<501;i++)
am[i]=a[i]+lambda[0];
for(int?i=0;i<501;i++) u[i]=1;
inversePowerIterationMethod(amu&lambda[1]501×);
lambdamin=min(lambda[1]-lambda[0]lambda[0]);
lambdamax=max(lambda[1]-lambda[0]lambda[0]);
cout.setf(ios::scientific);???????????//科學(xué)記數(shù)法?
cout.precision(12);???????????????????//精度為12,正常為6?
cout<<“λ1=“< cout<<“λ501=“< cout<<“λs=“< for(int?k=1;k<40;k++)
{
double?mu;
mu=lambdamin+k*(lambdamax-lambdamin)/40;
for(int?i=0;i<501;i++)
am[i]=a[i]-mu;
for(int?i=0;i<501;i++) u[i]=1;
inversePowerIterationMethod(amu&lambda[0]501×);
cout<<“k=“< }
cout<<“cond(A)=“< detMatrix(a&det501);
cout<<“det(A)=“< /*
cout<<“l(fā)ambda?1?or?lambda?501?=?“< cout<<“times?=?“< lambdamin=min(lambda[1]-lambda[0]lambda[0]);
lambdamax=max(lambda[1]-lambda[0]lambda[0]);
cout<<“l(fā)ambda?1?=?“< cout<<“l(fā)ambda?501?=?“< cout<<“l(fā)ambda?s?=?“<
for(int?k=0;k<40;k++)
{
int?mu;
mu=lambdamin+k*(lambdamax-lambdamin)/40;
for(int?i=0;i<501;i++)
am[i]=a[i]-mu;
inversePowerIterationMethod(amu&lambda[0]501×);
cout<<“k=“< }
cout<<“cond(A)=“< //printf(“k=%d?=>>>λi%d=%.13e\n“kkvalue_s);
//printf(“λs=%.13e\n“value_s);
double?det;
detMatrix(a&det501);
cout<<“det(A)=“< det=1;
for(int?i=0;i<501;i++)
det*=a[i];
cout<<“det(A)=“< */
return?0;
}
bool?detMatrix(double?*aodouble?*detint?N)
?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件??????53248??2014-10-31?18:35??Homework_1_1\debug\Homework_1.exe
?????文件?????655964??2014-10-31?18:35??Homework_1_1\debug\Homework_1.ilk
?????文件?????601088??2014-10-31?18:35??Homework_1_1\debug\Homework_1.pdb
?????文件???????6394??2014-10-31?18:35??Homework_1_1\Homework_1\Debug\BuildLog.htm
?????文件????????403??2014-10-27?10:45??Homework_1_1\Homework_1\Debug\Homework_1.exe.em
?????文件????????468??2014-10-27?10:45??Homework_1_1\Homework_1\Debug\Homework_1.exe.em
?????文件????????385??2014-10-31?18:35??Homework_1_1\Homework_1\Debug\Homework_1.exe.intermediate.manifest
?????文件??????69613??2014-10-31?18:35??Homework_1_1\Homework_1\Debug\Homework_1.obj
?????文件????1048576??2014-10-27?10:45??Homework_1_1\Homework_1\Debug\Homework_1.pch
?????文件?????????67??2014-10-31?18:35??Homework_1_1\Homework_1\Debug\mt.dep
?????文件??????10785??2014-10-27?10:45??Homework_1_1\Homework_1\Debug\stdafx.obj
?????文件?????232448??2014-10-31?18:35??Homework_1_1\Homework_1\Debug\vc80.idb
?????文件?????241664??2014-10-31?18:35??Homework_1_1\Homework_1\Debug\vc80.pdb
?????文件??????19598??2014-10-31?18:35??Homework_1_1\Homework_1\Homework_1.cpp
?????文件???????4560??2015-05-18?20:18??Homework_1_1\Homework_1\Homework_1.vcproj
?????文件???????4496??2014-10-27?10:24??Homework_1_1\Homework_1\Homework_1.vcproj.8.00.old
?????文件???????1427??2015-05-18?20:28??Homework_1_1\Homework_1\Homework_1.vcproj.woshidiannao-PC.woshidiannao.user
?????文件???????1427??2014-10-31?18:55??Homework_1_1\Homework_1\Homework_1.vcproj.WWW-A86CE6CB1B0.Administrator.user
?????文件????????968??2014-10-27?10:24??Homework_1_1\Homework_1\ReadMe.txt
?????文件????????215??2014-10-27?10:24??Homework_1_1\Homework_1\stdafx.cpp
?????文件????????276??2014-10-27?10:24??Homework_1_1\Homework_1\stdafx.h
?????文件????1289216??2015-05-18?20:28??Homework_1_1\Homework_1.ncb
?????文件????????896??2015-05-18?20:18??Homework_1_1\Homework_1.sln
?????文件????????895??2014-10-27?10:24??Homework_1_1\Homework_1.sln.old
????..A..H.?????17408??2015-05-18?20:28??Homework_1_1\Homework_1.suo
????..A..H.?????16384??2014-10-31?18:55??Homework_1_1\Homework_1.suo.old
????..A..H.?????16384??2014-10-31?18:55??Homework_1_1\Homework_1.v11.suo
?????文件???????2349??2015-05-18?20:18??Homework_1_1\UpgradeLog.xm
?????文件???????3348??2015-05-18?20:18??Homework_1_1\_UpgradeReport_Files\UpgradeReport.css
?????文件??????12579??2007-06-27?17:04??Homework_1_1\_UpgradeReport_Files\UpgradeReport.xslt
............此處省略10個(gè)文件信息
評(píng)論
共有 條評(píng)論
相關(guān)資源
- 數(shù)值分析算法程序設(shè)計(jì) C++實(shí)現(xiàn)
- 數(shù)值分析(Numerical Recipes)3rd Edition原
- Newton插值法求解一元高次方程的根
- 《數(shù)值分析》(Numerical Recipes)3rd E
- 算法經(jīng)典VC++數(shù)值分析
- 數(shù)值分析中的插值法包括拉格朗日n次
- 數(shù)值分析 牛頓插值法c++代碼
- 數(shù)值分析迭代法解非線性方程組c語(yǔ)言
- 北航數(shù)值分析大作業(yè)一
- Jacobi迭代法和Gauss-Seidel迭代法求解線
- 矩陣的三角分解c程序
- C語(yǔ)言 Newton迭代法解非線性方程組
- 數(shù)值分析實(shí)驗(yàn)之改進(jìn)的歐拉公式程序
- 數(shù)值分析三次樣條兩種邊界matlab代碼
- 數(shù)值分析牛頓迭代、三次樣條等五個(gè)
- 最小二乘法數(shù)值分析作業(yè)
- 數(shù)值分析c語(yǔ)言追趕法解方程