这是弄的一个牛顿迭代的程序,之前都没出现问题,这一下就出现了无法转换的问题,有点懵,麻烦各位给点意见,谢谢。
1条回答 默认 最新
- xuelanghanbao 2022-11-09 20:01关注
x是符号,不是值,所以根据x计算出的T、c也是一个表达式而不是一个具体的值,所以没法和一个具体的值1e-6来比较,要在使用前先把x转换为求解的值xk并计算才能到的一个具体的值。另外,注意到程序中在while里并没有更新c的值,我也添加了一下。
clear; clc; syms x; y = (exp((-x^2)/2))*((2*pi)^(-1/2)); xk=0.5; i=0; T=int(y,x,0,xk); Tk = eval(subs(T,x,xk)); c=abs(Tk-0.45); while c>=1e-6 xk=xk-(T-0.45)/((exp((-xk^2)/2))*((2*pi)^(-1/2))); i=i+1; T=int(y,x,0,xk); Tk = eval(subs(T,x,xk)); c=abs(Tk-0.45); end
本回答被题主选为最佳回答 , 对您是否有帮助呢?解决 5无用 2