今年学习数值分析时候,网上找了一大堆代码都是使用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;
}