tensorflow随笔-不动点迭代求一元方程

#!/usr/bin/env python2
# -*- coding: utf-8 -*-
"""
Created on Thu Sep  6 10:16:37 2018不动点法求解f(x)=x
"""import tensorflow as tf
import numpy as npdef f(x):y=tf.cos(x)+tf.sin(x)return ydef body(x,fx,tol,i,n):x=f(x)fx=f(x)return (tf.Print(x,[x],"x:"),tf.Print(fx,[fx],"fx:"),tf.Print(tol,[tol],"tol:"),tf.Print(i+1,[i],"i:"),tf.Print(n,[n],"n:"))def c(x,fx,tol,i,n):t1=tf.greater(tf.abs(tf.subtract(fx,x)),tol)t2=tf.less(i,n)return tf.logical_and(t1,t2)x = tf.placeholder(tf.float32,shape=(),name="myx")
tol= tf.placeholder(tf.float32,shape=(),name="mytol")
fx = tf.constant(0,dtype=tf.float32,name="myfx")
i = tf.constant(0,dtype=tf.int32,name="myi")
n = tf.constant(0,dtype=tf.int32,name="myn")input_dict={x:0.,fx:np.cos(0)+np.sin(0),tol:1e-8,i:0,n:100}
res = tf.while_loop(c, body, loop_vars=[x,fx,tol,i,n])with tf.Session() as sess:    y=sess.run(res,feed_dict=input_dict)print y

tol:[1e-08]
x:[1.25872827]
fx:[1.25872803]
n:[100]
i:[98]tol:[1e-08]x:[1.25872803]
fx:[1.25872827]

n:[100]
x:[1.25872827]i:[99]tol:[1e-08]
fx:[1.25872803]

(1.2587283, 1.258728, 1e-08, 100, 100)

tensorflow随笔-不动点迭代求一元方程相关推荐

不动点迭代（Fixed Point Iteration）
题目:不动点迭代(Fixed Point Iteration) 本篇介绍不动点迭代(Fixed Point Iteration).之所以学习不动点迭代是由于近来看到了FPC算法,即Fixed Poin ...
FPI(Fixed-point Iteration)不动点迭代法——迭代求方程的方法
一不动点是什么? 不动点,其实定义比较简单,对于一些方程,例如f(x)=cosx,那么令cosx=x的点就是函数的不动点,说白了,就是y=x这条直线与函数曲线的交点.这个不动点有什么用呢?请继续往下 ...
不动点迭代法（Fixed Point Iteration）迭代求根的python程序
迭代法的作用许多复杂的求解问题,都可以转换成方程f(x)=0的求解问题.这一系列的解叫做方程的根.对于非线性方程的求解,在自变量范围内往往有多个解,我们将此变化区域分为多个小的子区间,对每个区间进行 ...
数值分析——求方程解的不动点迭代法和斯特芬森法(Python实现)
一.不动点迭代法求方程的解 import sympy #迭代的方程为 f(x)=pow(x,3)-x-1 def psi(x):return pow(x+1,1/3)def dif(x0): #求导函 ...
用matlab求不动点迭代,matlab实现不动点迭代、牛顿法、割线法
不动点迭代 function xc = fpi( g, x0, tol ) x(1) = x0; i = 1; while 1 x(i + 1) = g(x(i)); if(abs(x(i+1) - ...
其他算法-不动点迭代
不动点迭代不动点迭代法(Fixed Point Iteration)又叫简单迭代法,对于一个非线性方程 f(x)=0f(x)=0f(x)=0 将其转换成以下形式: x=φ(x)x=\varphi ( ...
非线性方程(组)：一维非线性方程（一）二分法、不动点迭代、牛顿法 [MATLAB]...
1. 二分法(Bisection) 1) 原理 [介值定理] 对于连续的一元非线性函数,若其在两个点的取值异号,则在两点间必定存在零点. [迭代流程] 若左右两端取值不同,则取其中点,求其函数值,取中 ...
不动点迭代法求函数根（非线性方程求解）
问题描述在[a,b]区间内寻找方程x**5-2*x-1=0的根的初始近似值位置,确定不动点迭代的初始点(可能有多个),然后使用不动点迭代法求方程的根(可能有多个根).前后两次迭代的差的绝对值小于de ...
牛顿迭代法（Newton’s Method）迭代求根的Python程序
迭代法的作用许多复杂的求解问题,都可以转换成方程f(x)=0的求解问题.这一系列的解叫做方程的根.对于非线性方程的求解,在自变量范围内往往有多个解,我们将此变化区域分为多个小的子区间,对每个区间进行 ...

tensorflow随笔-不动点迭代求一元方程

tensorflow随笔-不动点迭代求一元方程相关推荐

最新文章

热门文章