#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
eeworm.com VIP专区 单片机源码系列 3资源包含以下内容:1. N97.rar2. TMS320C2000.rar3. wince.rar4. C8051F.rar5. ARM.rar6. nichelite.for.lpc.zip7. 12864.rar8. 神龙卡SDK (DirectShow版).zip9. READ_STC12C5A60S2.rar10. DS18B20.rar11. PowerLogic5.0.1.rar12. 51.rar13. 2410.rar14. UCOS2.52.rar15. arm9.rar16. yitaiwang.rar17. JLDLED.rar18. 51单片机方波发生器.rar19. efsl-0.3.6.zip20. MyGui.rar21. Zigbee2007.rar22. ZIGBEE.rar23. LCD12864.rar24. uCOS-II.rar25. uc-gui.rar26. wifi-of-88w8686.rar27. S3C2410_TEST.rar28. H-JTAG-V0.6.1.zip29. ZigBee_xbee.rar30. S7200Sysm.rar31. Led_font_extraction.zip32. IAR_example_EasyARM8962.zip33. siemensPLCjiaocheng.rar34. MFRC500antenna.rar35. LCD1602.rar36. ucos.rar37. s3c2410.rar38. 20086238534741.rar39. ISE7.1lesson.rar40. zlg_fs.rar41. DM9000.rar42. shuzicunchu.rar43. lgui_0.3.0.rar44. 单片机的40个实验.pdf45. STM32_fatfs_shell_SDHC.rar46. CC2430_chinese_pdf.rar47. cypressusb2.0.zip48. PIC.rar49. UCOS-2.rar50. pmc.rar51. boot.rar52. mp3.rar53. STdemo.rar54. ARM9.rar55. ARM7Proteus.rar56. 240128.rar57. AVRkaifashili.rar58. motor_control.zip59. yaoyaobang.rar60. wince_serialport.rar61. Document.rar62. zlggui.rar63. freescale.rar64. NIOS_UART_TEST.rar65. LPC2290.rar66. 硬件cpu&rom课程设计.rar67. ucfs.rar68. TL431.rar69. 单片机读写usb、sd卡技术参考资料.rar70. TCP_IPprotocol.rar71. TS8001K.rar72. nucleus_arm.rar73. 12864串行通信驱动.rar74. stc12c5608ad.zip75. ATmega16读写SD卡程序.rar76. Keil Cx51 V7.0单片机高级语言编程与uVision2应用实践.rar77. freescale单片机应用实例.rar78. 周立功 开发包 ARM9.rar79. 智能小车.rar80. 串行下载线的原理图.pdf81. C51入门.rar82. C51入门2.rar83. C51入门3.rar84. C51入门4.rar85. C51入门5.rar86. C51入门6.rar87. C51入门7.rar88. 单片机 漏电检测原理.doc89. nrf24z1 代码.rar90. 舞蹈机器人设计与制作.pdf91. C8051系列单片机开发与C语言编程!.rar92. PIC单片机C语言应用例程.rar93. 基于89c52的12864液晶程序.txt.rar94. 温度控制系统——毕业设计.rar95. 力天电子LPC210X例程代码.zip96. 新型楼道声控灯(毕业设计).rar97. stm32数码相框.rar98. 宏晶单片机的AD演示程序.rar99. 在VB下PC机与MCS.doc.doc100. 单片机技术应用教程.pdf
上传时间: 2013-05-15
上传用户:eeworm
第一章 有关数论的算法 1.1最大公约数与最小公倍数 1.2有关素数的算法 1.3方程ax+by=c的整数解及应用 1.4 求a^b mod n 第二章 高精度计算 2.1高精度加法 2.2高精度减法 2.3高精度乘法 2.4 高精度除法 练习 第三章 排列与组合 3.1加法原理与乘法原理 练习 3. 2 排列与组合的概念与计算公式 练习 3.3排列与组合的产生算法 练习 第四章 计算几何 4.1 基础知识 4.2 线段的相交判断 4.3寻找凸包算法 练习 第五章 其它数学知识及算法 5.1 鸽巢原理 5.2 容斥原理及应用 5.3 常见递推关系及应用
上传时间: 2016-01-05
上传用户:frank1234
专辑类-电子基础类专辑-153册-2.20G 21世纪大学新型参考教材系列-集成电路B-荒井-159页-2.8M.pdf
上传时间: 2013-05-16
上传用户:pkkkkp
里面有相应的hpunix(HP-UX hpl1000 B.11.00 U 9000/800 (tb)),linux(Red Hat Linux release 9 Kernel 2.4.20-8),windows的头文件、库文件,还有相应的demo程序
标签: release hpunix Kernel HP-UX
上传时间: 2015-01-06
上传用户:cursor
2次B样条曲线算法
标签: 算法
上传时间: 2014-01-27
上传用户:来茴
注塑机生产过程的微机控制系统设计本系统设计是通过8255A并行端口A、B传输数据控制12个LED发光二级管的亮灭来模拟控制两台注塑机的生产过程。其中端口A和B分别对应1号和2号注塑机,每一个LED代表一个工序。编程语言采用BORLAND C++语言。
上传时间: 2013-12-25
上传用户:叶山豪
使用的是API编程,可格式化、校验和读写特殊扇区。可用作Windows下的磁盘加密。本函数还有以下两个缺点以待改进: 1.本函数还只能读能读 A: 和 B:,即只能对软盘操作 2.不能改变磁盘扇区大小,只能是标准的 512 个字节。 参数说明: command 操作: 0 重置磁盘 2 读扇区 3 写扇区 4 校验磁道 5 格式化磁道 8 得到设备参数 (int 1EH) drive 驱动器 A:=0 B:=1 head 磁头号,范围 0 - 1 track 磁道号,范围 0 - 84 ( 80 - 84 为特殊磁道,通常用来加密 ) sector 扇区号,范围 0 - 255 ( 19 - 255 为非标准扇区编号,通常用来加密) nsectors 每次读或写的扇区数,不能超出每磁道的最大扇区数 buffer 数据写入或读出的缓冲区,大小为 512 个字节 返回值 ( 同 Int 13H ): 0x0 成功 0x1 无效的命令 0x3 磁盘被写保护 0x4 扇区没有找到 0xa 发现坏扇区 0x80 磁盘没有准备好
上传时间: 2013-12-05
上传用户:moerwang
书 名:Programming Windows程式开发设计指南 出版日期:2000/6/2 书 号:957-8239-73-4 I S B N:957-8239-73-4 原 作 者:Charles Petzold 译 者:余孟学
标签: 8239 Programming 957 Windows
上传时间: 2015-04-26
上传用户:xinyuzhiqiwuwu
CAN-bus 规范V2.0 版本 这本CAN 规范技术规范由两部分组成: • A 部分:CAN 的报文格式说明(按CAN1.2 规范定义)。 • B 部分:标准格式和扩展格式的说明。
上传时间: 2015-06-04
上传用户:天诚24