十年网站开发经验 + 多家企业客户 + 靠谱的建站团队
量身定制 + 运营维护+专业推广+无忧售后,网站问题一站解决
1、用C 语言求雅克比迭代法、高斯赛德尔迭代法(急求 拜谢) 20 1)、问题提出试分别选用Jacobi迭代法,Gauss-Seidol迭代法计算其解。
创新互联公司专注于江门企业网站建设,成都响应式网站建设公司,商城网站建设。江门网站建设公司,为江门等地区提供建站服务。全流程按需求定制网站,专业设计,全程项目跟踪,创新互联公司专业和态度为您提供的服务
2、if(fabs(y[i]-x[i])0.00001) j++;if(j==n)break;} 这儿,除非fabs(y[i]-x[i])0.00001 每次都成立 ,不然 j 是不可能 等于n 的,这样,你的for(k=0;k++)就没有break的机会 ,就会一直循环下去。
3、(2)考察迭代法的收敛性,若均收敛,对比两种方法的收敛速度。
4、但必须指出,上述所得的向量序列{x(k)}是否收敛于Ax=b的解是有条件的,而且即使同样是收敛的,还有收敛速度快慢的问题。
1、[33nodes.rar] - 配电网潮流计算源程序,可直接运行的,修改输入数据即可 [hhu_pfcp.rar] - 基于MATPOWER修改而成的潮流计算软件,可以采用直流潮流法、牛顿拉夫逊法、PQ分解法与高斯-塞代尔法解算IEEE9-300节点的潮流。
2、[33nodes.rar] - 配电网潮流计算源程序,可直接运行的,修改输入数据即可 [hhu_pfcp.rar] - 基于MATPOWER修改而成的潮流计算,可以采用直流潮流法、牛顿拉夫逊法、PQ分解法与高斯-塞代尔法解算IEEE9-300节点的潮流。
对于求平方根,变成方程模式为f(x)=x^2-a,即求此方程的实根;下面编写了两个function函数,可以直接调用。
if p = 1 迭代法不收敛return end while r e x0 = x;x = B*x0 + f;k = k + 1;r = norm (x-x0,inf);end 所求解为x 迭代次数为k 自己以前编的。。
楼主太坑了,一个死循环,害的我好多原来打开的东西都没有保存,只能强制关掉。
你的x没有初始化,导致当你迭代不收敛时x没有赋值,肯定x作为返回值报错。
b取zeros(20,1)的时候,由x=zeros(size(b)),即x初值为0和 x=D\((D-A)*x+b),即x每次迭代后仍为0可知,x始终为0,没有问题。