發問區
會員登入
新使用者？立即註冊
.
服務首頁|服務說明|Yahoo!奇摩.知識+ 首頁
知識分類
電腦網路
科學常識
醫療保健
煩惱心事
生活資訊
手機通訊
休閒嗜好
視聽娛樂
運動體育
社會人文
商業金融
教育學習 .如何做
煩惱
主題知識
.我要發問 發表 我要發問 ..熱門： 蛇年運勢 新年賀詞 壓歲錢由來 流感 金蛇年 主題 金蛇報喜，大過好年 用手機上知識+ .知識搜尋 ....知識+ 首頁> 電腦網路> 程式設計> 其他 加入追蹤 轉寄朋友 友善列印 .知識問題|
matlab程式...20點...急. 發問者： 龍 ( 初學者 5 級)
發問時間： 2005-06-29 23:23:04
解決時間： 2005-07-04 00:28:52
解答贈點： 20 ( 共有 0 人贊助 )
回答： 2 評論： 0 意見： 0
[ 檢舉 ]
網友正面評價
100%
.共有 6 人評價.請幫我解釋每一行代表什麼意思??愈詳細愈好...
尤其是mloop的..
function [fnew,xnew]=asaq(func,x,maxstep,qf,lb,ub,tinit)
xold=x;
fold=feval(func,x);
n=length(x);
lk=n*10;
q=qf*n;
nv=log(maxstep*ones(n,1));
mv=2*ones(n,1);
c=mv.*exp(-nv/n);
t0=tinit*ones(n,1);
tk=t0;
a=lb*ones(n,1);
b=ub*ones(n,1);
k=1;
for mloop = 1:maxstep
for tempkloop=1:lk
fold=feval(func,xold);
u=rand(n,1);
y=sign(u-0.5).*tk.*((1+ones(n,1)./tk).^(abs((2*u-1))-1));
xnew=xold+y.*(b-a);
fnew=feval(func,xnew);
if fnew rand
xold=xnew;
end
end
f(mloop)=fnew;
tk=t0.*exp(-c.*k^(q/n));
k=k+1;
end
tf=tk; 2005-07-01 01:26:51 補充
for mloop = 1:maxstep
這行是表示 mloop 為1到maxstep的所有數
這我知道...我要知道裡面的內容...@@
最佳解答發問者自選
.. 回答者： quy1qk ( 初學者 5 級 )
回答時間： 2005-07-02 20:11:53
[ 檢舉 ] .function [fnew,xnew]=asaq(func,x,maxstep,qf,lb,ub,tinit) %宣告一個function,名為asaq
% -----輸入-------
% func 為字串,代表所要執行的函數
% x 為矩陣
% qf 為數字
% maxstep 為矩陣
% tinit 為數字
% lb 為數字
% ub 為數字
% -----輸出-------
% xnew 為矩陣
% xnew 為矩陣
xold=x; % x 放到 xold
fold=feval(func,x);
% 利用 feval 指令,將 x 代入 func(字串) 函數,把計算結果儲存於 fold
n=length(x); % 取得矩陣 x 的行數,放到 n
lk=n*10; % 將 n*10 放到 lk
q=qf*n; % 將數字 qf*n 結果放到 q
nv=log(maxstep*ones(n,1));
% ones(n,1) 產生一個 n*1 的矩陣,其內容皆為1
% 將 maxstep 矩陣 * ones(n,1) 矩陣,將運算後的矩陣內的 元素 全部取 log 放到 nv
mv=2*ones(n,1); % 將 ones(n,1) 矩陣內的 元素 全部 *2 放到 mv
c=mv.*exp(-nv/n);
% exp(-nv/n) 將 nv 矩陣內的 元素 全部 * -1/n 再將每個 元素 取指數
% 將 n*1 的矩陣 mv 乘以 exp(-nv/n) 產生的 n*1 矩陣的 相同位置 的元素
% ex.mv = [ 3;3;3], exp(-nv/n) = [ 0.5 ; 0.5;0.5]
% mv.*exp(-nv/n)=[ 1.5 ; 1.5 ; 1.5 ] #有"."跟沒"."會差很多
% 再將結果儲存到 c 矩陣
t0=tinit*ones(n,1); % 將 ones(n,1) 矩陣內的 元素 全部 * 數字 tinit 放到 t0 矩陣
tk=t0; % 將 t0 矩陣給 tk 矩陣
a=lb*ones(n,1); % 將 ones(n,1) 矩陣內的 元素 全部 * 數字 lb 放到 a 矩陣
b=ub*ones(n,1); % 將 ones(n,1) 矩陣內的 元素 全部 * 數字 ub 放到 b 矩陣
k=1; % 將 k 的內容放 1
for mloop = 1:maxstep % mloop 由 1 到 maxstep 矩陣中 "最小" 的元素
for tempkloop=1:lk % tempkloop 由 1 到 數字 lk
fold=feval(func,xold); % 利用 feval 指令,將 xold 代入 func(字串) 函數,把計算結果儲存於 fold
u=rand(n,1); % rand(n,1) 產生一個 n*1 的亂數矩陣,其元素為 0 - 1 之間的隨機數字,儲存到 u 矩陣
y=sign(u-0.5).*tk.*((1+ones(n,1)./tk).^(abs((2*u-1))-1));
% (abs((2*u-1))-1) 將 (2*u-1) 矩陣內的元素全部取 絕對值 ,並且 全部 元素 -1
% (1+ones(n,1)./tk) 將 ones(n,1) 矩陣內的元素 相同位置 除於 tk 矩陣,然後全部元素 +1
% (1+ones(n,1)./tk).^(abs((2*u-1))-1) 將 (1+ones(n,1)./tk) 矩陣內的元素 取 (abs((2*u-1))-1) 矩陣內 相同位置 元素 的次方
% sign(n) 函數: n > 0 , sign=1 ; n = fnew 矩陣 相同位置 的元素 則 條件成立
xold=xnew; % 將 xnew 給 xold 矩陣
elseif exp((fold-fnew)/norm(tk))>rand
% norm(tk) 為一函數 得出數字
% exp((fold-fnew)/norm(tk)) 將矩陣(fold-fnew)內元素全部 / norm(tk) 數字,再求出每個元素的exp
% 將 exp((fold-fnew)/norm(tk)) 矩陣的每個元素跟 0-1 範圍的隨機亂數比較
xold=xnew; % 將 xnew 給 xold 矩陣
end
end
f(mloop)=fnew; % 給值
tk=t0.*exp(-c.*k^(q/n)); % 運算
k=k+1; % k=k+1
end
tf=tk; % 給值
% 註: 有些變數可能不是矩陣,只是這裡我用矩陣解釋而已
% 檔案 http://hk.geocities.com/quy1qk/asaq.zip
% flashget 下載同鏡向~~
參考資料
腦袋 + 經驗 相關詞：
矩陣相乘,矩陣運算,bcg矩陣,矩陣乘法,轉置矩陣,matlab 矩陣,矩陣行列式,矩陣組織,矩陣公式,管理矩陣
矩陣,元素,數字,位置,相同,matlab,flashget 下載,函數,儲存,abs[ 快速連結 ] 其它回答( 1 ) | 意見( 0 ) | 評論( 0 )
.發問者評價 謝謝你啦... .發表你的評價
你的評價 發表評價： 正面 普通 負面
評價內容： 發表 取消
.
加入追蹤 轉寄朋友 友善列印 .馬上按讚 加入 Yahoo! 奇摩 知識+ 粉絲團
•免費索取商英光碟+小書
•多益700分線上測驗題庫
•立即免費測試你的多益等級
•英文email超實用金句！
•測你在旁人眼中的英文力？
•多益700分線上測驗題庫
相關問答
[ C&C++ ]MATLAB矩陣問題
.
[ 數學 ]矩陣有重復的特徵值則無法對角化
.
[ C&C++ ]matlab矩陣的問題
.
[ 其他 ]請問 什麼是矩陣
.
[ 其他 ]矩陣的除法
.
[ 數學 ]有關向量和矩陣的問題!!!
.
更多
.其他回答(1)
意見(0)
相關評論(0)
.
001
回答者： ~小鐘~ ( 初學者 5 級 )
回答時間： 2005-06-30 01:01:27
[ 檢舉 ]
..for mloop = 1:maxstep
這行是表示 mloop 為1到maxstep的所有數 1
目前沒有資料 發表意見發表意見字數已達上限，要改成發表評論嗎？. 發表 取消 .
目前沒有資料 我要評論 註冊
會員登入 .公告： 知識團員轉粉絲全數完成 .
HOT! 拍賣 | 奇美19吋 Vista電腦 熱門3C .
.刊登贊助網站•優質半導體元件代理商-上巨電子
www.sb-ele.com.tw
專營Maxim、Dallas、DS等半導體IC元件代理，規格齊全，歡迎洽詢！
www.sb-ele.com.tw
•全民免費學電腦 8大熱門職訓課
www.pccenter.com.tw
優惠僅到月底！學承職訓電腦課程免費學，快速提升職場競爭力，立即搶先申請！
www.pccenter.com.tw
•免費線上電腦課程 再送職涯分析
www.lccnet.com.tw
搶先索取免費5000堂線上電腦課程，聯成再提供熱門廠商媒合，一次找到好工作！
www.lccnet.com.tw
•巨匠電腦八大主題課程免費學到飽
www.pcschool.com.tw
簡單三步驟，巨匠八大主題電腦課程讓你選，再送好禮6選1。免費線上課程學到飽。
www.pcschool.com.tw
•IT資訊人才聯合徵才培訓申請中
www.pcschool.com.tw
響應政府資訊人才培訓計畫，APP職缺大增，求職必備電腦證照，快線上申請培訓。
www.pcschool.com.tw
•GOMAJI軟體進修團購超值價
www.gomaji.com
學技能也要省！揪團學時下超夯電腦課程，GOMAJI要你省上加省！
www.gomaji.com.最新其他
發問中
已解決
.如何用QB語言寫一個簡型的文書軟體環境?
FPGA程式
如何在UBUNTU12.10下安裝舊版MYSQL 和GCC
更多
系統外包的實務問題 更多
精選關鍵字
..ASP.net 命令 函數 分析 Push 陣列 C++ RFID VB.net 原始碼 遊戲 型別 TABLE PHP XML 程式 8051 C 參數 TELNET 宣告 設計 字串 VB DOS 組合語言 ASP 條件 HTML 連結
.知識搜尋 ...雅虎資訊 版權所有 (c) 2013 Yahoo! Taiwan. All Rights Reserved. 「本服務設有管理員」
服務條款隱私權政策..知識+ 之問答內容是由參與Yahoo!奇摩知識+ 之網友提供，僅供參考，Yahoo!奇摩不保證其正確性。
...
.
|
phibrainno15 發表在 痞客邦 留言(0) 人氣()
留言列表
留言列表
