东南大学数值分析上机题（c语言版）

今年学习数值分析时候，网上找了一大堆代码都是使用matlab的，采用c语言编写的代码少之又少，加之大部分还有缺陷，
本人就根据网上收集加之利用chatgpt写的代码分享出来。*

第一章

一、问题

二、程序

for(i=2;i<=N;i++){
		Sn1=Sn1+1/(i*i-1);
}
for(j=N;j>=2;j--){
		Sn2=Sn2+1/(j*j-1);
}

第二章

一、问题

二、程序

#include<stdio.h>
#include<math.h>
#define eps 0.000001 // 给定容许误差
float f(float x) // 定义函数 f(x)
{
 float f;
 f=x*x*x/3-x; //f(x) 的表达式 ;
 return(f);
}
float df(float x) //df(x)
{
 float df;
 df=x*x-1; //f(x) 导函数的表达式 ;
 return df;
}
void main(void)
{
 float x0,x1,a;
 int k=0;
 printf("请输入初值x0=");
 scanf("%f",&x0);
 do
 {
   a=-f(x0)/df(x0);
   x1=x0+a;
   k++;
   x0=x1;
 }while(fabs(a)>eps);
 printf("迭代次数k=%d,根值是x0=%f",k,x0);
}

#include<stdio.h>
#include<math.h>
#define eps 0.01 // 给定容许误差
float f(float x) // 定义函数 f(x)
{
 float f;
 f=x*x*x/3-x; //f(x) 的表达式 ;
 return(f);
}
float df(float x) //df(x)
{
 float df;
 df=x*x-1; //f(x) 导函数的表达式 ;
 return df;
}
float ddf(float x) //ddf(x)
{
 float ddf;
 ddf=2*x; //f(x) 二阶导函数的表达式 ;
 return ddf;
}
void main(void)
{
 float x0,x1,a,dt=0;
 int k=0;
 do{
 	dt=dt+eps;
 	k++;
 }while((f(-dt)*f(dt)<0)&&(df(dt)!=0)&&(2*dt>=-f(-dt))&&(2*dt>=f(dt)/df(dt)));
 printf("请输入初值x0=");
 scanf("%f",&x0);
 if(x0<-1)
 do
 {
   a=-f(x0)/df(x0);
   x1=x0+a;
   k++;
   x0=x1;
 }while(fabs(a)>eps);
 else if(x0>-1&&x0<-dt)
 do
 {
   a=-f(x0)/df(x0);
   x1=x0+a;
   k++;
   x0=x1;
 }while(fabs(a)>eps);
 else if(x0>-dt&&x0<dt)
 do
 {
   a=-f(x0)/df(x0);
   x1=x0+a;
   k++;
   x0=x1;
 }while(fabs(a)>eps);
 else if(x0>dt&&x0<1)
 do
 {
   a=-f(x0)/df(x0);
   x1=x0+a;
   k++;
   x0=x1;
 }while(fabs(a)>eps);
 else if(x0>1)
 do
 {
   a=-f(x0)/df(x0);
   x1=x0+a;
   k++;
   x0=x1;
 }while(fabs(a)>eps);
 printf("收敛域=%f,迭代次数是=%d\n",dt,k);
 printf("根＝%f\n",x0);
}

三、结果

第三章（1）

一、问题一

二、程序

#include<stdio.h>
#include<math.h>
void main(void)
{
	int n,i,j,k,t;
	float a[10][11],s[10],x[10],s2,max,sum,b;
	printf("请输入n= ");
	scanf("%d",&n);
	printf("请输入矩阵A= \n");
	for(i=0;i<=n;i++)
	{
//		printf("第",i,"行");
		for(j=0;j<=(n+1);j++)
			scanf("%f",&a[i][j]);	
				
		printf("第%d行 \n",(i+1));
		for(j=0;j<=(n+1);j++)
			printf("%f",a[i][j]);

	}
	//变换高斯消去矩阵
	for(k=0;k<n;k++){
		j=0;
		for(i=k;i<=n;i++)
		{
			s[j]=a[i][k];
			j++;
		}
		max=fabs(s[0]);
		t=k;
		for(i=1;i<j;i++)
		{
			if(fabs(s[i])>max)
			{
				max=fabs(s[i]);
				t=k+i;
			}
		}
		
		for(j=0;j<=n+1;j++)
		{
			b=a[k][j];
			a[k][j]=a[t][j];
			a[t][j]=b;
		}
		
		for(i=(k+1);i<=n;i++)
		{
			s2=a[i][k]/a[k][k];
			for(j=k;j<=(n+1);j++)
				a[i][j]=a[i][j]-s2*a[k][j];
		}
		printf("\n");
	}
	printf("\n");
	for(i=n;i>=0;i--)
	{
		sum=0;
		x[n+1]=a[n][n+1]/a[n][n];
		for(j=i;j<=n;j++)
		{
			sum+=a[i-1][j]*x[j+1];	
		}
		x[i]=(a[i-1][n+1]-sum)/a[i-1][i-1];	
	}
	printf("方程的解为 \n");	
	for(i=0;i<=n;i++)
		printf("%f, ",x[i+1]);
	
}

三、结果

第三章（2）

一、问题

36、逐次超松弛迭代法
（1）编制解n 阶线性方程组Ax=b 的SOR方法的通用程序；
（2）对于35题中所给的线性方程组，取松弛因子w= i/50(i = 1,2…容许误差0.000005”，打印松弛因子、迭代次数、最佳松弛因子及解向量。

二、程序

#include<stdio.h>
#include<math.h>
#define eps 0.5e-5 
void main(void)
{
	int i,j,l,k=999;
	float w,t,p=0;
	float m[9];
	float sum;
	float a[9][9]={{31,-13,0,0,0,-10,0,0,0},{-13,35,-9,0,-11,0,0,0,0},{0,-9,31,-10,0,0,0,0,0},{0,0,-10,79,-30,0,0,0,-9},{0,0,0,-30,57,-7,0,-5,0},
				      {0,0,0,0,-7,47,-30,0,0},{0,0,0,0,0,-30,41,0,0},{0,0,0,0,-5,0,0,27,-2},{0,0,0,-9,0,0,0,-2,29}};
	float b[9]={-15,27,-23,0,-20,12,-7,7,10};
	float max(float m[9]);
	
	for(t=1;t<=99;t++)
	{
		l=0;
		float x0[9]={1,1,1,1,1,1,1,1,1};
		float x1[9]={1,1,1,1,1,1,1,1,1};
		w=t/50;
		do{
			for(i=0;i<9;i++)
			{
				x0[i]=x1[i];
			}
			for(i=0;i<9;i++)
			{
				sum=0;
				for(j=0;j<i;j++)
				{
					sum=sum+a[i][j]*x1[j];
				}
				for(j=i+1;j<9;j++)
				{
					sum=sum+a[i][j]*x0[j];	
				}
				x1[i]=(1-w)*x0[i]+w*(b[i]-sum)/a[i][i];
			}
			for(i=0;i<9;i++){
				m[i]=x1[i]-x0[i];
			}
			l++;
		}while(max(m)>=eps);
		if(max(m)<=eps)
		{
   			printf("\n 松弛因子：%f", w);
			printf("迭代次数：%d", l);
//			printf("最佳松弛因子：%f", w);
			printf("解向量：");
			for (i = 0; i < 9; i++) {
			    printf("%f ", x1[i]);
			}
		}
		//寻找最佳松弛因子
		if(l<k)
		{
			k=l;
			p=w;
		}
		
	}
	printf("最佳松弛因子：%f", p);
	printf("最佳松弛因子的迭代次数为：%d次", k);
}
float max(float m[9]){
	float k;
	k=(fabs(m[0]));
	for(int i=0;i<9;i++)
	{
		if(fabs(m[i])>k)
			k=fabs(m[i]);
	}
	return k;
}

三、结果