数值分析作业

发布 2022-08-27 20:10:28 阅读 1792

#include<>

#include<>

float max(float m, float n);

void main()

float f,x[50],y[50];

int i;

x[0]=1;

y[0]=0;

f=1;for(i=1;f>0.00005;i++)

float max(float m, float n)//求2者最大值函数。

float z;

if (m>n)

z=m;else

z=n;return (z);

输出结果:迭代27次后:x=0.900540 y=-0.869626

无穷范数为0.000040

心得:此题并不难,但修改了很多次才得到结果,主要一下注意事项:

1 .c语言中sin函数是弧度标示不是角度。

2 .数组的维数不能太少,开始设置迭代20次,结果没计算出来。

3.开始的时候x2开根号没有取正负,后来直接除以个x避免开根号。

第二题:#include<>

#include<>

float x[4],d[4];/定义为全局变量。

void jfc(float a[4][4],float b[4]);求解方程组函数声明。

void cj(float t[4][4],float r[4]);

void main()

float f[4]=,h[4],z[4],q[4];

float a[4][4]=,

int p,m;

jfc(a,f

for (p=0;p<4;p++)

z[p]=x[p];

cj(a,x);

for (p=0;p<4;p++)

h[p]=f[p]-d[p];

jfc(a,h);

for (p=0;p<4;p++)

z[p]=x[p]+z[p];

for (p=0;p<4;p++)

printf(" z[%d]=%f ",p,z[p]);

void jfc(float a[4][4],float b[4]) 求解方程组函数。

int n=4,i,j,k,mi;

float tmp,mx;

for(i=0;i

for(j=i+1;j

x[n-1]=b[n-1]/a[n-1][n-1];

for(i=n-2;i>=0;i--)

void cj(float t[4][4],float r[4]) 矩阵与向量乘。

int i,j,k;

for(i=0;i<4;i++)

for(j=0;j<4;j++)

d[i]+=t[i][j]*r[j];

return;

结果:迭代2次后,运行结果为:x1= 1.5442+0.2253 i 和 x2=0.1467+0.3898 i

心得:此题设x1 和x2的实部和虚部,然后代换进去成为一个4元方程,然后按照迭代改善方法处理。

第三题:#include<>

#include<>

float x[4];/定义为全局变量。

void jfc(float a[4][4],float b[4]);求解方程组。

void main()

float f[4]=;使用作第一次迭代。

float d[100],e[4];

float mx;

float a[4][4]=,初始矩阵a

int k,p,h;

for(k=0;k<4;k++)对初始矩阵处理求的a-7e

a[k][k]-=8.1;

d[0]=0;

d[1]=1;

for(k=2;fabs(d[k-1]-d[k-2])>0.01;k++)多次迭代,直到fabs(d[k-1]-d[k-2])《0.01条件满足。

printf("%d次迭代,特征值为%f",k,7+mx);

for (p=0;p<4;p++)

printf(" f[%d]=%f ",p,f[p]);

void jfc(float a[4][4],float b[4]) 求解方程组函数。

int n=4,i,j,k,mi;

float tmp,mx;

for(i=0;i

for(j=i+1;j

x[n-1]=b[n-1]/a[n-1][n-1];

for(i=n-2;i>=0;i--)

迭代6次后,运行结果为:

特征值7.46159

特征向量为。

心得:此题迭代方式并不复杂,由于迭代的初始矩阵始终不变因此设置为全局变量,可以反复使用,每次只需要代入新的vk即可。

数值分析作业

12.求在 0,1 上的一次最佳平方逼近多项式与二次最佳平方逼近多项式。函数 function s zjpfbj n,a,b 创建一个函数,里面填入次数,和区间范围。base inline x j 1 x j 定义多项式。quan inline 1 x 权函数。a zeros n 1 y zeros...

数值分析作业

实验2.2算法设计与比较。实验目的 编制不同方法的matlab程序,这些方法的计算效果和特点。问题提出 非线性方程的数值解法很多,不同的方法效果如何,要靠计算的实践来分析 比较。实验内容 考虑下列算法 1 牛顿法 2 弦割法 3 抛物线法。分别编写它们的matlab程序。牛顿法程序 function...

数值分析作业

数值分析上机作业。zhh mushan 班级 土木工程一班。姓名 赵惠惠。班级 2010510359 二零一一年六月 舍入误差与有效数。一 问题说明。1 设,其精确值为。1 编制按从小到大的顺序计算通用程序。2 编制按从小到大的顺序计算通用程序。3 按两种顺序分别计算并指出有效位数。编制程序用单精度...