#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 按两种顺序分别计算并指出有效位数。编制程序用单精度...