多項式的案例
ZEMAX | 如何創建ZOS-API自定義擴展將切比雪夫多項式轉換為擴展多項式
有關切比雪夫多項式的更多信息，請參閱： 基于切比雪夫多項式的新型自由曲面的設計與實現。
對該系統進行優化和公差分析之后，在將自由曲面反射鏡的圖紙發送給制造商之前，將切比雪夫多項式轉換為擴展多項式，這樣設計的系統就可以通過計算機輔助制造方程、模具校正和注塑過程中的模具收縮補償等工具實現。
OpticStudio有內置的非球面轉換工具，但是沒有將自由曲面轉換為其他面型的工具。幸運的是，可以使用ZOS-API構建工具。
將切比雪夫多項式轉換為擴展多項式
切比雪夫多項式
切比雪夫多項式由包含X和Y的方程式表示，這使得它們作為直角正交多項式特別實用。
擴展多項式
擴展多項式和切比雪夫多項式的定義之間的主要區別是：
圓錐常數 k
多項式的每個系數包含一個歸一化因子，而切比雪夫包含 x0 和 y0
系數隨 “i ”變化
轉換
利用Mathematica和上述定義，可以計算得出每個擴展多項式的項等于一個包含切比雪夫多項式的項的方程。結果如下：
根據上述結果，生成用戶自定義擴展，它將通過讀取切比雪夫多項式表面的系數，并計算擴展多項式表面的系數來自動轉換。
用戶自定義擴展將在切比雪夫多項式之后添加具有計算出的系數的擴展多項式表面，以及包含兩個多項式表面之間矢高差的網格矢高 ( Grid Sag ) 表面。得出的結果將滿足要求。
展開 ZEMAX軟件技術應用教程：如何創建ZOS-API自定義擴展將切比雪夫多項式轉換為擴展多項式
有關切比雪夫多項式的更多信息，請參閱： 基于切比雪夫多項式的新型自由曲面的設計與實現。
對該系統進行優化和公差分析之后，在將自由曲面反射鏡的圖紙發送給制造商之前，將切比雪夫多項式轉換為擴展多項式，這樣設計的系統就可以通過計算機輔助制造方程、模具校正和注塑過程中的模具收縮補償等工具實現。
OpticStudio有內置的非球面轉換工具，但是沒有將自由曲面轉換為其他面型的工具。幸運的是，可以使用ZOS-API構建工具。
將切比雪夫多項式轉換為擴展多項式
切比雪夫多項式
切比雪夫多項式由包含X和Y的方程式表示，這使得它們作為直角正交多項式特別實用。
擴展多項式
擴展多項式和切比雪夫多項式的定義之間的主要區別是：
圓錐常數 k
多項式的每個系數包含一個歸一化因子，而切比雪夫包含 x0 和 y0
系數隨 “i ”變化
轉換
利用Mathematica和上述定義，可以計算得出每個擴展多項式的項等于一個包含切比雪夫多項式的項的方程。結果如下：
根據上述結果，生成用戶自定義擴展，它將通過讀取切比雪夫多項式表面的系數，并計算擴展多項式表面的系數來自動轉換。
用戶自定義擴展將在切比雪夫多項式之后添加具有計算出的系數的擴展多項式表面，以及包含兩個多項式表面之間矢高差的網格矢高 ( Grid Sag ) 表面。得出的結果將滿足要求。
展開 用 Mathematica 求解多項式
多項式是由一組常數系數，a、b、c、……（數值）確定的。
TableForm[{a x + b, a x^2 + b x + c, a x^3 + b x^2 + c x + d, ". . ."}] // T
多項式求解問題就是找到一個值 x，使這些項的總和等于 0. 根據 x 的最高次數分別稱為線性、二次、三次、四次、五次、六次、七次、八次...... 多項式。我們稱 y = a x + b 為線性，是因為它的圖線是一條直線. 比如令 a = 2，b = 3，
Plot[2 x + 3, {x, -2, 1}]
2 x + 3 = 0 的（唯一）解是 x = -3/2. 一般情況下，有 x = -b/a. 由于含有 x 的平方項，y = a x^2 + b x + c 是二次的. 你會記得一元二次方程有兩個通解：
Solve[a x^2 + b x + c == 0, x]
這樣的表達式被稱為不盡根式. 最常見的應用是在幾何上. 圓、拋物線和雙曲線通常由二次多項式指定。當我們想知道一個二次多項式與已知直線何時相交時，我們就得到一個二次方程.
展開 用 Mathematica 求解多項式
多項式是由一組常數系數，a、b、c、……（數值）確定的。
TableForm[{a x + b, a x^2 + b x + c, a x^3 + b x^2 + c x + d, ". . ."}] // TraditionalForm
多項式求解問題就是找到一個值 x，使這些項的總和等于 0. 根據 x 的最高次數分別稱為線性、二次、三次、四次、五次、六次、七次、八次...... 多項式。我們稱 y = a x + b 為線性，是因為它的圖線是一條直線. 比如令 a = 2，b = 3，
Plot[2 x + 3, {x, -2, 1}]
2 x + 3 = 0 的（唯一）解是 x = -3/2. 一般情況下，有 x = -b/a. 由于含有 x 的平方項，y = a x^2 + b x + c 是二次的. 你會記得一元二次方程有兩個通解：
Solve[a x^2 + b x + c == 0, x]
這樣的表達式被稱為不盡根式. 最常見的應用是在幾何上. 圓、拋物線和雙曲線通常由二次多項式指定。當我們想知道一個二次多項式與已知直線何時相交時，我們就得到一個二次方程.
展開 
基于Mathematica的多項式系數的位數之和的并行計算
本文所介紹的算法可以幫助您探索某種NP問題蘊含在多項式內的潛在規律。
多項式是一種常見的數學對象，多項式系數（又稱組合數）是指多項式的n次方的展開式的各項系數，如下圖所示，這是一個3項式的3次方的展開式。
有多種算法可以計算多項式系數（組合數）之和，但本文要求解的問題是求這個展開式的各項系數的位數之和（十進制），這就不得不把每一項系數的位數都計算出來，對于上圖所示的多項式來說，因為展開式的每一項的系數都是一位十進制整數，所以位數之和是10，如下圖所示，可以簡單地用幾個函數在不到1毫秒的時間內計算出來。
在這個簡潔的問題中，蘊含了一個NP問題，由于所有NP問題都可以互相轉換，故不再贅述是哪一個NP問題，只需知道本文要求解的問題目前還沒有多項式級別的時間復雜度的求解算法，這也是NP（非確定性多項式）問題的含義。
若要求解7項式的30次方的多項式系數的位數之和（以下簡稱：位數之和），計算時間就來到了20秒，如下圖所示：
考慮到多項式的結構的對稱性，可以計算出系數相同的項的個數，再乘以對應的系數的位數，經過一系列計算就能在20毫秒內得到位數之和，比上圖所示的快了近1000倍。這種算法將本文的問題提煉成一個NP問題，復雜度取決于有多少種不同的系數，對于每一種不同的系數，算法只需要一次計算就能求解出這種系數的位數之和與帶有這種系數的項的個數。通過觀察多項式展開式的每一項的冪，易知，對于3項式的3次方的展開式的項來說，只有三種系數。
但這只不過是把階乘級的復雜度降低到指數級，而且這種算法的空間復雜度與其時間復雜度一樣高！非常消耗內存，因為要窮舉每一種系數。
展開 ZEMAX | 如何使用 Zernike 多項式模擬黑盒光學系統
OpticStudio 將計算系統的波前，并匹配對應的澤尼克多項式。
波前的采樣率和澤尼克多項式的項數可以在參數設置菜單中設置。波前采樣和多項式系數的關鍵參數是 RMS 擬合誤差（ RMS fit error ）和最大擬合誤差（Maxium fit error ）。在本例中使用默認采樣和多項式項數參數得到如下結果：
這表示當我們在澤尼克系數表示的波前上減去實際波前時，其殘留的誤差在百萬分之一個波前量級。這已經非常接近了！然而在實際應用時您需要對波前采樣率和多項式系數進行調整以保證多項式較高的匹配度。
我們現在需要將澤尼克多項式系數轉移到近軸等效系統中。我們可以輸出澤尼克數據并將其輸入到近軸系統中，但這一過程非常繁瑣。我們可以使用宏來完成這一操作。
下面這個宏（包含于示例文件中）名為 Zernike Readout.zpl ，它可以提取系統的澤尼克數據并保存成 .DAT 格式的文件，該格式文件可以在近軸等效系統中的表面屬性(Surface Properties)->導入( Import )->導入數據文件( Import Data File )中導入。宏執行該操作的過程如下所示：
首先定義所有需要的變量。
（需要注意的是，ZPL 中設置的采樣率和最大澤尼克項數應該與您在澤尼克分析中設置的數據相同）隨后，宏會提取出瞳直徑和澤尼克數據：
注意澤尼克表面的歸一化半徑表示出瞳直徑的一半。
展開 澤尼克多項式鏡面擬合
有會使用澤尼克多項式進行鏡面擬合的嗎？請有償指導
六自由度機械臂五次多項式軌跡規劃（Fivejtraj_Function自編寫函數） ￥56
，定義插補次數n,根據驅動器支持的最大插補次數設定
Time=[0,0,0,0,0,0;5,5,5,5,5,5;10,10,10,10,10,10];
StopTime=Time(3,1);
FixedStep=0.2;
SimpleNum=0:FixedStep:StopTime;Cunt=length(SimpleNum);
for i=1:n
%% 求每個關節角的五次多項式插補軌跡點
[q(:,i),qd(:,i),qdd(:,i)] = Fivejtraj_Function(Theta(:,i),Time(:,i), ...
ZEMAX | 在 OpticStudio 中使用自由曲面進行設計
第一類的切比雪夫多項式如下：
前十個切比雪夫多項式參數如下：
T0 (x) = 1
T1 (x) = x
T2 (x) = 2x2 - 1
T3 (x) = 4x3 - 3x
T4 (x) = 8x4 - 8x2 + 1
T5 (x) = 16x5- 20x3 + 5x
T6 (x) = 32x6 - 48x4 + 18x2 - 1
T7 (x) = 64x7 - 112x5 + 56x3 - 7x
T8 (x) = 128x8 - 256x6 + 160x4 - 32x2 + 1
T9 (x) = 256x9 - 576x7 + 432x5 - 120x3 + 9x
T10 (x) = 512x10 - 1280x8 + 1120x6 - 400x4 + 50x2 - 1
我們可以使用下方的tij（x，y）關系，將上式簡化為二維切比雪夫多項式：
將切比雪夫多項式的有限項次進行加總，產生如下的自由曲面矢高公式：
這里的aij代表多項式總和的系數；x和y則是歸一化表面的坐標；N和M則是多項式在x和y方向上的最大項次；c則代表基本球體的曲率半徑，球體位于多項式的最上方。
由于多項式的最大階數僅受到運算速度的限制，因此就一般的情況而言，10階的多項式就已經能充分的描述大部分的表面了。
使用切比雪夫多項式表面設計離軸拋物線
我們使用切比雪夫多項式產生離軸的拋物面作為一個簡單的示例。系統的規格如下方的圖1。
展開 Ansys Zemax | 在 OpticStudio 中使用自由曲面進行設計
第一類的切比雪夫多項式如下：
前十個切比雪夫多項式參數如下：
T0（x） = 1
T1（x） = x
T2（x） = 2x2- 1
T3（x） = 4x3- 3倍
T4（x） = 8x4- 8倍2+ 1
T5（x） = 16x5- 20倍3+ 5 倍
T6（x） = 32x6- 48倍4+ 18 倍2- 1
T7（x） = 64x7- 112倍5+ 56 倍3- 7倍
T8（x） = 128x8- 256倍6+ 160 倍4- 32倍2+ 1
T9（x） = 256x9- 576倍7+ 432 倍5- 120倍3+ 9 倍
T10（x） = 512x10- 1280倍8+ 1120 倍6- 400倍4+ 50 倍2- 1
我們可以使用下方的tij （x，y）關系，將上式簡化為二維切比雪夫多項式：
將切比雪夫多項式的有限項次進行加總，產生如下的自由曲面矢高公式：
這里的一個ij 代表多項式總和的系數; x和y則是歸一化表面的坐標; N和M則是多項式在x和y方向上的最大項次; c則代表基本球體的曲率半徑，球體位于多項式的最上方。
由于多項式的最大階數僅受到運算速度的限制，因此就一般的情況而言，10階的多項式就已經能充分的描述大部分的表面了。
使用切比雪夫多項式表面設計離軸拋物線
我們使用切比雪夫多項式產生離軸的拋物面作為一個簡單的范例。系統的規格如下方的圖 1。
展開 Ansys Zemax | 在 OpticStudio 中使用自由曲面進行設計
第一類的切比雪夫多項式如下：
前十個切比雪夫多項式參數如下：
T0（x） = 1
T1（x） = x
T2（x） = 2x2- 1
T3（x） = 4x3- 3倍
T4（x） = 8x4- 8倍2+ 1
T5（x） = 16x5- 20倍3+ 5 倍
T6（x） = 32x6- 48倍4+ 18 倍2- 1
T7（x） = 64x7- 112倍5+ 56 倍3- 7倍
T8（x） = 128x8- 256倍6+ 160 倍4- 32倍2+ 1
T9（x） = 256x9- 576倍7+ 432 倍5- 120倍3+ 9 倍
T10（x） = 512x10- 1280倍8+ 1120 倍6- 400倍4+ 50 倍2- 1
我們可以使用下方的tij （x，y）關系，將上式簡化為二維切比雪夫多項式：
將切比雪夫多項式的有限項次進行加總，產生如下的自由曲面矢高公式：
這里的一個ij 代表多項式總和的系數; x和y則是歸一化表面的坐標; N和M則是多項式在x和y方向上的最大項次; c則代表基本球體的曲率半徑，球體位于多項式的最上方。
由于多項式的最大階數僅受到運算速度的限制，因此就一般的情況而言，10階的多項式就已經能充分的描述大部分的表面了。
使用切比雪夫多項式表面設計離軸拋物線
我們使用切比雪夫多項式產生離軸的拋物面作為一個簡單的范例。系統的規格如下方的圖 1。
展開 
知識分享 | 應變四分之一橋溫度補償的理論與實踐
在catman中進行熱補償
在某些情況下，多項式還包括在溫度變化時影響應變信號的因素：
引線的影響(εl)：在這個具體示例中，還添加了測量不確定度和連接到測量柵絲應變片引線的影響。一般來說，必須考慮引線的影響，但在應變片類型和制造商之間可能會有所不同。如果您使用的是我們 HBM 3- 或 4-線專利技術，可補償所有電纜電阻的影響，但在某些情況下，2線部分不能自動補償。
測量不確定度 (εu)：測量不確定度是總計算中應考慮的一般部分。甚至多項式也有一些導致了不確定性的分散項。
適用的多項式如下：
在本例中，多項式如下：
假設應變測量期間的溫度恒定為100°C（T=100°C），引線長度為10 mm（L=10 mm）。請注意引線的長度可能導致結果會有所不同。從多項式的角度看，它表明熱應變對結果有很大的影響，因為它高于100μm/m。
為了計算熱輸出應變，將溫度和引線長度部署到多項式中：
結果與數據表上顯示的多項式非常吻合。通過多項式，考慮了最顯著的影響，而只有非常長的2線電纜才會產生額外的影響。因此我們需要使用3線和4線技術來補償電纜電阻的影響。
將熱應變疊加在測量應變上，即可得到應變校正值，該值僅考慮機械應變：
02 應變系數調整
大多數制造商的多項式是用k=2的系數（或k系數）來測量的，而實際測量中，k系數通常是不同的。這種效應主要在極端條件下出現，如高溫和低溫或高應變。
展開 網格階數詳解：高階網格生成
回歸分析被用于確定多項式模型或者等效樣條模型的系數然后用插值方法給出兩端點間的數據點，并將這些數據點賦予曲線網格以符合多項式模型。高階網格生成要將類似的過程用于線性網格（不管是結構化網格還是混合網網格），以便提取與多項式曲線相關的連續點多項式曲線。
讓我們來看下面的圖示，線性網格被用于描述有曲度變化的渦輪葉片表面。在對線性網格設定邊界條件后，利用算法將線性網格的節點與多項式曲線匹配關聯。CFD 工程師可以自行選擇最適合模擬需求的多項式網格階數。針對一些具有特殊多項式曲率的曲面，生成的多項式曲線網格也可以很好的符合葉片表面的曲度變化，且不需要線性網格那樣高密度網格節點分布。
完成相關表面的多項式曲線定義后，可以用插值法高效生成任意密度的網格。網格的精度可以通過調整插值后的網格密度或不同的插值方法來進一步優化。下圖左可見插值后高階網格的示例。下圖右可以看到一些插值法可能在生成的插值網格中產生偽影，所以選擇正確的插值方法也是生成高精確曲線網格的關鍵。
插值后的多項式曲線網格與插值法導致偽影的線性網格
Cadence Pointwise 網格生成工具可以幫助 CFD 工程師創建復雜幾何模型高精度模擬所需要的高階網格，且不會顯著增加計算復雜性。
文章來源Cadence楷登PCB及封裝資源中心
展開 MATLAB的數據結構
function p=plus(a,b)
a=polynom(a); b=polynom(b);
k=length(b.c)-length(a.c);
p=polynom([zeros(1,k) a.c]+[zeros(1,-k) b.c]);
同理，還可以重載定義多項式的減法運算：
function p=minus(a,b)
a=polynom(a); b=polynom(b);
k=length(b.c)-length(a.c);
p=polynom([zeros(1,k) a.c]-[zeros(1,-k) b.c]);
乘法運算：多項式的乘法實際上可以表示為系數向量的卷積，可以由 conv() 函數直接獲得。 故可以如下重載定義多項式的乘法運算。
function p=mtimes(a,b)
a=polynom(a); b=polynom(b); p=polynom(conv(a.c,b.c));
乘方運算： 多項式的乘方運算只限于正整數乘方的運算，其 n 次方相當于將該多項式自乘 n 次。若 n=0,則結果為 1。 這樣我們就可以重載定義多項式的乘方運算為：
function p=mpower(a,n)
if n>=0, n=floor(n); a=polynom(a); p=1;
if n>=1,
for i=1:n, p=p*a; end
end
else, error('Power should be a non-negative integer.')
end
多項式求值問題：可以對多項式求值函數 polyval() 進行重載定義。
function y=polyval(a,x)
a=polynom(a); y=polyval(a.c,x);
定義了此類之后，我們就可以方便地進行多項式處理了。
展開 ZEMAX | 旋轉對稱不規則性（RSI）簡介
誤差由Zernike多項式表示，具有三階球差和高階球差的形式。
RSI多項式形式
RSI的多項式形式是：
根據OpticStudio中定義的Zernike標準系數，RSI涉及的Zernike多項式為：
為什么RSI是一個問題？
在任何具有大量表面的系統（例如內窺鏡、投影透鏡或光刻透鏡）中，RSI很快就會成為問題。通常假設此類系統中的表面誤差是隨機組合的，因此對總誤差的估計是每種誤差類型的和均方根 (RSS)。但這不適用于 RSI 多項式，RSI誤差不直接相加，但總誤差大于RSS計算得出的值。例如，四個RSI多項式項如下所示，每個項的系數值為0.01。
上述四個多項式之和如下所示。多項式的RSS為 0.1，但真正的表面誤差為0.14。
請注意，對于上圖，RSI多項式可用于捕獲表面上的卷邊。卷邊是常見的制造缺陷。
根據ISO 10110在圖紙上指定RSI
RSI的定義在ISO 10110第5部分：表面形狀公差中給出。
光學圖紙上使用的最基本的RSI形式是：
其參數定義如下:??
A 表面的 P-V 功率誤差，以環數或納米為單位
B 總 P-V 表面不規則度，以環數或納米為單位
C P-V RSI，以環數或納米為單位
按照ISO10110-5的規定，這種基本形式有許多變體。可以指定測試波長。還可以包括對表面總RMS的誤差限制。
展開 