小佳的天地: 11月 2012

2012年11月22日星期四

101年度-北科大-有限元素法期中考題

解：

1. Discrefize fhe problem info nodes and elemenls.

Element	Node i	Node j	L(m)	θ
(1)	1	2	4	0
(2)	2	4	2	90
(3)	1	4	4.47	26.57
(4)	3	4	4	0

2. Assume a solution that approximates the behavior of an element.

3. Develop equations for elements

4. Assemble elements.

5. Apply the boundary conditions and loads.

7. Obtain other information.

善意像顆球丟出還會再回來^^

有一次家人一起到一家西餐廳用餐。
那天客人滿多的，菜出得慢，不過還好，陸續送了來，大家開始用餐，吃到甜點
時，我點的餐還沒送來，因為是一家人，所以我吃一點你的前餐，你吃一點他的主
菜，大家都沒注意到這事，但我沒忘記。
送飲料來時，服務小姐很不好意思地詢問，是不是還有一份餐未送來？
我說是，服務小姐很客氣的說了一聲對不起就離開了。
不到三十秒，老闆娘帶著大廚到我旁邊，連聲抱歉，說剛把主菜放下去而已，還要一
二十分鐘，問我能不能等一下，或者要退掉。
我答：「沒關係，我知道你們今天忙，難免會慢一些，我等一下好了，等好的時候幫
我打包起來，我帶回家。」
老闆娘跟大廚連聲道謝離開。
然後，原來沒點餐的小朋友都有了甜點、水果與飲料。
埋單時我驚訝地問說是否算錯了，老闆娘在旁邊解釋道：「因為你的諒解與客氣，所
以餐點打八折，沒上桌的那一份免費，小孩子的附餐也免費招待。」
我笑笑說：「妳太客氣了。」

老闆娘回了一句話：「因為你的客氣，所以不得不讓我們更客氣。」
我笑著離開，不因少花錢，只因客氣也可以傳染給別人……
我的工作讓我常常有機會介紹想裝潢的客戶給做室內設計的朋友，按照行規，或多或
少總會有些介紹費，但我從來沒接受。
大概是從高中開始，當朋友要回報我對他們的幫助時，我總是拒絕。

我只認真的告訴朋友，那天我需要幫助，拉我一把就好了。
因此，當我須全家外出離家數天時，我不用擔心家裡的魚會餓死、花會枯死，當我需
要搬運東西時，我不用找搬家公司﹔當我無車可用時，不必擔心沒人載……
善意與幫助像是一顆球，當你毫不遲疑的將球丟向對你招手的人後；有天，當你也在
招手時，也會有顆球飛到你手中，或者，還不只一顆。
即使你的善意被誤解與踩在地上了，有一天，老天也會找機會補償你的!
喜歡這篇文章嗎....
如果喜歡 ~~~ ，就把它當顆球讓別人也感受吧！

狗尾續貂：只是這顆球有時彈性好，因此彈回來得快些；有時它的彈性卻又好像差了一點，因此彈回來得慢一些。有時候大概是因為擲錯了地方吧，球就沒再彈回來了。此時要告誡自己：只要是顆球，除非它是用金屬做成的，否則一定都有彈性，我們不要對「球一定有彈性」這個法則懷疑或失去信心！擲出去的球沒再彈回來，可能是因為「彈性疲乏」的緣故吧？名為「善意」的這顆球，價格並不貴，我們就再花點小錢買一顆新的回來，在往後的一生中，繼續作個擲出「善意球」的投手吧。
*************************************************************************
*************************************************************************

--
This message has been scanned for viruses and
dangerous content by MailScanner, and is
believed to be clean.

2012年11月21日星期三

臺灣二、三、四、五年級生的成長史！

臺灣二、三、四、五年級生的成長史！

~~~換句話說，也就是民國20至50年間出生的呆胞成長史！~~~

作者：四年級生

靠著這一群呆胞，臺灣才有穩定、安寧的今天！

靠著這一群呆胞，臺灣才有今天的好日子過！

首先，提醒大家：

台灣沒石油、沒鐵、沒煤、沒黃金、沒牧場、沒資源、沒邦交。

可是今天：

台灣GNP21000美元
GNP 購買力37500美元
外匯存底世界前4名
電子業世界前4名
鋼鐵業世界前12名
石化業世界前12名
進出口世界前12名
高鐵捷運品質世界前5名
全民健保品質世界前5名
水費世界最低
電費世界最低
汽油世界最低
學費世界最低
繳稅率世界最低
罵總統世界最大聲
反對黨世界最兇悍

這全都是因為當年二三四五年級生當年
都很乖，都很認真聽講總統的訓話，都很努力讀書、都很努力考試的結果！
每週上班五天半 薪水比現在少 學生比現在多

不會也不敢跳樓

不知道如何罷工只會努力工作

不敢丟豬糞

不可能去總統府前種稻米

只會聽鳳飛飛白嘉莉唱愛國正經歌！

然後，讓你我重溫下列的史實：

「二、三、四、五年級生」：

★童年時期只有三家電視台可以選擇，台語節目每天限播時數與時段，布袋戲甚至要「說國語」。

★為了政令宣導，九點鐘到九點半必須是「三台聯播」，轉過來轉過去都是「寒流」「海棠血淚」。堅信大陸同胞都活在「水深火熱」之中，都在啃樹皮（長大之後第一次在新加坡某研討會遇到一位大陸同胞，才聽說他們的教科書也說台灣同胞在啃樹皮。樹皮好無辜∼）

★小孩時期都曾經被叔叔阿姨舅舅姑姑帶去冰果室跟情人約會充當電燈泡。

★晚飯之後搶不到電視，因為阿嬤都在看楊麗花歌仔戲。

★最紅的卡通片應該是[太空飛鼠]「雷鳥神機隊」「小英的故事」「小甜甜」「科學小飛俠」「無敵鐵金剛」. ★最愛的漫畫書是[漫畫大王]「小叮噹」「老夫子」「四郎與真平」，第一次投稿是「王子月刊」和「國語日報」。

★「包青天」「保鏢」演好久，好比「夜市人生」和「娘家」。

★很期待中秋節，因為可以蒐集中秋月餅的玻璃紙。

★三不五時就有「反共義士」駕機投奔自由，可以領黃金還能上電視綜藝節目。

★一二三自由日，要聽谷正綱演講，但是聽不懂說什麼。

★大多數的人都知道「宋能爾牧師佈道大會」。

★ 四月四日 兒童節會收到健素糖當禮物，長大之後才知道那是給豬吃的。

★地理課本歷史課本關於台灣的篇幅可能只有半頁，一定要背黃河跟長江流經哪些省分，但不知道濁水溪在哪裡？

★小時候最大的恐懼就是「共匪打過來」！

★在學校講台語要被罰錢（侯文詠比較倒楣，他要掛牌子！）

★ 10月31日 會發動全國募捐，籌建「介壽號」戰機向偉大的蔣公祝壽

★蔣總統死掉的時候，電視機會從彩色變黑白，學生制服要別黑紗，升旗典禮要「默哀三分鐘」，從此畢業旅行都要去慈湖跟頭寮謁靈。

★但「畢業旅行」不可以是「畢業旅行」，而是「公民訓練活動」。

★每到國定假日就要穿制服去體育館開「團結自強大會」，不知道是什麼道理，但也是男女校互窺的好時機。

★演講比賽不管什麼題目，最後都要「復興中華文化，堅守民主陣容，消滅萬惡共匪，拯救大陸同胞」「以三民主義統一中國」「把青天白日滿地紅國旗插滿秋海棠的每一個角落」等口號做結尾，這是最早的內容自動產生器嗎？

★選舉賄選最流行送味精與彩色鍋。

★最冷的笑話是：「請問下一屆蔣總統是誰？」

★寫作文寫到「總統」，前面沒有空格就死定了！

★朝會聽演講要是聽到「蔣總統」三個字就要立正！

★超過五成的人會唱兩種版本的「總統蔣公紀念哥」 （噗，寫錯啦，是紀念歌……還有，忘了空一格，警總不要來抓我）

★高中男女在校外走在一起，被教官發現一定記過！

★制服沒有紮進褲子裡也要記過！

★高中女生頭髮長度不得超過耳下一公分或不能超過衣領（所以脖子短的人就倒楣了），要穿軍訓裙，要戴船型帽，軍歌比賽要踢木蘭步，每學期打靶還會被後座力撞到內傷

★大學很難考，俗稱「窄門」，還有一部賣座電影叫做「門裡門外」。

★大學在校外偷開舞會如果遇到警察臨檢，必須開燈關音響，立刻變成包餃子或湯圓聚會。

★幾乎每個人都玩過小天使遊戲。

★暑假不能出國遊學，只好參加救國團活動，因為不是國民黨員，所以報不到最熱門的「溪阿縱走」和「虎嘯戰鬥營」。

★每個班級都有國民黨指派的小組長。

★要看藝術片就要去小隔間的MTV，想約會就去西門町冰宮，要改褲子就去中華商場。

★最時髦的手機是黑金剛。

★有BBCall就很了不起，用數字就可以發簡訊，譬如520就是我愛你！

★看「天天開心」看到主題曲一字不漏記下來。

★剛畢業的時候，平均有四個工作機會可以選擇（好囂張）！

★領到第一份薪水的年頭，定存利率還有10％（神話∼∼）

★見識過股市破萬點的榮景，而且破很多（人生少數驕傲）

★所謂懂電腦的人，就是學過COBOL、BASIC、Fortran。沒有PC，沒有Internet，只有半夜登記排隊上機房。3.5吋磁碟片還沒出現，用那種軟軟薄薄，容量 1.2M 的5.25吋軟碟。

★剛開始上班的時候，使用過Telex電報機，用過手動打字機（打錯字要用刀片刮。什麼意思？ 唉，聽不懂就算了！）

★沒有ATM，沒有號碼牌，提款存款都要在櫃臺前面排隊。

★星期六要上半天班，除夕要到五點鐘才能下班。

★看病之前要先去公司管理部開一張手寫「勞保單」。

★因為沒有手機，所以才有日劇「東京愛情故事」莉香與完治最後在車站擦身而過的淒美梗。

★沒有交友網站，只有「愛情青紅燈」與「姊妹雜誌」。

★還沒有「壹週刊」跟「蘋果日報」，最猛的是「獨家報導」跟「自立晚報」。

★以前怕人口爆炸，大力鼓吹「家庭計畫」「3321：結婚三年生一個，再過三年生一個， 兩個孩子恰恰好，男孩女孩一樣好」；現在憂心少子化與高齡化，高喊「助你好孕」,發獎金催生！

★以前要「保密防諜」，不得與匪勾結；現在要ECFA，沒有中國就沒有世界觀！

★以前聽費玉清唱「中華民國頌」好激動；現在重新把歌詞看一遍，覺得好KUSO。

★以前罵總統要被抓去關，現在罵總統有時候被鼓勵，有時候被教育部警告，但還不至於抓去關。

★以前最紅的是綜藝節目，現在最紅的是call in 節目，但兩者基本功能相同

★以前說「當選」，現在說「凍蒜」（吼，講台語，罰錢。指∼∼）

這篇故事真是生動、忠實！ 真是越看越有味啊！

如果你/妳就是這群呆胞之一，拍拍手吧！讓你/妳/我肯定自我一下吧！

最後更要謝謝寫這這篇故事的作者，讓你/妳/我都能溫故知新！

2012年11月12日星期一

matlab格式化數值顯示

格式化數值顯示
當MATLAB顯示數值的結果時，它遵守幾條規則。預設是：如果結果是整數，那麼MATLAB會以整數顯示它。同樣的，當結果是實數，MATLAB會以大約以小數點後4位小數來顯示它。如果結果中有很大的數字超過了這個範圍，MATLAB會將結果以類似科學記算機的科學符號來表示它。你可以在File選單中的preference選單項目中的Command window preferences（命令視窗偏好）中指定不同的數值格式，或是在提示符號後打上適當的MATLAB format命令來覆寫（override）這個預設行為。以特殊變數pi為例，以不同的format選項產生的數值顯示格式如下：

MATLAB命令 pi 註解
format short 3.1416 5位的數字
format long 3.14159265358979 15位的數字
format short e 3.1416e+000 5位的數字加上指數
format long e 3.141592653589793e+000 15位的數字加上指數
format short g 3.1416 以format short或format short e來佳最表示
format long g 3.14159265358979 以format long或format long e來佳最表示
format hex 400921fb54442d18 16進位的浮點數
format bank 3.14 小數點後取2位
format + + 正（＋），負（－），或是零（0）
format rat 355/113 近似的有理數
format debug Structure address = 2937df8 m = 1 n = 1pr = 1b238c08 pi = 0 3.1416 除了short g外，內部儲存的資訊

範例：原預設為顯示小數至第4位

>> K\F

ans =

0.0005

0.0010

0.0016

0.0022

0.0029

0.0036

0.0044

0.0053

修改preference選單項目中的Command window preferences（命令視窗偏好）顯示為short g後，即會顯示最佳化的數值格式

>> K\F

ans =

0.00049628

0.0010268

0.0015966

0.002212

0.0028809

0.0036135

0.0044232

0.0053282

matlab 線性聯立方程式矩陣解法

線性聯立方程式矩陣解法

所謂線性方程式係指方程式中之變數僅屬一階者，其幕次不能大於一，且任何項中不得有兩變數相乘或相除之情形。在矩陣表示法上，通常採用[A][x]=[b]之型式，其解為[x]=[A]\[b]或[x]=[b]/[A]。

對於線性聯立方程式目前已經發展出許多種解法，其中包括變數迭代消去法及克雷蒙法(Creamer's Method)。目前MATLAB 所使用之解法則以變數迭代消去法為主，或稱為高斯(Gauss Ｅlimination)消去法。這是利用兩線性方程式分別乘以某特定常數，使其與另一方程式之同一變數係數相同，因而兩式相減得以消除該變數。如此展轉消除，最後可以得聯立方程之解。茲舉例說明如下：


  3x +4y =10
  5x -2y =8

上述聯立方程式中，先將第一式乘兩邊5、第二次兩邊乘3，兩式相減，即可得：


   5(4y) - 3(-2y) = 5* 10 - 3*8

由上式整理之後，可得果為y = 1；代入原式任一式可得x = 2，終得解。如果利用第九章之繪圖指令可以繪出此兩條曲線，由其交點即可得到相同的解。


   ezplot('3*x+4*y-10',[-10,10]);hold on 
   ezplot('5*x-2*y-8',[-10,10])

圖10.1

利用這樣的解法並不是解線性方程式之重點，因為線性代數方程式可用矩陣的型式表示，甚至以矩陣之乘除法可以得解。以上面之聯立線性方程式為例，可以化成ＡＸ＝Ｃ之型式，設x=x1，y=x2：

左除法求解

利用MATLAB求解時，只要用矩陣倒除即可，或稱為左除法。例如一組聯立方程式以矩陣表示為[A][X]=[C]時，其未知數項[X]可以利用MATLAB倒除的指令求得，即[X]＝[A]\[C]：


A=[3 4;5 -2];
C=[10;8];
X=A\C 

X =
    2
    1

結果立即得到x1=2、x2=1。

反矩陣法求解

解上式矩陣有時可用反矩陣法。反矩陣有倒數的意義，但在矩陣中必須為方矩陣，且須為非特異矩陣(non-singular)。通常用(^-1)次方表示之，如A^-1 。其基本特質為與原矩陣相乘後，將得到單位矩陣Ｉ：


  A^-1A=AA^-1=I

以此乘於ＡＸ＝Ｃ之兩側，可得：


  A^-1AX=A^-1C 　或A^-1AX=IX=X= A^-1C

在MATLAB反矩陣A^-1 之求法為inv(A)；Ｉ或稱為單位矩陣，表示為eye(A)。就上面之聯立方程式之解，可以用MATLAB演算如下：


A=[3 4;5 -2];
C=[10;8];
X=inv(A)*C 
X =
   2.0000
   1.0000

所得的結果與前述之分析相同。理論上雖然如此，在演算時求反矩陣較費時間，使用左除法是
為較佳的操作法。

範例：有一組聯立方程式如下：


   6x - 3y +5z = 12
  10x + 4y -8z =-20
  -6x + 2y +3z = 15

解：先安排Ａ及Ｃ矩陣，再利用左除法求得[X]。


A=[6 -3 5;10 4 -8;-6 2 3];C=[12 -20 15]';
x=A\C
C=A*x 

x =
   0.0479
   2.1737
   3.6467
C =
  12.0000
 -20.0000
  15.0000

由上述的計算，可以迅速得到結果，即　x= 0.0479, y= 2.1737, z= 3.6467。通常利用MATLAB解聯立方程式，亦會碰到無解的情形，此時或稱為特異狀況（或 |A| 為零的情形）。前述之左除結果，可用以反求Ｃ值，Ｃ＝Ａx，其結果也正確。

作業04_屬性資料查詢(SQL)

國立台北科技大學土木與防災研究所

地理資訊系統

作業04_屬性資料查詢(SQL)

學生：鄭佳邦學號：101429005

一、問題與討論：

1. 學習SQL查詢每平方公里人口密度大於500之國家時，利用RT SQL layer SQL查詢，指令如下：SELECT * FROM "public"."world" AS "t_0" WHERE "t_0"."Density" > '500' ，結果發生圖層無效之錯誤訊息(訊息如下)

答：經研判可能係 PostGIS support not enabled!，重新再次安裝設定PostGIS連線

二、作業成果

1.每平方公里人口密度大於500之國家計有24國，成果如下：

2.查詢全世界人口總合、面積總合及人口密度

利用DB Manager 查詢，指令如下：
select sum("pop_1994") as people,sum("area") as area,sum("pop_1994")/sum("area") as density from "public"."world"

結果如下：

3.查詢各洲大陸之國家數、人口數及面積：

利用DB Manager 查詢，指令如下：
select "continent",count(*) as country_no,sum("pop_1994") as people ,sum( "area") as area from "public"."world" group by continent

結果如下：