原创力文档- 1、本文档共28页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
C程序设计实践结果解析
计算机程序设计实践报告
学生姓名： 李登峰
学生学号： 0306100218
所在班级： 生物工程1002班
任课教师： 杨英杰 刘新星
2011年6月24日
1. 某选矿车间某班取样化验数据如表1，利用矿量平衡和金属平衡可以列出含四个方程的线性方程组，请利用行线性变换编程求出各精矿及尾矿的理论产率
表1 各矿流的化验品位
--------------------------------
矿 产 化验品位%
率 ------------------
流 % Pb Zn S
--------------------------------
原 矿 100.00 3.08 11.92 14.58
铅精矿 x1 57.16 4.08 21.19
锌精矿 x2 0.79 56.56 31.71
硫精矿 x3 2.06 2.90 40.24
尾 矿 x4 0.25 0.73 4.93
--------------------------------
解题思路：利用矿量和三种组分的进出平衡可列出四个方程的线性方程组：
x1+x2+x3+x4=100
57.16x1+0.79x2+2.06x3+0.25x4=308
4.08x1+56.56x2+2.09x3+0.73x4=1192
21.19x1+31.71x2+40.24x3+4.93x4=1458
由于独立方程个数等于未知数个数属于唯一解的线性方程组。
令，X= x1 D=1 1 1 Y=100
x2 57.16 0.79 2.06 0.25 300
x3 4.08 56.56 2.09 0.73 1171
x4 21.19 31.71 40.24 4.93
X是未知数向量，D是系数矩阵，Y是观测值向量，则方程组可以表示为矩阵方程：
DX=Y
则可以求出X
通过对输入的矩阵的对角线上的数据和其下方的数据相比较，若对角线上的数的绝对值大于其所在列的数据，不交换；若对角线上的数据小于其所在列下方的数，则进行行交换。然后对对角线下方的数进行消元。使原矩阵变成上三角阵，接着用牛顿迭代法进行求解。即可求的线性方程组的解。
解：#includeiostream
#includeiomanip
using namespace std;
int main()
{ /*构建方程组输入系数矩阵*/
float a[4][5]={ {1,1,1,1,100},{57.16,0.79,2.06,0.25,308},
{4.08,56.56,2.90,0.73,1192},{21.19,31.71,40.24,4.93,1458}};
float t00=a[0][0],t10=a[1][0],t20=a[2][0],t30=a[3][0],t;
for(int m=0;m5;m++)
{
a[0][m]/=t00;
a[1][m]-=(t10*a[0][m]);
a[2][m]-=(t20*a[0][m]);
a[3][m]-=(t30*a[0][m]);
}
if(a[1][1]==0)
if(a[2][1]!=0)
for(int r=0;r5;r++)
{
t=a[1][r];
a[1][r]=a[2][r]; /*进行行变换求出逆矩阵*/
a[2][r]=t;
}
else
for(int s=0;s5;s++)
{
t=a[1][s];
a[1][s]=a[3][s];
a[2][3]=
文档评论（0）