网页资讯视频图片知道文库贴吧地图采购
进入贴吧全吧搜索

 
 
 
日一二三四五六
       
       
       
       
       
       

签到排名:今日本吧第个签到,

本吧因你更精彩,明天继续来努力!

本吧签到人数:0

一键签到
成为超级会员,使用一键签到
一键签到
本月漏签0次!
0
成为超级会员,赠送8张补签卡
如何使用?
点击日历上漏签日期,即可进行补签。
连续签到:天  累计签到:天
0
超级会员单次开通12个月以上,赠送连续签到卡3张
使用连续签到卡
07月28日漏签0天
c++吧 关注:626,284贴子:2,112,813
  • 看贴

  • 图片

  • 吧主推荐

  • 游戏

  • 0回复贴,共1页
<<返回c++吧
>0< 加载中...

这些怎么改呢?

  • 取消只看楼主
  • 收藏

  • 回复
  • 卷卷是小公主
  • ,
    1
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
# include<stdio.h>
# include<math.h>
# define N 6
/*
*使用雅可比迭代法和高斯-赛德尔迭代法 求解线性方程组
*/
void main(){
float NF2(float *x,float *y);
float A[N][N],b[N],sum=0;
float x[N],y[N]={0},x0[N];
int i,j,n=0;
//输入系数矩阵
for(i=0;i<N;i++){
for(j=0;j<N;j++){
scanf("%f",&A[i][j]);
}
}
//输入常数矩阵
for(i=0;i<N;i++){
scanf("%f",&b[i]);
}
//输入解的初值
for(i=0;i<N;i++){
scanf("%f",&x0[i]);
}
//输出系数矩阵
printf("输出该方程组的系数矩阵:\n");
for(i=0;i<N;i++){
for(j=0;j<N;j++){
printf("%3.1f ",A[i][j]);
}
printf("\n");
}
//输出成数矩阵
printf("输出该方程组的常数矩阵:\n");
for(i=0;i<N;i++){
printf("%3.1f\n",b[i]);
}
//输出解的迭代初值
printf("解该方程组的的迭代初值是:\n");
for(i=0;i<N;i++)
{
printf("%3.1f\n",x0[i]);
}
/*
*利用雅可比迭代法求解线性方程组
*/
for(i=0;i<N;i++)
{
x[i]=x0[i];
}
for(n=0;;n++){
//计算下一个值
for(i=0;i<N;i++){
sum=0;
for(j=0;j<N;j++){
if(j!=i){
sum=sum+A[i][j]*x[j];
}
}
y[i]=(1/A[i][i])*(b[i]-sum);
//sum=0;
}
//判断误差大小
if(NF2(x,y)>0.01){
for(i=0;i<N;i++){
x[i]=y[i];
}
}
else
break;
}
printf("经过%d次雅可比迭代解出方程组的解:\n",n+1);
for(i=0;i<N;i++){
printf("%f ",y[i]);
}
/*
*利用高斯-赛德尔迭代法求解线性方程组
*/
for(i=0;i<N;i++)
{
x[i]=x0[i];
y[i]=0;
}
for(n=0;;n++){
//计算下一个值
for(i=0;i<N;i++){
sum=0;
for(j=0;j<i;j++){
sum=sum+A[i][j]*y[j];
}
for(j=i+1;j<N;j++){
sum=sum+A[i][j]*x[j];
}
y[i]=(1/A[i][i])*(b[i]-sum);
//sum=0;
}
//判断误差大小
if(NF2(x,y)>0.01){
for(i=0;i<N;i++){
x[i]=y[i];
}
}
else
break;
}
printf("\n经过%d次高斯-赛德尔迭代解出方程组的解:\n",n+1);
for(i=0;i<N;i++){
printf("%f ",y[i]);
}
}
//求两个向量差的二范数函数
float NF2(float *x,float *y){
int i;
float z,sum1=0;
for(i=0;i<N;i++){
sum1=sum1+pow(y[i]-x[i],2);
}
z=sqrt(sum1);
return z;


登录百度账号

扫二维码下载贴吧客户端

下载贴吧APP
看高清直播、视频!
  • 贴吧页面意见反馈
  • 违规贴吧举报反馈通道
  • 贴吧违规信息处理公示
  • 0回复贴,共1页
<<返回c++吧
分享到:
©2025 Baidu贴吧协议|隐私政策|吧主制度|意见反馈|网络谣言警示