C# 定積分求周長&面積原理 代碼實現

来源:https://www.cnblogs.com/u3ddjw/archive/2018/05/01/8945782.html
-Advertisement-
Play Games

前言: 前些日子,因為工作原因,接觸到了求解曲線周長,真的是搞了很久,學生時代真的很簡單,但是如今的我來說,忘記了....很多人跟我應該一樣。 所以來鞏固加強一下記憶。一開始的時候,求周長嘛,找公式唄,什麼matlab呀,亂七八糟的,暈,最後找到了可能還不能滿足項目的需求,因為可能計算量過大。(我就 ...


前言:

         前些日子,因為工作原因,接觸到了求解曲線周長,真的是搞了很久,學生時代真的很簡單,但是如今的我來說,忘記了....很多人跟我應該一樣。

          所以來鞏固加強一下記憶。一開始的時候,求周長嘛,找公式唄,什麼matlab呀,亂七八糟的,暈,最後找到了可能還不能滿足項目的需求,因為可能計算量過大。(我就是這樣子的,靈活性相對較低)

          還有就是明明自己可以用代碼實現,為什麼非要插件,工具吶,這麼不自信的?

          所以,“一怒之下”,自己去看了一下定積分求周長的原理,自己還是用代碼來實現吧。

(以下內容純是個人這段時間的理解,如果有錯誤的,歡迎指正出來。)

 

首先需要說說兩個概念,曲線和周長,因為我們要求他們嘛。

曲線:

  這個世界,有曲線嗎?我的回答是,沒有。那...這...曲線是由無數個直接拼接而成。再準確的說無數個很短的曲線拼接而成。

(如果您徹底理解了這句話,後面就不用看了,基本就沒了。)

 

面積:

  與周長的概念類似,沒有正方形,沒有圓形。只有三角形,所有的圖形都是三角形拼接而成。而兩個三角形拼成長方形,而我們的面積是由無數個長方形,拼接而成。

  

                (定積分原理的參考圖)

 

源碼實現: 

double GetLength(float start,float end)
        {
            double sumLength = 0;

            float eachX = (end - start) / testCount; 
            for (int i = 1; i < testCount;i++ )
            {
              
                double curY =ArcFunction(start+eachX*i);
                double previousY = ArcFunction(start+eachX*(i-1));
                //根據c²=a²+b²
                double curLength = Math.Sqrt(Math.Pow(eachX, 2) + Math.Pow(curY - previousY, 2));
                sumLength += curLength;
            }
                return sumLength;
        }


   解釋:

             testCount,即自定義的測試數量,可以理解為精細度,值越大,計算量越大,數據越準確,這個可看你項目需求精細度,通過該變數

可在實現最少的計算量情況,實現你要的效果。

    eachX,就是你的曲線被分成N份,每份的長度。

    curY,當前點的y軸分量

    previousY,上一個點的y軸分量,

    curLength,即如圖   

    

                                                      

                  (剩餘的部分,代碼裡面含解釋,個人喜歡放在源碼裡面,原生的,純24k原創)   

 1     int testCount = 1000; //所謂的測試細緻度吧,可動態調控,你自己掌握。
 2         /// <summary>
 3         /// 通過已知周長,獲取x軸的分量
 4         /// </summary>
 5         /// <param name="length"></param>
 6         /// <returns></returns>
 7         double GetRateXByLength(double length)
 8         {
 9             float eachX = 1.0f;
10             for (int i = 1; i < testCount; i++)
11             {
12                 double curY = ArcFunction(eachX * i);
13                 double previousY = ArcFunction(eachX * (i - 1));
14                 double curLength = Math.Sqrt(Math.Pow(eachX, 2) + Math.Pow(curY - previousY, 2));
15                 length -= curLength;
16                 if(length<=0)
17                 {
18                     return i * eachX;
19                 }
20             }
21             return testCount * eachX;
22         }
23             
24         
25       
26 
27 
28         /// <summary>
29         /// start到end範圍內的面積
30         /// </summary>
31         /// <param name="start"></param>
32         /// <param name="end"></param>
33         /// <returns></returns>
34         double GetArea(float start,float end)
35         {
36             double sumAera = 0;
37 
38             float eachX = (end - start) / testCount;
39             for (int i = 1; i < testCount; i++)
40             {
41                 double curY = ArcFunction(start + eachX * i);
42                 //面積 = 長*寬
43                 double curAera = curY * eachX;
44                 sumAera += curAera;
45             }
46             return sumAera;
47         }
48 
49 
50 
51         /// <summary>
52         /// 通過已知面積,獲取x軸分量
53         /// </summary>
54         /// <param name="aera"></param>
55         /// <returns></returns>
56         double GetRateXByAera(double aera)
57         {
58             float eachX = 1.0f;
59             for (int i = 1; i < testCount; i++)
60             {
61                 double curY = ArcFunction(eachX * i);
62                 double curAera = curY * eachX;
63                 aera -= curAera;
64                 if(aera<=0)
65                 {
66                     return i * eachX;
67                 }
68             }
69             return testCount*eachX;
70         }
71 
72 
73 
74         /// <summary>
75         /// 通過x分量,得出y的值。(好像意義不大,但是好像可能有些人不是很理解,寫給某些人看的,一目瞭然)
76         /// </summary>
77         /// <param name="x"></param>
78         /// <returns></returns>
79         double GetYByX(float x)
80         {
81             return ArcFunction(x);
82         }
83 
84 
85         /// <summary>
86         /// 核心控制函數。
87         /// </summary>
88         /// <param name="x"></param>
89         /// <returns></returns>
90         double ArcFunction(float x)
91         {
92             return Math.Pow(x, 2); //這邊我用冪函數來測試。各位爺可以換其他函數啊。
93 
94             //註:如果對曲線靈活性要求很高,推薦使用貝塞爾曲線。
95             //詳情可參考:
96         }

 

 

補充:這邊,我對貝塞爾曲線做一下補充吧。因為我因為什麼冪函數,指數函數遇到的肯,因為這些函數畢竟還不是那麼靈活,都具有一定“規律”。

貝塞爾曲線的靈活,受6個參數控制,三個點嘛。(二維空間)

 詳情可參考這篇:http://www.sohu.com/a/118656687_466876


 


您的分享是我們最大的動力!

-Advertisement-
Play Games
更多相關文章
  • 二分查找也稱為折半查找,是對有序元素查找的一種演算法,在查找的過程中,不斷的將搜索長度減半,因此效率不錯。Java的JDK提供了二分法查找的演算法,使用的方法是Arrays.binarySearch()。binarySearch()方法提供了多種數據類型的二分查找,比如實現了int、float、doub ...
  • 1 學習計劃 1、在realm中進行授權 2、使用shiro的方法註解方式許可權控制 n 在spring文件中配置開啟shiro註解支持 n 在Action方法上使用註解 3、★使用shiro的標簽進行許可權控制 n 在頁面引入shiro的標簽庫 n 在頁面中使用標簽 4、總結shiro提供的許可權控制方 ...
  • 本文主要介紹MyBatis映射mapper.xml中,輸入映射和輸出映射、動態sql、關聯查詢等的寫法。 ...
  • 1.1 常見對象及含義 1.2 worksheets方法一覽 1.3 事件分類 1.4 字元處理函數 實例 sub strcomp運算() msgbox strcomp("ABCD","abcd",1) '參數1以文本比較方式,不區分大小寫,返回0,表示相等 msgbox strcomp("ABCD ...
  • DozerBeanMapper是JavaBean的映射工具,可以進行對象之間相同屬性名賦值 關於PO、DTO、VO在分層模型之間的關係:首先在持久層由DAO訪問資料庫將數據對象封裝成PO,然後在service層將PO映射成DTO對象併進行進一步處理,處理完之後DTO傳到表現層,接著在表現層映射成VO ...
  • ThreadLocal 稱為線程本地存儲,它為每一個使用它的線程提供一個其值(value)的副本。可以將 ThreadLocal<T> 理解成 Map<Thread, T>,即使用當前線程為 key 的一個 Map,ThreadLocal 的 get() 方法從 Map 里取本地變數(本地 valu ...
  • EF Core (EntityFramework Core)是實體關係映射(O/RM)資料庫訪問框架。這個模式的好處就是讓開發人員可以用對象模型來操作資料庫,這是一種對開發人員較為友好的方式。O/RM牛逼在哪?非對象模型操作資料庫,就要直接用sql語言來做大量的CURD操作(Creat 創建;upd... ...
  • StackExchange.Redis 擴展,更簡單的泛型操作,並提供一些的適用於業務場景中的擴展 ...
一周排行
    -Advertisement-
    Play Games
  • 示例項目結構 在 Visual Studio 中創建一個 WinForms 應用程式後,項目結構如下所示: MyWinFormsApp/ │ ├───Properties/ │ └───Settings.settings │ ├───bin/ │ ├───Debug/ │ └───Release/ ...
  • [STAThread] 特性用於需要與 COM 組件交互的應用程式,尤其是依賴單線程模型(如 Windows Forms 應用程式)的組件。在 STA 模式下,線程擁有自己的消息迴圈,這對於處理用戶界面和某些 COM 組件是必要的。 [STAThread] static void Main(stri ...
  • 在WinForm中使用全局異常捕獲處理 在WinForm應用程式中,全局異常捕獲是確保程式穩定性的關鍵。通過在Program類的Main方法中設置全局異常處理,可以有效地捕獲並處理未預見的異常,從而避免程式崩潰。 註冊全局異常事件 [STAThread] static void Main() { / ...
  • 前言 給大家推薦一款開源的 Winform 控制項庫,可以幫助我們開發更加美觀、漂亮的 WinForm 界面。 項目介紹 SunnyUI.NET 是一個基於 .NET Framework 4.0+、.NET 6、.NET 7 和 .NET 8 的 WinForm 開源控制項庫,同時也提供了工具類庫、擴展 ...
  • 說明 該文章是屬於OverallAuth2.0系列文章,每周更新一篇該系列文章(從0到1完成系統開發)。 該系統文章,我會儘量說的非常詳細,做到不管新手、老手都能看懂。 說明:OverallAuth2.0 是一個簡單、易懂、功能強大的許可權+可視化流程管理系統。 有興趣的朋友,請關註我吧(*^▽^*) ...
  • 一、下載安裝 1.下載git 必須先下載並安裝git,再TortoiseGit下載安裝 git安裝參考教程:https://blog.csdn.net/mukes/article/details/115693833 2.TortoiseGit下載與安裝 TortoiseGit,Git客戶端,32/6 ...
  • 前言 在項目開發過程中,理解數據結構和演算法如同掌握蓋房子的秘訣。演算法不僅能幫助我們編寫高效、優質的代碼,還能解決項目中遇到的各種難題。 給大家推薦一個支持C#的開源免費、新手友好的數據結構與演算法入門教程:Hello演算法。 項目介紹 《Hello Algo》是一本開源免費、新手友好的數據結構與演算法入門 ...
  • 1.生成單個Proto.bat內容 @rem Copyright 2016, Google Inc. @rem All rights reserved. @rem @rem Redistribution and use in source and binary forms, with or with ...
  • 一:背景 1. 講故事 前段時間有位朋友找到我,說他的窗體程式在客戶這邊出現了卡死,讓我幫忙看下怎麼回事?dump也生成了,既然有dump了那就上 windbg 分析吧。 二:WinDbg 分析 1. 為什麼會卡死 窗體程式的卡死,入口門檻很低,後續往下分析就不一定了,不管怎麼說先用 !clrsta ...
  • 前言 人工智慧時代,人臉識別技術已成為安全驗證、身份識別和用戶交互的關鍵工具。 給大家推薦一款.NET 開源提供了強大的人臉識別 API,工具不僅易於集成,還具備高效處理能力。 本文將介紹一款如何利用這些API,為我們的項目添加智能識別的亮點。 項目介紹 GitHub 上擁有 1.2k 星標的 C# ...