function?y=wtfusion（x1x2Nwname）

%函數功能：
%?????函數x=wtfusion（x1x2Nwname）將兩幅原圖像x1x2進行基于小波變換的圖像融合，得到融合后的圖像y
%?????近似分量采用加權平均的融合規則，各細節分量采用基于區域特性量測的融合規則
%輸入參數：
%?????x1----輸入原圖像1
%?????x2----輸入原圖像2
%?????N----小波分解的層數
%?????wname----小波基函數
%輸出參數：
%?????y----原圖像融合后得到的圖像
%-----------------------------------------------------------------%

x1=double（x1）;???????????????????%將uint8的圖像數據類型轉換成double型進行數據處理
x2=double（x2）;

?%將原圖像x1x2分別進行N層小波分解，wname為小波基函數，
?%C為各層分解系數S為各層分解系數長度（也就是大小）.
?%C的結構:c=[A（N）|H（N）|V（N）|D（N）|H（N-1）|V（N-1）|D（N-1）|H（N-2）|V（N-2）|D（N-2）|...|H（1）|V（1）|D（1）]
?%A（N）代表第N層低頻系數（可見只有第n層含有低頻信息）H（N）|V（N）|D（N）代表第N層高頻系數分別是水平垂直對角高頻
?%S（N+2行2列）的結構是儲存各層分解系數長度的即第一行是A（N）的長度（其實是A（N）的原矩陣的行數和列數）
?%第二行是H（N）|V（N）|D（N）|的長度第三行是H（N-1）|V（N-1）|D（N-1）的長度
?%倒數第二行是H（1）|V（1）|D（1）長度最后一行是X的長度（即大小）

[C1S1]=wavedec2（x1Nwname）;?
[C2S2]=wavedec2（x2Nwname）;??

A1=appcoef2（C1S1wnameN）;???????