#include "iostream" using namespace std; class Matrix { private: double** A; //矩阵A double *b; //向量b public: int size; Matrix(int ); ~Matrix(); friend double* Dooli(Matrix& ); void Input(); void Disp(); }; Matrix::Matrix(int x) { size=x; //为向量b分配空间并初始化为0 b=new double [x]; for(int j=0;j<x;j++) b[j]=0; //为向量A分配空间并初始化为0 A=new double* [x]; for(int i=0;i<x;i++) A[i]=new double [x]; for(int m=0;m<x;m++) for(int n=0;n<x;n++) A[m][n]=0; } Matrix::~Matrix() { cout<<"正在析构中~~~~"<<endl; delete b; for(int i=0;i<size;i++) delete A[i]; delete A; } void Matrix::Disp() { for(int i=0;i<size;i++) { for(int j=0;j<size;j++) cout<<A[i][j]<<" "; cout<<endl; } } void Matrix::Input() { cout<<"请输入A:"<<endl; for(int i=0;i<size;i++) for(int j=0;j<size;j++){ cout<<"第"<<i+1<<"行"<<"第"<<j+1<<"列:"<<endl; cin>>A[i][j]; } cout<<"请输入b:"<<endl; for(int j=0;j<size;j++){ cout<<"第"<<j+1<<"个:"<<endl; cin>>b[j]; } } double* Dooli(Matrix& A) { double *Xn=new double [A.size]; Matrix L(A.size),U(A.size); //分别求得U，L的第一行与第一列 for(int i=0;i<A.size;i++) U.A[0][i]=A.A[0][i]; for(int j=1;j<A.size;j++) L.A[j][0]=A.A[j][0]/U.A[0][0]; //分别求得U，L的第r行，第r列 double temp1=0,temp2=0; for(int r=1;r<A.size;r++){ //U for(int i=r;i<A.size;i++){ for(int k=0;k<r-1;k++) temp1=temp1+L.A[r][k]*U.A[k][i]; U.A[r][i]=A.A[r][i]-temp1; } //L for(int i=r+1;i<A.size;i++){ for(int k=0;k<r-1;k++) temp2=temp2+L.A[i][k]*U.A[k][r]; L.A[i][r]=(A.A[i][r]-temp2)/U.A[r][r]; } } cout<<"计算U得:"<<endl; U.Disp(); cout<<"计算L的:"<<endl; L.Disp(); double *Y=new double [A.size]; Y[0]=A.b[0]; for(int i=1;i<A.size;i++ ){ double temp3=0; for(int k=0;k<i-1;k++) temp3=temp3+L.A[i][k]*Y[k]; Y[i]=A.b[i]-temp3; } Xn[A.size-1]=Y[A.size-1]/U.A[A.size-1][A.size-1]; for(int i=A.size-1;i>=0;i--){ double temp4=0; for(int k=i+1;k<A.size;k++) temp4=temp4+U.A[i][k]*Xn[k]; Xn[i]=(Y[i]-temp4)/U.A[i][i]; } return Xn; } int main() { Matrix B(4); B.Input(); double *X; X=Dooli(B); cout<<"~~~~解得:"<<endl; for(int i=0;i<B.size;i++) cout<<"X["<<i<<"]:"<<X[i]<<" "; cout<<endl<<"呵呵呵呵呵"; return 0; }
标签： 道理特分解法
上传时间： 2018-05-20
上传用户：Aa123456789
FC162是一款低功耗，高速，高噪声容限，EPROM/ROM基于8位CMOS工艺制造的单片机，采用RISC指令集，共有42条指令， 除分支指令为两个周期指令以外其余为单周期指令。这种易用、易记的指令集大大缩短了开发时间。 FC162包含了上电复位(Power-on Reset POR)，掉电复位(Brown-out Reset BOR)， 上电复位计数器（Power-up Reset Timer PWRT），振荡启动计数器 (Oscillator Start-up Timer OST)， 看门狗定时器(Watchdog Timer)， EPROM/ROM， SRAM，双向三 态I/O口，（可以设置为上拉/下拉）， 省电睡眠模式， 一个带8位预置器的8位定时/计数器，独立中断，睡眠唤醒模式和可靠的代码保 护，有两个振荡源可供用户配置选择，包含省电振荡源和低功耗振荡器。 FC162可访问256×13的程序存储空间。 FC162能直接或间接访问寄存器以及数据存储区，所有的特殊功能寄存器分布在数据存储区同时包含特定的程序指针。
标签： fc162
上传时间： 2021-11-13
上传用户：qingfengchizhu
B＋树实现源码，JAVA语言，数据库索引技术
上传时间： 2013-12-22
上传用户：xc216
b+樹源碼，b+樹結構，刪除，插入，等值搜索，範圍搜索等功能
标签：
上传时间： 2015-03-23
上传用户：tianjinfan
无线数传电台串口通信C程序源码 电台型号为深圳友迅达公司的F201/B无线数传电台
上传时间： 2014-01-24
上传用户：王小奇
非均匀有理b样条曲线nurbs的C＋＋源码库。 3.0.11版
上传时间： 2015-04-24
上传用户：Miyuki
基于b/s的酒店管理系统asp源码，大家免费下载，功能还行，就是有一些小bug
上传时间： 2014-01-02
上传用户：qlpqlq
生产排程计划的一份源码.JAVA.B/S结构的。希望对一部分人有帮助
上传时间： 2015-08-30
上传用户：小草123
这是Kenneth C.Loudden所著<<编译原理与实践>>附录B中源码.但原书公布的下载地址已失效.故上传
上传时间： 2016-01-23
上传用户：x4587
比较牛b的一个计算器程序,源码在word文档里面,
上传时间： 2016-02-14
上传用户：z1191176801
