哪位高手看一下，这是输出什么？

skycncomp 2004-05-08 06:39:53
哪位帮忙看看，这个程序想要输出什么，这个程序是什么意思？
被弄得太晕了。变理乱七八糟的。都不知道输出什么。请指点，给程序加上注释.
谢谢了。
#include "aband.c"
#include "stdio.h"
main()
{ int i;
static double b[8][5]={ {3.0,-4.0,1.0,0.0,0.0},
{-2.0,-5.0,6.0,1.0,0.0},{1.0,3.0,-1.0,2.0,-3.0},
{2.0,5.0,-5.0,6.0,-1.0},{-3.0,1.0,-1.0,2.0,-5.0},
{6.0,1.0,-3.0,2.0,-9.0},{-4.0,1.0,-1.0,2.0,0.0},
{5.0,1.0,-7.0,0.0,0.0}};
static double d[8][3]={ {13.0,29.0,-13.0},
{-6.0,17.0,-21.0},{-31.0,-6.0,4.0},{64.0,3.0,16.0},
{-20.0,1.0,-5.0},{-22.0,-41.0,56.0},{-29.0,10.0,-21.0},
{7.0,-24.0,20.0}};
if (aband(b,d,8,2,5,3)>0)
for (i=0;i<=7;i++)
printf("x(%d)=%13.7e, %13.7e, %13.7e\n",
i,d[i][0],d[i][1],d[i][2]);
}

#include "math.h"
#include "stdio.h"
int aband(b,d,n,l,il,m)
int n,l,il,m;
double b[],d[];
{ int ls,k,i,j,is,u,v;
double p,t;

if (il!=(2*l+1))
{ printf("fail\n"); return(-2);}
ls=l;

for (k=0;k<=n-2;k++)
{
p=0.0;
for (i=k;i<=ls;i++)
{
t=fabs(b[i*il]);
if (t>p) {p=t; is=i;}
}
if (p+1.0==1.0)
{
printf("fail\n"); return(0);
}

for (j=0;j<=m-1;j++)
{
u=k*m+j; v=is*m+j;
t=d[u]; d[u]=d[v]; d[v]=t;
}

for (j=0;j<=il-1;j++)
{ u=k*il+j; v=is*il+j;
t=b[u]; b[u]=b[v]; b[v]=t;
}
for (j=0;j<=m-1;j++)
{ u=k*m+j; d[u]=d[u]/b[k*il];}
for (j=1;j<=il-1;j++)
{ u=k*il+j; b[u]=b[u]/b[k*il];}
for (i=k+1;i<=ls;i++)
{ t=b[i*il];
for (j=0;j<=m-1;j++)
{ u=i*m+j; v=k*m+j;
d[u]=d[u]-t*d[v];
}
for (j=1;j<=il-1;j++)
{ u=i*il+j; v=k*il+j;
b[u-1]=b[u]-t*b[v];
}
u=i*il+il-1; b[u]=0.0;
}
if (ls!=(n-1)) ls=ls+1;
}
p=b[(n-1)*il];
if (fabs(p)+1.0==1.0)
{ printf("fail\n"); return(0);}
for (j=0;j<=m-1;j++)
{ u=(n-1)*m+j; d[u]=d[u]/p;}
ls=1;
for (i=n-2;i>=0;i--)
{ for (k=0;k<=m-1;k++)
{ u=i*m+k;
for (j=1;j<=ls;j++)
{ v=i*il+j; is=(i+j)*m+k;
d[u]=d[u]-b[v]*d[is];
}
}
if (ls!=(il-1)) ls=ls+1;
}
return(2);
}

...全文
80 21 打赏 收藏 转发到动态 举报
AI 作业
写回复
用AI写文章
21 条回复
切换为时间正序
请发表友善的回复…
发表回复
improgrammer 2004-05-11
  • 打赏
  • 举报
 回复
知道“矩阵”、“矩阵乘法”等概念吗？

系数矩阵B：
{
{ 3.0,-4.0, 1.0, 0.0, 0.0},
{-2.0,-5.0, 6.0, 1.0, 0.0},
{ 1.0, 3.0,-1.0, 2.0,-3.0},
{ 2.0, 5.0,-5.0, 6.0,-1.0},
{-3.0, 1.0,-1.0, 2.0,-5.0},
{ 6.0, 1.0,-3.0, 2.0,-9.0},
{-4.0, 1.0,-1.0, 2.0, 0.0},
{ 5.0, 1.0,-7.0, 0.0, 0.0}
};


得到的解X：
+3.000 +5.000 +0.000
-1.000 -3.000 +3.000
+0.000 +2.000 -1.000
-5.000 -0.000 +0.000
+7.000 -0.000 +2.000
+1.000 +1.000 -3.000
+2.000 -1.000 +0.000
+0.000 +4.000 -5.000

方程右边（矩阵D）：
{
{ 13.0, 29.0,-13.0},
{ -6.0, 17.0,-21.0},
{-31.0, -6.0, 4.0},
{ 64.0, 3.0, 16.0},
{-20.0, 1.0, -5.0},
{-22.0,-41.0, 56.0},
{-29.0, 10.0,-21.0},
{ 7.0,-24.0, 20.0}
};

你做一下矩阵乘法，看是不是B*X=D。

improgrammer 2004-05-11
  • 打赏
  • 举报
 回复
告诉你了还不知道，看书去吧。
skycncomp 2004-05-11
  • 打赏
  • 举报
 回复
我顶
improgrammer 2004-05-10
  • 打赏
  • 举报
 回复
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <assert.h>

//判断p是否充分接近零
int near_zero(double p)
{
//——只要浮点数p与1.0对齐后没有有效数字，就判p充分接近零
return p+1.0==1.0;
}

//打印两矩阵
void my_print(double b[],double d[],const int N,const int IL,const int M)
{
int i,j;
for (i=0;i<N;i++)
{
printf("row[%d]",i);
for(j=0;j<IL;++j)
{
printf(" %+7.3f",b[i*IL+j]);
}
printf(" |");
for(j=0;j<M;++j)
{
printf(" %+7.3f",d[i*M+j]);
}
printf("\n");
}
printf("\n");
}

//向量点积
double vec_dm(const double *v1,const double *v2,int d1,int d2,int n)
{
double r=0.0;
while((n--) > 0)
{
r += (*v1)*(*v2);
v1+=d1; v2+=d2;
}
return r;
}

//v1 <--> v2
void vec_swap(double *v1,double *v2,int d1,int d2,int n)
{
double t;
while((n--) > 0)
{
t=*v1;*v1=*v2;*v2=t;
v1+=d1;v2+=d2;
}
}

// vec *= scale
void vec_scale(double *vec,int d,int n,double scale)
{
while((n--) > 0)
{
*vec *= scale;
vec+=d;
}
}

// v1 += v2 * scale
void vec_scale_add(double *v1,const double *v2,int d1,int d2,int n,double scale)
{
while((n--) > 0)
{
(*v1) += (*v2) * scale;
v1+=d1;v2+=d2;
}
}

//求向量中绝对值最大元素之索引
int vec_max_index(const double *v,const int d,const int n)
{
int i,r=0;
double t,p=0.0;
for(i=0;i<n;++i)
{
t=fabs(*v);
if(t>p)
{
r=i;
p=t;
}
v+=d;
}
return r;
}

//向量左移一个元素，最左元素被丢弃，最右元素写为0
void vec_lshift(double *v,int d,int n)
{
while((n--)>1)
{
*v = *(v+d);
v+=d;
}
*v=0.0;
}

//哪位高手看一下，这是输出什么？
//——应该是解方程
int aband(double b[],double d[],const int N,const int IL,const int M)
{
const int L=(IL+1)/2;
int i,j,is,_ls;
double p;
//第一步：好象是做一个消元法，完成后b数组的第一列为主对角线元素
for (i=0;i<N-1;i++)
{
_ls=min(i+L,N);
//从b[i][0],b[i+1][0]..b[_ls-1][0]找最大元素[绝对值最大]
//——最大元素的索引保存在is中
is=i+vec_max_index(b+i*IL,IL,_ls-i);
p=b[is*IL];
assert(!near_zero(p));//最大元素应不为0
//对调第i,is行
vec_swap(d+i*M,d+is*M,1,1,M);
vec_swap(b+i*IL,b+is*IL,1,1,IL);
//d[N][M]第i行除以b[N][IL]的第i行排头元素
vec_scale(d+i*M,1,M,1/p);
//b[N][M]第i行（排头除外）除以排头元素
vec_scale(b+i*IL+1,1,IL-1,1/p);
//
for (j=i+1;j<_ls;j++)
{
p=b[j*IL];
vec_scale_add(d+j*M,d+i*M,1,1,M,-p);
vec_scale_add(b+j*IL+1,b+i*IL+1,1,1,IL-1,-p);
vec_lshift(b+j*IL,1,IL);
}
}
//d[N-1][0..M-1] /= p;
p=b[(N-1)*IL];
assert(!near_zero(fabs(p)));
vec_scale(d+(N-1)*M,1,M,1/p);
//第二步：某种卷积运算
for (i = N - 2;i >= 0;i--)
{
const int R = min(N-i-1,IL-1);
for (j = 0;j < M;j++)
{
d[i*M+j] -= vec_dm(b+i*IL+1,d+(i+1)*M+j,1,M,R);
}
}
return(2);
}

void work()
{
static double b[8][5]=
{
{ 3.0,-4.0, 1.0, 0.0, 0.0},
{-2.0,-5.0, 6.0, 1.0, 0.0},
{ 1.0, 3.0,-1.0, 2.0,-3.0},
{ 2.0, 5.0,-5.0, 6.0,-1.0},
{-3.0, 1.0,-1.0, 2.0,-5.0},
{ 6.0, 1.0,-3.0, 2.0,-9.0},
{-4.0, 1.0,-1.0, 2.0, 0.0},
{ 5.0, 1.0,-7.0, 0.0, 0.0}
};
static double d[8][3]=
{
{ 13.0, 29.0,-13.0},
{ -6.0, 17.0,-21.0},
{-31.0, -6.0, 4.0},
{ 64.0, 3.0, 16.0},
{-20.0, 1.0, -5.0},
{-22.0,-41.0, 56.0},
{-29.0, 10.0,-21.0},
{ 7.0,-24.0, 20.0}
};
if (aband(b,d,8,5,3)>0)
{
my_print(b,d,8,5,3);
}
}
skycncomp 2004-05-10
  • 打赏
  • 举报
 回复
上面的是什么呀？
我实在不知道。
skycncomp 2004-05-10
  • 打赏
  • 举报
 回复
上面是两个程序，主函数调用的。
我只是没有分天而以。
自己分一下就OK了。
improgrammer 2004-05-10
  • 打赏
  • 举报
 回复

用列主消元法解线性方程组。

B*X=D
B是N*N方阵
X是N*M矩阵
D是N*M矩阵

B是对角方阵（宽度是IL），所以b[N][IL]采用了压缩存储法。

skycncomp 2004-05-09
  • 打赏
  • 举报
 回复
我知道结果。不过不知道是什么意思呀。
上面输出的什么呢？用科学计数法表示数字吗？
我不知道，哪位给解释一下。
holyfair 2004-05-09
  • 打赏
  • 举报
 回复
运行结果你运行一下就知道了。意思你自己都不懂你看这个程序干什么？
x(0)=3.0000000e+000, 5.0000000e+000, -1.2767565e-015
x(1)=-1.0000000e+000, -3.0000000e+000, 3.0000000e+000
x(2)=8.3266727e-016, 2.0000000e+000, -1.0000000e+000
x(3)=-5.0000000e+000, -2.6645353e-015, -1.7763568e-01
x(4)=7.0000000e+000, -8.8817842e-016, 2.0000000e+000
x(5)=1.0000000e+000, 1.0000000e+000, -3.0000000e+000
x(6)=2.0000000e+000, -1.0000000e+000, -8.8817842e-016
x(7)=1.5199238e-016, 4.0000000e+000, -5.0000000e+000


skycncomp 2004-05-09
  • 打赏
  • 举报
 回复
在上边了。
一看你就没好好看程序。：……》
上边的程序是两个。看好了。
易点互联 2004-05-09
  • 打赏
  • 举报
 回复
aband.c原文件在那里
liangenbo 2004-05-09
  • 打赏
  • 举报
 回复
在开头你用#include "aband.c"是包含头文件，而你后面又用的是函数的形式啊
我想应该会出错的哈
不知道我对不对哈
skycncomp 2004-05-09
  • 打赏
  • 举报
 回复
我晕，我顶。
非要顶出个所以然来.
improgrammer 2004-05-09
  • 打赏
  • 举报
 回复
还没看完，估计是某种解向量值方程组的算法。
或者与线性规划有关。
skycncomp 2004-05-09
  • 打赏
  • 举报
 回复
高手进来看看呀。我看得头晕了也不明白是什么东西。
运行结果不知道是什么。
julyclyde 2004-05-09
  • 打赏
  • 举报
 回复
不要到处找一些不知道干什么的程序来看
对提高水平没好处
cancer001 2004-05-08
  • 打赏
  • 举报
 回复
帮你顶！
wangnewton 2004-05-08
  • 打赏
  • 举报
 回复
帮你顶

ps：一个注释都没有
很难有看头，一行一行看下来，会死人的
期待高手出马
sten 2004-05-08
  • 打赏
  • 举报
 回复
运行一下不就得了吗？
tianyaht 2004-05-08
  • 打赏
  • 举报
 回复
确实，没有注释
真的是很难看的……
我是没有耐心，只有希望高手了
加载更多回复（1）
c语言字母意义,%C是什么意思？ c语言中？和：是什么意思
导航：网站首页 > %C是什么意思？ c语言中？和：是什么意思 %C是什么意思？ c语言中？和：是什么意思 相关问题: 匿名网友: %c单个字符输出的意思。 %s是输出字符串。 %d是输出整型。 %f是输出整型。 这是在scanf、printf这样的函数中，将参数类型与输入/输出句子中的相应位置相对应所使用的转义标志符。 其中常用的有%c--对应字符型char%s--对应字符串型char*(char[])%d--对应整形int%ld--对应长整形long%x--对应十六进制整数，大
为什么C语言输出结果总是1,请高手解答一下，为什么我打印出来的，总是两个相同的数据...
请高手解答一下，为什么我打印出来的，总是两个相同的数据为什我输入一个数据后，再打印，出现两个一模一样的数所，比如，我输入，学号：1001 姓名 aaa,打印出为居然是两个1001 aaa,1001 aaa这是为什么，我看讲课上老师打出来的是一个啊，同样是保存再加载。代码如下：#define_CRT_SECURE_NO_WARNINGS#include#include#includevoid...
【194期】HashSet 是无序的吗？为什么输出的结果有序？
点击上方“Java精选”，选择“设为星标”别问别人为什么，多问自己凭什么！下方有惊喜，留言必回，有问必答！每天08:15更新文章，每天进步一点点...什么是集合？集合简称集，是数学中一...
某位高手对高输入阻抗和低输出阻抗的解释
一般我们常耳闻的说法是：扩大机的输入阻抗是愈高愈好，而输出阻抗是愈低愈好。为什么呢？因为输入阻抗高了，从讯号源来的讯号功率强度就可以不必那么大。这么说也许还有读者不甚了解，让我们再回想一下欧姆定律；假设讯源输出不甚了解，让我们再回想一下欧姆定律；假设讯源输出一个固定电压，传送往下一级，如果这一级的输入阻抗高，是不是由讯源所提供的讯号电流就可以降低？如果输入阻抗非常非常的高，则几乎不会消耗讯号电流(
#NAME?_#NAME?:EXCEL表中出现#NAME是什么意思？
EXCEL表中出现#NAME是什么意思？这是因为公式中A、B和C+是文本，1、如图输入题主的公式。会弹出所键入的公式含有错误，会显示C这里有错，不止C这里有错，前面的A和B也有错，3、因为公式EXCEL表中出现#NAME是什么意思？这是因为公式中A、B和C+是文本，1、如图输入题主的公式。会弹出所键入的公式含有错误，会显示C这里有错，不止C这里有错，前面的A和B也有错，3、因为公式中字母部分外面应...
C语言

70,024

社区成员

243,263

社区内容

发帖
与我相关
我的任务
社区描述
C语言相关问题讨论
社区管理员
  • C语言
  • 花神庙码农
  • 架构师李肯
加入社区
  • 近7日
  • 近30日
  • 至今

加载中

查看更多榜单
社区公告
暂无公告

试试用AI创作助手写篇文章吧

+ 用AI写文章