hs演算法

㈠ 請問有人知道共軛梯度法的FR,PRP,HS三個演算法的Matlab程序嗎

%共軛梯度法 FR
% G為對稱正定矩陣，X是初始點，e為精度
%a是精確線搜索步長
function [m2,a,d,X,g1,f1] = conjgrad(G,b,c,X,e)
n=length(G);
if n==2
format long e %rat
syms x1 x2
f=1/2*[x1,x2]*G*[x1;x2]+b'*[x1;x2]+c;
g=[diff(f,x1);diff(f,x2)];
g1=subs(subs(g,x1,X(1,1)),x2,X(2,1));
d=-g1;
a=-(d'*g1)/(d'*G*d);% a=-((X(:,1)'*G*d+b'*d)/(d'*G*d)); a=g1(:,1)'*g1(:,1)/(d(:,1)'*G*d(:,1));
X(:,2)=X(:,1)+a*d;
g1=[g1 subs(subs(g,x1,X(1,2)),x2,X(2,2))];
m1=norm(g1(:,1));m2=norm(g1(:,2));
k=(m2/m1)^2;
i=2;
while m2>=e
d(:,i)=-g1(:,i)+k(i-1)*d(:,i-1);
a(i)=-(d(:,i)'*g1(:,i))/(d(:,i)'*G*d(:,i));
%a1(i)=-((X(:,i)'*G*d(:,i)+b'*d(:,i))/(d(:,i)'*G*d(:,i)));a(i)=g1(:,i)'*g1(:,i)/(d(:,i)'*G*d(:,i));
X(:,i+1)=X(:,i)+a(i)*d(:,i);
g1=[g1 subs(subs(g,x1,X(1,i+1)),x2,X(2,i+1))];
m1=m2;m2=norm(g1(:,i+1));
k(i)=(m2/m1)^2;
i=i+1;
end
f1=subs(subs(f,x1,X(1,i)),x2,X(2,i));
elseif n==3
format long
syms x1 x2 x3
f=1/2*[x1,x2,x3]*G*[x1;x2;x3]+b'*[x1;x2;x3]+c;
g=[diff(f,x1);diff(f,x2);diff(f,x3)];
g1=subs(subs(subs(g,x1,X(1,1)),x2,X(2,1)),x3,X(3,1));
d=-g1;
a=-((X(:,1)'*G*d+b'*d)/(d'*G*d));
X(:,2)=X(:,1)+a*d;
g1=[g1 subs(subs(subs(g,x1,X(1,2)),x2,X(2,2)),x3,X(3,2))];
k=(norm(g1(:,2))/norm(g1(:,1)))^2;
m=norm(g1(:,2));
i=2;
while m>=e
d(:,i)=-g1(:,i)+k*d(:,i-1);
a(i)=-((X(:,i)'*G*d(:,i)+b'*d(:,i))/(d(:,i)'*G*d(:,i)));
X(:,i+1)=X(:,i)+a(i)*d(:,i);
g1=[g1 subs(subs(subs(g,x1,X(1,i+1)),x2,X(2,i+1)),x3,X(3,i+1))];
k=(norm(g1(:,i+1))/norm(g1(:,i)))^2;
m=norm(g1(:,i+1));
i=i+1;
end
f1=subs(subs(subs(f,x1,X(1,i)),x2,X(2,i)),x3,X(3,i));
end

㈡ 主曲線的演算法研究

本程序定義了主曲線和確定主曲線的實際演算法。多邊形線演算法的基本運演算法則是首先確定一條直線段，然後在循環演算法中通過不斷加入新的頂點來增加線段的數量。在加入一個新的頂點以後，所有的頂點位置在一個內部的環中被更新。由演算法所產生的曲線如圖1，在這個例子中，PL運演算法則和HS運演算法則在計算的復雜程度和演算法的完善程度上作出了比較。4段和15段，由在半圓上任意兩個坐標點加入單獨的高斯誤差而產生。
PL演算法[12]在由NIST19號專有資料庫產生的單獨數據元構成的圖像中得到了測試。我們發現PL演算法可以有效的找出沒有環和分叉的圖像的中間軸。這個在圖2中有顯示。由於中間軸可能是一些曲線連接而成而不是只有一條曲線，所以在這里我們擴展了PL演算法，找出數據元的主曲線。擴展了的演算法也包含了實現分段線性骨架的兩個原則，一種獲取字元圖像近似輪廓的初始化方法和一系列用來改善由初始化方法獲得的骨架結構質量的更改結構工作。為了避免混淆，我們用術語「骨架」來表示近似性中間軸的二元圖像，我們把由PL演算法產生出的連接曲線看做模板骨架.
應用
主曲線以前被用在圖像處理領域。這種圖像用來描述在衛星拍攝下的冰川及其輪廓。其主程序用了一個閉合的曲線來估算冰川的輪廓。專家們除去了HS演算法[13]中不合適的部分，並且用更完善的粗略估算冰川輪廓的程序來取代HS演算法的初始化步驟。此外，採用數據元集合，而不是HS演算法所產生的點或線的集合，向中間軸集中的方式來擴展現有的聚合演算法。初始化的曲線由SOM演算法[12]的一個變數產生，在SOM演算法中相互關系被定義為字元圖像的一個最小二叉樹。HS演算法用來使曲線與模板相對應。下一步的要點與SOM演算法的更新規則類似。
利用主曲線實現分段線性骨架的方法被Mahmoud、Datta和Parui[14]等人所提出。同樣的方法，在SOM演算法中用來最優化分段線性骨架的頂點位置。演算法在構建骨架方面採用「自頂向下」的策略：近似性地從一個線性拓撲出發，然後逐步加入環和交叉到骨架中，骨架是由基於SOM演算法的當前幾何樣式得出的。SOM演算法涉及到一個獲取字元圖像分段線性骨架的運演算法則。這種運演算法則以ISODATA演算法[12]為基礎，ISODATA演算法近似於SOM演算法。
目的
主曲線演算法的目的是找出與字元圖像相對應的光滑的分段線性曲線。這些曲線在某個頂點彼此連接，因而在字元圖像的中心平面范圍內形成一個歐幾里德曲線。一個歐幾里德曲線G（V，S）在空間中由變數V和S確定，
主曲線演算法從一個基於傳統稀釋方法的初始化工作開始。初始曲線找出字元圖像的近似拓撲結構，然而,它不是平滑的，而且它通常包含許多假的分叉和不適合的結構元素。為了解決這兩個問題，在前兩步中間增加了更改結構的步驟(圖3)使得主曲線演算法更加有效。在更改結構步驟中，我們運用一些操作手段來調整曲線骨架結構的不完整的地方。(a)圖是在初始化步驟後由主曲線演算法生成的曲線骨架；（b）圖是經過首次擬合-光滑步驟後生成的曲線骨架；(c)圖是經過更改結構後生成的曲線骨架；(d)圖是第二次擬合-光滑步驟後產生的曲線骨架(演算法輸出)。我們重復使用PL演算法的擴展版本來構造光滑的骨架曲線，同時保持曲線與字元圖像的輪廓的距離近似相等。演算法建立在最小能量函數的基礎之上
研究動機與意義
自1904年Spearman[13]提出線性主成分分析方法以來，由於這種方法簡單且便於使用，至今還是數據統計分析的重要工具之一。線性主成分分析的原理是將數據集合投影到一個矢量，使得投影的均方差最大，由此，將這個矢量稱為數據集合的第一主成分。正是這個考慮，在均方差的意義下，這個方法有兩個重要的優點：其一，數據集合可以使用一個矢量來描述，從而達到減小信息描述長度的目的，其二，計算第一以及依次主成分，可以變換為求解基於數據自相關矩陣的特徵值方程。另外，第一與依次主成分矢量保持正交關系，這意味著，與主成分矢量平行的矢量具有與主成分相同的性質。正是這兩個原因，加上在統計上以均方差為保證，主成分分析得到廣泛的應用。 由於信息描述長度與信息保持性之間存在矛盾，相對較長的信息描述長度，較短描述長度的信息描述是以損失信息保持性為代價的，而主成分分析的本質是一種在均方差意義下的統計平均。盡管它可以獲得較短的信息描述長度，但是，信息保持性的代價相當大，由於信息保持性是對數據分布的規律性認識，因此，對某些問題，這個代價可接受的，但是，另外一些問題，可能就不能接受了。例如，對原始語音信號的分析，單純的主成分分析就不一定有效。 為了說明信息描述長度與信息保持性之間的關系，下圖是一個簡單的例子。圖5是由300個包含誤差的數據點構成的餘弦狀分布，圖5(a)的直線是數據的第一主成分線，其對餘弦數據的描述長度顯然較圖5(b)要短，因為這些數據點將使一個線段描述，因此，只需給出線段起點和終點即可，可以認為圖5(a)給出了一個短描述長度的關於數據集合的描述；顯然，圖5(b)對數據的信息保持性則比圖5(a)要好，一方面，它與數據間的距離均方差也要小，另一方面，它勾畫出原始信息的輪廓。圖5(b)恰恰是本文所討論的根據主曲線原理所獲得的結果。那麼，兩種描述哪一個更為合理呢？顯然，這沒有一個一般性的答案，這取決於所需解決問題的需求。
圖5 信息描述長度與信息保持之間的關系
圖5也說明無監督學習較監督學習困難的原因，問題本身的病態定義導致不得不引入復雜性思想，如統計學習理論中的風險結構最小、貝葉斯學派中的貝葉斯信息准則、Kolmogrov演算法[11]復雜度引出的最小描述長度等等，來尋求在信息保持性與數據描述長度之間的折衷。目前，盡管在主曲線的研究中，還存在著大量的數學問題，但是，由於其暗示的廣泛應用前景，已引起計算機科學家的關注，現在應用方面已有大量的報道，如線性對撞機中對電子束運行軌跡的控制、圖像處理中辨識冰原輪廓、手寫體的主曲線模板化、語音識別中對聲音數據的約簡建模和數據可聽化、生態學中尋找種群的有序分布及過程監控。同時，在認知領域Seung[14]提出感知以流形方式存在，並在神經生理學上發現整個神經細胞群的觸發率可以由少量的變數組成的函數來描述，如眼的角度和頭的方向，這隱含了神經元群體活動性是由其內在的低維結構所控制。主曲線本身是流形的一維形式。
原理、發展脈絡以及應用
為了說明主曲線的原理、發展脈絡以及應用，首先介紹其原始動機是必要的。Hastie在他關於主曲線的開創性論文中描述了其研究動機，Hastie認為主曲線與非線性回歸方法的發展上具有對稱性，分別是線性主成分分析與線性回歸分析的非線性推廣模型，Hastie寫到：類似於線性回歸方法的非線性推廣，使用光滑曲線代替線性主成分總結數據，以求出對稱變數之間的曲線，這樣的曲線從數據的中部光滑地通過，且不限制於對數據的光滑線性平均，甚至不限制於數據的中部是直線，只希望使得數據點集合到曲線的正交距離最小。這個陳述清晰地指出了他的研究與主成分分析和回歸分析的區別，並給出了建模的統計目標。同時，從這個陳述中也可以看出，基於直線的主成分分析只是主曲線的一個特例。因此，主曲線的提出，可以理解為基於線性的主成分分析方法向更精確描述實際世界的非線性分析方法的推廣。 應該指出的是，目前，「向非線性」推廣是數據統計分析的研究主流，但是，存在著不同的技術路線。
分類
最典型的研究大致可以分為兩類：
其一，根據統計學習理論中的核技術，將數據集合映射到特徵空間，然後，在特徵空間計算數據集合的主成分，稱為核主成分分析(KPCA)，這個技術路線的本質還是線性主成分分析。
其二，從數據本身的分布出發，希望找到能最好描述數據內在結構的概率分布模型，如生成式拓撲映射(Generative Topographic Mapping,縮寫為GTM)，矢量量化(VQ)及主曲線，以及流形擬合等。提出「主曲線的研究與回歸分析有何區別」是自然的，兩者的動機都是企望求出一個函數代替數據集合，但是，兩者有本質的差別：其一，一般地說，回歸分析是假設數據集合的變數有因果關系，換句話說，數據的變數可以表示為一個因果關系的函數，有些是自變數，有些則是因變數。其二，回歸分析一般需要給定一個數學基函數，回歸分析是根據數據集合中變數的因果關系，計算這個數學基函數待定的參數。這與主曲線的研究動機完全不同。對於前者，主曲線的研究目標是解決數據變數沒有必然因果關系的問題，更重要的是後者，認為事先假定數據服從某種分布(如正態分布)的方法(這與給定數學基函數，根據數據確定參數的方法類似)，對某些未知世界的解釋是不合理的，因為這個假設可能是錯誤的，為了保證數據分析不是假設在一個錯誤結構下，主曲線的研究強調非參數分析。當然，這並不是完全否定了參數法這一經典方法，事實上，參數法在先驗知識明確時存在相當大的好處。總之，根據上述討論的動機，主曲線是尋找一種幾何上直觀、理論上完備、就解決的問題而言，這個方法與主成分分析一致，但是，主曲線的目標不是僅僅尋找一個數據集合的脊樑骨，而是希望獲得數據集合的骨架。數據集合的脊樑骨可以使用線性的方法解決，但骨架就需要藉助非線性方法了。應該強調的是，主曲線繼承了主成分分析的眾多思想，它是主成分分析的非線性推廣。另外，盡管這個方法似乎與回歸分析的目標類似，都是試圖獲得對數據集合信息長度更短的表示，但是，這個方法與回歸分析完全不同，最大的差別是它不是事先給定一個基函數(或假定一個分布)，從而將問題變換為參數估計問題，而是一種非參數的方法

㈢ 求樓梯鋼筋的計算方法

樓梯鋼筋按情型應當分為：板下端負彎筋、板下端負彎筋(直錨式)、板上端負彎筋、板上端負彎筋(直錨式)、梯板分布筋、板下端負彎筋分布筋、板上端負彎筋分布筋、梯板面貫通筋。

其數量計算公式為：(BLT-2*CZ)/S+1（其中BLT：梯段板寬；CZ主筋保護層厚度；S：箍筋、板、牆鋼筋按間距排列布置的分布間距）

其長度計算公式為：L/4*sqrt(BS*BS+HS*HS)/BS+(T-CZ)+LA+2*WG180（其中BS：踏步寬；HS：踏步高度；T：梯板厚度；LA：鋼筋錨固長度；L：梯段水平長；sqrt：平方根函數）可直接在excel中將公式中的數據填寫完整，計算。

(3)hs演算法擴展閱讀

樓梯鋼筋布置方法：

一、先在土建基層上放線，確定立柱固定點的位置，然後在樓梯地面上用沖擊鑽鑽孔，再安裝膨脹螺栓，螺栓保持足夠的長度，在螺栓定位以後，將螺栓擰緊同時將螺母與螺桿間焊死，防止螺母與鋼板松動。扶手與牆體面的連接也同樣採取上述方法。

二、放線。由於上述後加埋件施工，有可能產生誤差，因此，在立柱安裝之前，應重新放線，以確定埋板位置與焊接立桿的准確性，如有偏差，及時修正。應保證不銹鋼立柱全部座落在鋼板上，並且四周能夠焊接。

三、安裝立柱焊接立柱時，需雙人配合，一個扶住鋼管使其保持垂直，在焊接時不能晃動，另一人施焊，要四周施焊，並應符合焊接規范。

四、扶手與立柱連接。立柱在安裝前，通過拉長線放線，根據樓梯的傾斜角度及所用扶手的圓度，在其上端加工出凹槽。然後把扶手直接放入立柱凹槽中，從一端向另一端順次點焊安裝，相鄰扶手安裝對接准確，接縫嚴密。相鄰鋼筋對接好後，將接縫用不銹鋼焊條進行焊接。焊接前，必須將沿焊縫每邊30―50的尺寸。