import numpy as np
成都創(chuàng)新互聯(lián)公司是一家以網(wǎng)站建設(shè)公司、網(wǎng)頁(yè)設(shè)計(jì)、品牌設(shè)計(jì)、軟件運(yùn)維、營(yíng)銷推廣、小程序App開發(fā)等移動(dòng)開發(fā)為一體互聯(lián)網(wǎng)公司。已累計(jì)為會(huì)所設(shè)計(jì)等眾行業(yè)中小客戶提供優(yōu)質(zhì)的互聯(lián)網(wǎng)建站和軟件開發(fā)服務(wù)。
def solve_quad(a,b,c):
if a == 0:
print('您輸入的不是二次方程!')
else:
delta = b*b-4*a*c
x = -b/(2*a)
if delta == 0:
print('方程有惟一解,X=%f'%(x))
return x
elif delta 0:
x1 = x-np.sqrt(delta)/(2*a)
x2 = x+np.sqrt(delta)/(2*a)
print('方程有兩個(gè)實(shí)根:X1=%f,X2=%f'%(x1,x2))
return x1,x2
else:
x1 = (-b+complex(0,1)*np.sqrt((-1)*delta))/(2*a)
x2 = (-b-complex(0,1)*np.sqrt((-1)*delta))/(2*a)
print(x1,x2)
return x1,x2
Python
是完全面向?qū)ο蟮恼Z(yǔ)言。函數(shù)、模塊、數(shù)字、字符串都是對(duì)象。并且完全支持繼承、重載、派生、多繼承,有益于增強(qiáng)源代碼的復(fù)用性。Python支持重載運(yùn)算符和動(dòng)態(tài)類型。相對(duì)于Lisp這種傳統(tǒng)的函數(shù)式編程語(yǔ)言,Python對(duì)函數(shù)式設(shè)計(jì)只提供了有限的支持。有兩個(gè)標(biāo)準(zhǔn)庫(kù)(functools, itertools)提供了Haskell和Standard ML中久經(jīng)考驗(yàn)的函數(shù)式程序設(shè)計(jì)工具。
#includeiostream.h
#includemath.h
#includeconio.h
const int N=200;
//帶入原函數(shù)后所得的值
double f(float x)
{
return (x*x*x-1.8*x*x+0.15*x+0.65);
}
//帶入一階導(dǎo)函數(shù)后所得的值
double f1(double x)
{
return (3*x*x-3.6*x+0.15);
}
//牛頓迭代函數(shù)
double F(double x)
{
double x1;
x1=x-1.0*f(x)/f1(x);
return (x1);
}
void main()
{
double x0,D_value,x1,y[4];
int k=0,count=0;
for(;;)
{
if(count==3)break;
cout"輸入初始值:";
cinx0;
do
{
k++;
x1=F(x0);
D_value=fabs(x1-x0);
x0=x1;
}
while((D_value0.000005)(k=N));
for(int j=0,flag=0;jcount;j++)
{
if(fabs(y[j]-x1)0.000005)
{ flag=1;
cout"該數(shù)值附近的根已經(jīng)求出,請(qǐng)重新?lián)Q近似值"endl;
break;
}
}
if(flag==1)
continue;
else
{
cout"方程的一個(gè)根:"x1","" 迭代次數(shù)為:"kendl;
y[count]=x1;
count++;
}
//else
//cout"計(jì)算失敗!"endl;
}
}
//你的程序其實(shí)沒問(wèn)題,牛頓迭代法本身循環(huán)一次只能找到一個(gè)答案,只要再建一個(gè)循環(huán)控制使
//用迭代法的次數(shù)和判斷根的個(gè)數(shù)就行。我又加了一個(gè)判斷是否有重復(fù)的根的循環(huán)。
//希望能對(duì)你有所幫助。
1:二分法
求根號(hào)5
a:折半:?????? 5/2=2.5
b:平方校驗(yàn):? 2.5*2.5=6.255,并且得到當(dāng)前上限2.5
c:再次向下折半:2.5/2=1.25
d:平方校驗(yàn):1.25*1.25=1.56255,得到當(dāng)前下限1.25
e:再次折半:2.5-(2.5-1.25)/2=1.875
f:平方校驗(yàn):1.875*1.875=3.5156255,得到當(dāng)前下限1.875
每次得到當(dāng)前值和5進(jìn)行比較,并且記下下下限和上限,依次迭代,逐漸逼近平方根:
代碼如下:
import math
from math import sqrt
def sqrt_binary(num):
x=sqrt(num)
y=num/2.0
low=0.0
up=num*1.0
count=1
while abs(y-x)0.00000001:
print count,y
count+=1
if (y*ynum):
up=y
y=low+(y-low)/2
else:
low=y
y=up-(up-y)/2
return y
print(sqrt_binary(5))
print(sqrt(5))
2:牛頓迭代
仔細(xì)思考一下就能發(fā)現(xiàn),我們需要解決的問(wèn)題可以簡(jiǎn)單化理解。
從函數(shù)意義上理解:我們是要求函數(shù)f(x) = x2,使f(x) = num的近似解,即x2 - num = 0的近似解。
從幾何意義上理解:我們是要求拋物線g(x) = x2 - num與x軸交點(diǎn)(g(x) = 0)最接近的點(diǎn)。
我們假設(shè)g(x0)=0,即x0是正解,那么我們要做的就是讓近似解x不斷逼近x0,這是函數(shù)導(dǎo)數(shù)的定義:
從幾何圖形上看,因?yàn)閷?dǎo)數(shù)是切線,通過(guò)不斷迭代,導(dǎo)數(shù)與x軸的交點(diǎn)會(huì)不斷逼近x0。