色哟哟视频在线观看-色哟哟视频在线-色哟哟欧美15最新在线-色哟哟免费在线观看-国产l精品国产亚洲区在线观看-国产l精品国产亚洲区久久

0
  • 聊天消息
  • 系統消息
  • 評論與回復
登錄后你可以
  • 下載海量資料
  • 學習在線課程
  • 觀看技術視頻
  • 寫文章/發帖/加入社區
會員中心
創作中心

完善資料讓更多小伙伴認識你,還能領取20積分哦,立即完善>

3天內不再提示

帶大家領略sqrt的神奇之處:開平方的7種算法介紹

算法與數據結構 ? 來源:未知 ? 2019-11-28 09:29 ? 次閱讀

作者:nash_

sqrt()函數,是絕大部分語言支持的常用函數,它實現的是開方運算;開方運算最早是在我國魏晉時數學家劉徽所著的《九章算術》被提及。今天寫了幾個函數加上國外大神的幾個神級程序帶大家領略sqrt的神奇之處。

1、古人算法(暴力法)

原理:從0開始0.00001,000002...一個一個試,直到找到x的平方根,代碼如下:

publicclassAPIsqrt{

staticdoublebaoliSqrt(doublex){

finaldouble_JINGDU=1e-6;
doublei;
for(i=0;Math.abs(x-i*i)>_JINGDU;i+=_JINGDU)
;
returni;
}

publicstaticvoidmain(String[]args){
doublex=3;
doubleroot=baoliSqrt(x);
System.out.println(root);
}

測試結果:

1、7320509999476947

2、牛頓迭代法

計算機科班出身的童鞋可能首先會想到的是《數值分析》中的牛頓迭代法求平方根。原理是:隨意選一個數比如說8,要求根號3,我們可以這么算:

(8 + 3/8) = 4.1875

(4.1875 + 3/4.1875) = 2.4519

(2.4519 + 3/2.4519) = 1.837

(1.837 + 3/1.837) = 1.735

做了4步基本算出了近似值了,這種迭代的方式就是傳說中的牛頓迭代法了,代碼如下:

publicclassAPIsqrt{

staticdoublenewtonSqrt(doublex){

if(x0){
System.out.println("負數沒事開什么方");
return-1;
}
if(x==0)
return0;

double_avg=x;
doublelast_avg=Double.MAX_VALUE;
finaldouble_JINGDU=1e-6;

while(Math.abs(_avg-last_avg)>_JINGDU){
last_avg=_avg;
_avg=(_avg+x/_avg)/2;
}
return_avg;
}

publicstaticvoidmain(String[]args){
doublex=3;
doubleroot=newtonSqrt(x);
System.out.println(root);
}
}

測試結果:

17320508075688772

3、暴力-牛頓綜合法

原理:還是以根號3為例,先用暴力法講根號3逼近到1.7,然后再利用上述的牛頓迭代法。雖然沒有用牛頓迭代好,但是也為我們提供一種思路。代碼如下:

publicclassAPIsqrt{

staticdoublebaoliAndNewTonSqrt(doublex){

if(x0){
System.out.println("負數沒事開什么方");
return-1;
}
if(x==0)
return0;

doublei=0;
double_avg;
doublelast_avg=Double.MAX_VALUE;
for(i=0;i*i0.1);
_avg=i;

finaldouble_JINGDU=1e-6;

while(Math.abs(_avg-last_avg)>_JINGDU){
last_avg=_avg;
_avg=(_avg+x/_avg)/2;
}
return_avg;
}

publicstaticvoidmain(String[]args){
doublex=3;
doubleroot=baoliAndNewTonSqrt(x);
System.out.println(root);
}
}

測試結果:

1、7320508075689423

4、二分開方法

原理:還是以3舉例:

(0+3)/2 = 1.5, 1.5^2 = 2.25, 2.25 < 3;

(1.5+3)/2 = 2.25, 2.25^2 = 5.0625, 5.0625 > 3;

(1.5+2.25)/2 = 1.875, 1.875^2 = 3.515625; 3.515625>3;

直到前后兩次平均值只差小于自定義精度為止,代碼如下:

publicclassAPIsqrt{

staticdoubleerfenSqrt(doublex){

if(x0){
System.out.println("負數沒事開什么方");
return-1;
}
if(x==0)
return0;

finaldouble_JINGDU=1e-6;
double_low=0;
double_high=x;
double_mid=Double.MAX_VALUE;
doublelast_mid=Double.MIN_VALUE;

while(Math.abs(_mid-last_mid)>_JINGDU){

last_mid=_mid;
_mid=(_low+_high)/2;
if(_mid*_mid>x)
_high=_mid;
if(_mid*_midreturn_mid;

}

publicstaticvoidmain(String[]args){
doublex=3;
doubleroot=erfenSqrt(x);
System.out.println(root);
}
}

測試結果:

1、732051134109497

5、計算 (int)(sqrt(x))算法

PS:此算法非博主所寫

原理:空間換時間,細節請大家自行探究,代碼如下:

publicclassAPIsqrt2{
finalstaticint[]table={0,16,22,27,32,35,39,42,45,48,50,53,
55,57,59,61,64,65,67,69,71,73,75,76,78,80,81,83,84,
86,87,89,90,91,93,94,96,97,98,99,101,102,103,104,
106,107,108,109,110,112,113,114,115,116,117,118,119,
120,121,122,123,124,125,126,128,128,129,130,131,132,
133,134,135,136,137,138,139,140,141,142,143,144,144,
145,146,147,148,149,150,150,151,152,153,154,155,155,
156,157,158,159,160,160,161,162,163,163,164,165,166,
167,167,168,169,170,170,171,172,173,173,174,175,176,
176,177,178,178,179,180,181,181,182,183,183,184,185,
185,186,187,187,188,189,189,190,191,192,192,193,193,
194,195,195,196,197,197,198,199,199,200,201,201,202,
203,203,204,204,205,206,206,207,208,208,209,209,210,
211,211,212,212,213,214,214,215,215,216,217,217,218,
218,219,219,220,221,221,222,222,223,224,224,225,225,
226,226,227,227,228,229,229,230,230,231,231,232,232,
233,234,234,235,235,236,236,237,237,238,238,239,240,
240,241,241,242,242,243,243,244,244,245,245,246,246,
247,247,248,248,249,249,250,250,251,251,252,252,253,
253,254,254,255};

/**
*Afasterreplacementfor(int)(java.lang.Math.sqrt(x)).Completely
*accurateforx2147483648(i.e.2^31)...
*/
staticintsqrt(intx){
intxn;

if(x>=0x10000){
if(x>=0x1000000){
if(x>=0x10000000){
if(x>=0x40000000){
xn=table[x>>24]<8;
}else{
xn=table[x>>22]<7;
}
}else{
if(x>=0x4000000){
xn=table[x>>20]<6;
}else{
xn=table[x>>18]<5;
}
}

xn=(xn+1+(x/xn))>>1;
xn=(xn+1+(x/xn))>>1;
return((xn*xn)>x)?--xn:xn;
}else{
if(x>=0x100000){
if(x>=0x400000){
xn=table[x>>16]<4;
}else{
xn=table[x>>14]<3;
}
}else{
if(x>=0x40000){
xn=table[x>>12]<2;
}else{
xn=table[x>>10]<1;
}
}

xn=(xn+1+(x/xn))>>1;

return((xn*xn)>x)?--xn:xn;
}
}else{
if(x>=0x100){
if(x>=0x1000){
if(x>=0x4000){
xn=(table[x>>8])+1;
}else{
xn=(table[x>>6]>>1)+1;
}
}else{
if(x>=0x400){
xn=(table[x>>4]>>2)+1;
}else{
xn=(table[x>>2]>>3)+1;
}
}

return((xn*xn)>x)?--xn:xn;
}else{
if(x>=0){
returntable[x]>>4;
}
}
}

return-1;
}
publicstaticvoidmain(String[]args){
System.out.println(sqrt(65));

}
}

測試結果:8

6、最快的sqrt算法

PS:此算法非博主所寫

這個算法很有名,大家可能也見過,作者是開發游戲的,圖形算法中經常用到sqrt,作者才寫了一個神級算法,和他那神秘的0x5f3759df,代碼如下

#include
floatInvSqrt(floatx)
{
floatxhalf=0.5f*x;
inti=*(int*)&x;//getbitsforfloatingVALUE
i=0x5f375a86-(i>>1);//givesinitialguessy0
x=*(float*)&i;//convertbitsBACKtofloat
x=x*(1.5f-xhalf*x*x);//Newtonstep,repeatingincreasesaccuracy
returnx;
}

intmain()
{
printf("%lf",1/InvSqrt(3));

return0;
}

測試結果:

感興趣的朋友可以參考http://wenku.baidu.com/view/a0174fa20029bd64783e2cc0.html 是作者解釋這個算法的14頁論文《Fast Inverse Square Root》

7、一個與算法6相似的算法

PS:此算法非博主所寫

代碼如下:

#include
floatSquareRootFloat(floatnumber){
longi;
floatx,y;
constfloatf=1.5F;

x=number*0.5F;
y=number;
i=*(long*)&y;
i=0x5f3759df-(i>>1);
y=*(float*)&i;
y=y*(f-(x*y*y));
y=y*(f-(x*y*y));
returnnumber*y;
}

intmain()
{
printf("%f",SquareRootFloat(3));

return0;
}

測試結果:


聲明:本文內容及配圖由入駐作者撰寫或者入駐合作網站授權轉載。文章觀點僅代表作者本人,不代表電子發燒友網立場。文章及其配圖僅供工程師學習之用,如有內容侵權或者其他違規問題,請聯系本站處理。 舉報投訴
  • 算法
    +關注

    關注

    23

    文章

    4607

    瀏覽量

    92840
  • 函數
    +關注

    關注

    3

    文章

    4327

    瀏覽量

    62573

原文標題:開平方的 7 種算法

文章出處:【微信號:TheAlgorithm,微信公眾號:算法與數據結構】歡迎添加關注!文章轉載請注明出處。

收藏 人收藏

    評論

    相關推薦

    單片機開平方的快速算法

    本帖最后由 scan0123 于 2012-8-13 14:52 編輯 單片機開平方的快速算法為工作的需要,要在單片機上實現開根號的操作。目前開平方的方法大部分是用牛頓迭代法。我在查了一些資料
    發表于 08-02 14:40

    VeriLog 開平方怎么做?

    VeriLog 開平方怎么做?
    發表于 10-11 14:09

    關于開平方算法

    電機算法中會用到開平方算法,是個相對復雜和耗時些的數學算法C math庫的double sqrt耗時太長,不能用我沒有用過IQmath庫,自
    發表于 09-27 11:27

    程序中定義了一個變量進行sqrt()開根號處理與預想的結果不對

    我在程序中定義了一個變量對其進行sqrt()開根號處理,結果與預想的不對。在TMS320F28335里面,執行開平方有什么注意的嗎?
    發表于 11-05 11:45

    圖像算法有什么神奇之處

    吳軍在《數學之美》中,用通俗易懂的語言,讓非專業讀者領略了一把數學的魅力。數學作為自然科學的基礎,它的應用遍及物理、工程、生物、化學和經濟等多個領域,在拍照這樣一個似乎不太相關的領域,數學算法再一次
    發表于 09-17 09:05

    STC12C5A60S2無法實現開平方算法怎么回事

    求助!STC12C5A60S2無法實現開平方算法sqrt函數),以及atan2和asin怎么辦?我已經包含了相關的頭文件了,但是編譯通不過。
    發表于 05-20 09:07

    在單片機中開平方會用到哪些算法

    C語言中開平方算法中要開平方的話,可以在頭文件中加#include .然后調sqrt(n);函數即可.但在單片機中要開平方.可以用到下面
    發表于 07-15 07:03

    基于定點DSP的浮點開平方算法的實現

    本文提出了基于TMS320C2XX 定點DSP 的浮點開平方算法,給出了實現方法及程序清單。實踐證明該方法具有精度高、運算速度快、程序簡單等特點。以美國 TI 公司的TMS320C2XX 為代
    發表于 07-31 08:11 ?42次下載

    基于逐次逼近的VHDL開平方算法

    Cyclone 器件為在FPGA 上實現低成本的數字信號處理(DSP)系統提供了一個理想的平臺,成為大規模生產時的最佳方案.。本文介紹了在以EP1C6 芯片為核心的硬件系統中,運用VHDL 實現
    發表于 11-30 14:14 ?32次下載

    MCS-51單片機實用子程序庫實驗(六)

    標號: FSQR 功能:浮點數開平方(快速逼近算法)入口條件:操作數在[R0]中。出口信息:OV=0時,平方根仍在[R0]中,OV=1時,負數開平方
    發表于 01-08 10:11 ?54次下載

    電源變壓器設計與計算鐵心面積公式

    電源變壓器設計與計算鐵心面積公式:鐵心面積=容量開平方;以100W為例;100開平方=10平方厘米.電子管輸出變壓器的鐵心面積=容量開平方*2.以15W輸出變壓器為例:15
    發表于 10-23 14:39 ?7876次閱讀
    電源變壓器設計與計算鐵心面積公式

    電力系統微機保護中開平方運算的一新的快速算法

    電力系統微機保護中開平方運算的一新的快速算法
    發表于 11-02 11:03 ?8次下載

    單片機快速開平方算法資料和代碼免費下載

    C語言中開平方算法中要開平方的話,可以在頭文件中加#include 。然后調sqrt(n);函數即可。但在單片機中要開平方。可以用到下面
    發表于 09-24 17:18 ?1次下載
    單片機快速<b class='flag-5'>開平方</b>的<b class='flag-5'>算法</b>資料和代碼免費下載

    單片機快速開平方算法

    C語言中開平方算法中要開平方的話,可以在頭文件中加#include .然后調sqrt(n);函數即可.但在單片機中要開平方
    發表于 11-11 13:36 ?1次下載
    單片機快速<b class='flag-5'>開平方</b>的<b class='flag-5'>算法</b>

    NI Multisim 10經典教程分享--除法與開平方運算電路

    NI Multisim 10經典教程分享--除法與開平方運算電路
    的頭像 發表于 02-08 09:18 ?2167次閱讀
    主站蜘蛛池模板: 午夜神器18以下不能进免费| 亚洲电影二区| 亚洲人成无码久久久AAA片| 69国产精品成人无码视频| vidosgratis tv少女| 国产精品一区二区亚瑟不卡| 九九热精品免费观看| 青柠在线观看视频在线| 亚洲精品视频免费观看| adc影院在线| 国语自产精品一区在线视频观看| 男女床上黄色| 亚洲黄色成人| 百性阁综合社区| 久久精品国产亚洲AV麻豆欧美玲| 秋霞三级理伦免费观看| 亚洲欧美日韩在线观看一区二区三区| 99久久久精品免费观看国产| 国产精品爽爽久久久久久竹菊 | av色天堂2018在线观看| 国产精品外围在线观看| 男人都懂www深夜免费网站| 羞羞答答的免费视频在线观看| 99精品久久久久久久| 国语自产精品一区在线视频观看 | 久久久国产精品免费A片蜜芽广| 日本男女动态图| 中国人泡妞xxxxxxxx19| 国产精品久久久久久久人人看| 男人插曲女人的视频| 亚洲乱码日产精品BD在线下载| xiao776唯美清纯| 快播官方网站| 亚洲精品AV一二三区无码| 被老头下药玩好爽| 伦理电影v男人天堂| 亚洲国产精品一区二区第一页| 囯产精品一区二区三区线| 墨西哥美女主播| 正在播放久久| 精品人妻一区二区三区视频53|