2012 一月號 Engineering Structures

人生中的第一篇SCI文章終於排定於2012年的一月首發了!

期間的甘苦等待算是走完了一個行程，下一個行程看可不

可以再半年內將工不二字出個頭呢? 要好好排定進度了!

但手邊的事還真多!! @#$%^&*......

看樣子是要花一點功夫小心答覆了新工作的近對應退了!

最佳化有感!

最近因為需要使用MATLAB來跑位移的最佳化，但是對MATLAB的概念根本是小兒科，除了一堆基本的問題沒搞清楚外更別提要如何做最佳化了!:(

但是看了半天，終於釐清最佳化的頭緒了，第一個就是需要目標方程式，第二個就是限制方程式，最後就是要"我猜我猜我猜猜了!"也就是初始值的GUESS!

因為我在看的一篇論文其中他使用了FMINCON的SOLVER方法，他說了兩個準則。

1) the ratios of the measured damage displacements from the column database to the calculated damage displacements had a mean value equal to 1.0

2) the coefficient of variation (COV) of the ratios was minimized.

我從兩年前就看過這句話，但是看很明白但是如何用?就像武俠小說中的口訣，背得很熟，但是如何使勁發招，完全不靠譜呀!!! 更別提如何用最佳化，因為根本就搞不懂。但是人就是如此卑賤，被逼到梁山了，不得不發狠，於是好好找了一堆文件，看了半天，沒用就是沒用!不得已只好發問了原始作者，沒答!只好發問題於網路:

問題如下:

若已知一資料庫如下:

Y a b c
.4 1 3 4
.6 2 4 7
.9 3 9 3
1 4 3 3

回歸方程式假設如下:

y=x1*(1+x2*a)*(1+x3*b)^-2*(1+x4*c)^1.5

條件式需求如下:
平均的Y需等於計算後的平均的y
計算後的變異係數需最小(minimize)
其中的x1~x4為待求之係數

試問 利用MATLAB中的optimization 的中fmincon目標方程式與條件方程式該如何撰寫?

在中國的小木蟲與知識+發問，小木蟲的回答不切題，台灣的延長了一次，心想唉!完了，天要滅我!要轉路了

沒想到善心人士-->Elisha，彷如白衣大士來臨回了文，當時我覺得不會又不切題唄!

一看，怪怪回得有模有樣!

ANS:

主要有兩個m檔, 一個為你的程式, 一個為fmincon的執行程式

-----------------------------------------

function f = model_1(x)

x

yg = mean([0.4, 0.6, 0.9, 1.0]);

a = [1, 2, 3, 4];

b = [3, 4, 9, 3];

c = [4, 7, 3, 3];

ys = mean(x(1)*(1 + x(2)*a).*(1 + x(3)*b).^-2.*(1 + x(4)*c).^1.5);

f = abs(yg - ys);

-----------------------------------------

clear all

clc

% Initial Guess

x0 = [0; 0; 0; 0];

A = []; b = [];

Aeq = []; beq = [];

lb = [-10;-10;-10;-10];

ub = [10;10;10;10];

nonlcon = [];

[x,fval,exitflag] = fmincon(@(x) model_1(x),x0,A,b,Aeq,beq,lb,ub,nonlcon);

-----------------------------------------

我找出來的x如下所示

x =

1.7471

-0.1154

0.3006

0.2143

誤差值

fval =

7.1505e-007

他點醒我兩件事，一是我當初自己也寫過但是他明白告訴我當時我寫的方式是大錯，應該用行列式的格式，我寫成矩陣，笨!第二個是目標方程式的寫法是很單純的，而不是y=x1*(1+x2*a)*(1+x3*b)^-2*(1+x4*c)^1.5，蠢!

經過嘗試後，發現該答案仍無法達到我的需要，但是卻點醒我格式的基本問題。於是又在苦思兩天，忽然清晨中腦海靈光一現，恍然大悟!!!!這老美那兩點點出了核心，但未察覺，這個死老美，害我想破頭，原來:

2) the coefficient of variation (COV) of the ratios was minimized.-->這個就是目標方程式

1) the ratios of the measured damage displacements from the column database to the calculated damage displacements had a mean value equal to 1.0-->這個就是限制方程式

於是將ELISHA的內容改寫，重新再跑! YES! 跑出來的不僅合理且也符合標準!

但也發現老美的問題

(1)根本不相等! (騙我)

(2)COV也不是最小化! (又騙我)

但是經過此一訓練，我終於明瞭一些事，就是刻骨銘心才是硬功夫!!!!!!

答案應該如下:

x1=	0.1948613
x2=	0.9824801
x3=	-0.0092008
x4=	-0.0035009

The value of COV is very close to zero.

期刊論文接受了!!

終於第一篇的SCI的文章，在今天晚上(2011/9/9,20:17)主編寄信通知正式接受了! 謝謝葉勇凱博士的耐心與黃世建教授的提攜，與當時一起辛苦該試驗的人員(當然包括學妹"盧錦嫻")!!，也謝謝阿正的討論與九儀的英文協助!真的謝謝您們!!
題目:Out-of-Plane Seismic Behavior and CFRP Retrofitting of RC Frames Infilled with Brick Walls.

不過這條路還真漫長，摘要過程:

2010/5/30 投稿

2010/6/1 投稿編號

2010/6/4 主編確定 (陳錦順)

2010/8/5 退回第一次大修

2010/10/9 修完上傳

2011/1/21 退回第二次大修

2011/2/21 修完上傳

....

....

....

漫長的等待

....

等到沒感覺了

....

....

2011/9/6 退回小小修改 GOOD! 如荒漠甘霖!

2011/9/9 修完上傳

2011/9/9 接受，YEAH!!!!

心跳聲音可以聽的到!

鬱積了半年，竟然在今天晚上，出現了既怕又等了半天的
關鍵寄信人的的字眼-->EES ，心中竟然心跳加速到自己
都可以聽到 :(，只有自己控制自己的呼吸，慢慢端詳他的判決
-->Accept pending minor changes(小修後接受)，
真是令我雀躍!!!!!!
這等於是博士畢業的重要門檻的至大關鍵!!!!!!
好了! 後面要繼續加油了!!高興一晚就好了! 繼續苦工中，

工不青蛙繼續加油唄!

引自聯合報:

http://udn.com/NEWS/NATIONAL/NATS1/6353063.shtml

若是台九線是珍珠的話，台一線還真是個平民呢!同樣都是從行政院出發，但是最後與台九線相接於楓港，但是台一線的風景呢? 沒有，它只有默默地提供台灣大多數人行的方便，有機會可以從台一線連到台九線的環島旅程中，見識到台灣這塊土體的繁華與美景!

微軟新注音 2003!

話說輸入法應該要有學習的能力，但微軟新注音2003好像智障還是腦袋有洞，怎麼學也不容易學!氣得我在寫論文時又氣又急! 尤其在打"位一十" -->位移時! 唉! 不說了!

下面的內容轉自網路:

這是網路上某個女孩子的網誌，是真實的笑話喔~

因為它造成我與朋友對話時不少的誤會

你們看看……這是我某天與一個新朋友在網路上的對話

『夜來襄，妳說妳教音樂，是哪方面？』朋友問

『就平時交肛琴，禮拜六日在教會工作…』我答

『妳..妳..是說教鋼琴嗎?』對方有點遲疑

『對阿…』我說

『喔..那~在教會是在做什麼工作阿？』對方又問

『就袋屍斑，領屍之類的…』我回答

『對….對不起…..妳…妳是說…帶……帶詩班…領詩是嗎?』對方語氣有點怪怪的

(唉呦妖壽…怎麼沒檢查就送出，幸好對方對教會還有點認知)

『呵呵~對阿對阿~不好意思打錯字…』我很不好意思地回

『哈….沒關係的，我懂。那妳在教會工作很已經久囉？』對方又問

『其實我是這一年半才開始焦屍斑的…….以前是做行政方面的事』我仔細算算這樣告訴他

『這…….』對方(見鬼了…….這是葬儀社輸入法嗎?)

『呵呵~我是說這一年半開始教詩班，以前是當幹事…..』我趕緊解釋

『沒..沒關係….我了解。那~聊聊妳的興趣吧，妳平時喜歡做什麼？』對方問

『我啊？最多時候就看VCD囉！』我說

『這樣啊？我也是耶！那你會喜歡韓劇嗎？最近很熱門喔，還是比較喜歡西洋的？』對方問

『含巨根吸陽具一樣插進，我都不喜歡。』我回

(啊啊啊…..來不及了….又按ENTER太快了….)

『對不起…我是說…韓劇跟西洋劇一樣差勁，我都不喜歡啦…』我連忙解釋企圖挽回

(我看現在代誌大條了，我的慾女形象毀於一旦，又來了….是玉女才對)

『…………….』對方

( 我想…..他已經嚇昏了吧)

抖掉包袱，可不是相聲的抖包袱!

摘自網路:

這個世界上有太多煩惱，所以有許多人不約而同地問了佛祖一樣的問題：「我該怎麼做，才能不再煩憂？」
佛祖給的答案都相同：「只要放下，你就能不再煩惱。」
有個自以為聰明的人很不服氣，便專程去找佛祖，挑釁地問：「世上有千千萬萬個人，就有千千萬萬種煩惱。但是您給他們的解決方式都完全相同，那豈不是太可笑了？」
佛祖沒有生氣，只是反問男子：「你晚上睡覺的時候，會做夢嗎？」
「當然會！」男子回答。
「那麼，你每天晚上做的夢，都是一樣的嗎？」佛祖又問。
「當然是不一樣的。」
「你睡了千千萬萬次，就做了千千萬萬個夢。」佛祖微笑地說：「但是要結束夢的方法，卻都是一樣的。那就是『醒過來』！」
男子聽到佛祖的回答，啞口無言。
我很喜歡這則禪宗故事。這則故事把煩惱比喻成「夢」，只要我們願意放下，就能從夢中清醒來。無論你的煩惱是什麼，方法都是一樣的。
我當義工的時候，分別探訪過兩個老婦人。
她們的人生際遇有驚人的雷同：都自小喪母，成為別人的養女；都嫁了一個酗酒的老公；兒子都很不孝順；兩人身體都不好，晚年都在洗腎。然而，她們兩人的心境卻大不相同。
第一個老婦人談及往事，她處之泰然。她說種種的不順遂，是她人生的功課，現在她年紀大了，該吃的苦都吃過了，所以她覺得自己真好運，人生像倒吃甘蔗，愈來愈順利。
這名老婦人，現在是某公益團體的義工。
第二個老婦人談及往事，仍然咬牙切齒。她指天怪地，罵兒子詛咒媳婦，她說老天爺對她真不公平！
這名老婦人，也跟那個公益團體有關，但她不是義工，而是受到義工關懷的憂鬱症患者。
讓我印象深刻的，是她們兩個臉上的表情、散發的氣息形成了強烈的對比。
當下讓我感覺：中國人說的「相由心生」，真的非常有道理。
當義工的老太太，因為已經放下了煩惱，所以給人的感覺是祥和又親切，讓人很想親近，表情有如菩薩；得憂鬱症的老婦人，還把多年來的不順遂、痛苦扛在肩上，她的氣息尖銳又暴烈，讓人避之唯恐不及，表情如同厲鬼。
然而，對於後者，其實我是深刻地同情的。她的一輩子都在吃苦，但她不但不願意放下那些苦，還要抱著那些其實已經遠離的苦難不放，讓它們繼續折磨自己。
也許，人生真的有太多無法逆轉的苦難、無法挽回的遺憾。
也許在夜闌人靜的時候，我們會因此哭泣，覺得蒼天真是無情，竟然讓我承受如此巨大的痛苦……
但是，你我都不應該忘記，能將我們從痛苦的地獄拯救出來的絕對不是別人，只有我們自己。
人生有太多包袱，我們可以選擇將它一一扛在肩上，也能選擇瀟灑地放下它們，而能否放下，則會決定我們未來的人生將繼續痛苦，還是可以享有幸福。
想擺脫惡夢的糾纏，唯一的方式是「醒過來」；想忘卻人生的不順遂，唯一的方法就是「放下」。
人生有太多包袱，我們可以選擇將它一一扛在肩上，也能瀟灑地放下它們。

少將共諜案之核心問題

　羅賢哲被吸收，不必怪中共，這是它情報工作該做的，全世界也都這樣做，我們該檢討的是他潛伏七年為何沒發現，最後還是靠美方提供情資破獲。其實台灣的反情報未受到重視，幾達門戶洞開地步。

by 丁渝洲

我也是這麼認為的!
再者有三個問題值得當局深思
1.財產申報制的真實用處與此案產生極大的反諷!
2.我國少將的資料美國竟然一清二楚，這對於自家的情報或反情報是個警訊
3.國內情報與軍情單位的首長應該換人坐了!