Re: [討論] 男生到底都多久打一次手槍呢?
0.序(引用)
泌尿科醫師 王起杰 提供一個小算式,只要把年齡的第一位數字乘上9
,一周該打幾次,一目瞭然。
快來算一算,這一套簡易公式→即年齡第一位數乘上9。
例如:30多歲者,3×9=27,為每2周7次;40多歲者,4×9=36,每3周6次;70多歲者,7
×9=63,每6周3次,依此類推。
----------------------------------------------------------------
1.前言:
這個公式當然只是個約略估計
它存在的弊病就是函數的不連續
假如有一個人他明天滿三十歲
那這個禮拜平均大約是1周8次
可是就在隔周他就瞬間變成大約1周3.5次
說實在如果年齡太接近這個臨界值那會很不準
所以現在我們來做一個更準的函數來估算
-----------------------------------------------------------------
2.估算考量因素
(1)因為醫師是說"二十幾歲"一周8次,所以我們取組中點也就是25歲
(2)只取"二十幾歲"、"三十幾歲"、"四十幾歲"、"五十幾歲"
六十歲以上的人就不納入考慮,這樣比較合情合理
(3)由上,可以取出(25,8)、(35,3.5)、(45,2)、(55,1.25)
-----------------------------------------------------------------
3.決定分析方法
(1)為方便計算,使用n次曲線做最小平方法之迴歸分析
(2)略畫曲線會發現十分不線性,故不採用一次曲線來做迴歸分析
因此,方便起見,用二次迴歸分析
-----------------------------------------------------------------
4.計算過程
(1)首先假設 ax^2 + bx + c = 0 是這條回歸二次曲線
(2)將題目寫成如下之表示方法:
625 25 1 a 8
令 A = [ 1225 35 1 ]、v = [ b ]、b = [ 3.5 ]
2025 45 1 c 2
3025 55 1 1.25
則題目變成:求Av和b這兩個四維度向量的最短距離
625 25 1
(3)由於Av屬於三向量[ 1225 ]、[ 35 ]、[1]張開的空間(即三向量的線性組合,記做S)
2025 45 1
3025 55 1
因此考慮的最短距離就是b對S的距離(廢話)
亦即b - (b在S上的投影)
(4)因此現在我們的目標是算出b在S上的投影
方便起見,先將b在S上的投影記為Proj(b)
則因為b-Proj(b)垂直於A的三個行向量
所以(A^T)(b-Proj(b)) = 0 ......(*) *註1:A^T為A的轉置矩陣,文末說明
又因為Proj(b)是S上的一個向量
所以存在一個三維向量u(即為我們要求的v)使得Proj(b) = Au
將以上結果代入(*)
得(A^T)(b-Au) = 0
化簡得
(A^T)b = (A^T)Au
所以u = [(A^T)A]^(-1)(A^T)b *註2:[(A^T)A]^(-1)為(A^T)A的反矩陣
*註3:[(A^T)A]有反矩陣
625 25 1 8
(5)所以將A = [1125 35 1 ]、b = [ 3.5 ]代入
2025 45 1 2
3025 55 1 1.25
0.01
得 u ≒ [ -1 ]
10.8
結論:
按照醫師給定之公式,年齡x歲的人會有每周大約0.01x^2-x+26.8的使用次數
註1:轉置矩陣就是和原矩陣行列互換的矩陣
如:
1 2 3
A = [ ]
4 5 6
1 4
則A^T = [ 2 5 ]
3 6
註2:反矩陣為和原矩陣相乘會等於單位矩陣(對角線為1,其餘為0)的方陣
如:
1 -1
B = [ ]
1 1
0.5 0.5
則B^(-1) = [ ]
-0.5 0.5
1 0
因為B*B^(-1) = B^(-1)*B = [ ]
0 1
註3(較為專業):
(1)因為A的三個行向量為線性獨立
因此rank(A) = 3,N(A) = 0
(2)假如存在向量u使得(A^T)Au = 0
則Au屬於N(A^T)也屬於R(A)
但是N(A^T) = V-R(A)(這裡的V就是vector space over R)
所以Au = 0,即u屬於N(A)
(3)但由(1)知道N(A)只包含0,所以N(A^T*A)也只有0
故det(A^T*A)≠0
即反矩陣存在
--
看完之後大家來做個練習吧
假如某星球中有這些人:
21歲乖學生是個不折不扣的處男,沒有和任何女生發生過性關係
35歲的某運動員雖然已婚,但是和315位女性發生婚外情,其中有三分之一為人所知
49歲的某夜店份子身經百戰,一共和438位女性發生婚外情
若某馬姓政客的年齡對婚外情對象數關係落在上述三人的年齡對婚外情對象數之迴歸直線
則現年63歲的他應該和多少個人發生過婚外情?
推文第一位答對的人可以獲得本篇的批幣
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.112.245.138
※ 文章網址: http://www.ptt.cc/bbs/sex/M.1398973560.A.19C.html
推
05/02 03:50, , 1F
05/02 03:50, 1F
推
05/02 03:53, , 2F
05/02 03:53, 2F
推
05/02 03:54, , 3F
05/02 03:54, 3F
推
05/02 04:09, , 4F
05/02 04:09, 4F
→
05/02 04:13, , 5F
05/02 04:13, 5F
正解!!可是我上面說要推文不好意思囉
推
05/02 04:15, , 6F
05/02 04:15, 6F
推
05/02 04:18, , 7F
05/02 04:18, 7F
錢是你的
噓
05/02 04:18, , 8F
05/02 04:18, 8F
噓
05/02 04:18, , 9F
05/02 04:18, 9F
噓
05/02 04:21, , 10F
05/02 04:21, 10F
推
05/02 04:22, , 11F
05/02 04:22, 11F
推
05/02 04:27, , 12F
05/02 04:27, 12F
推
05/02 04:43, , 13F
05/02 04:43, 13F
噓
05/02 04:47, , 14F
05/02 04:47, 14F
推
05/02 05:01, , 15F
05/02 05:01, 15F
推
05/02 05:15, , 16F
05/02 05:15, 16F
推
05/02 05:21, , 17F
05/02 05:21, 17F
噓
05/02 05:23, , 18F
05/02 05:23, 18F
不會很意外吧
畢竟只有二次
而且取點只有取到55
為了平衡過大的(25,8)這是可以預料的
噓
05/02 06:31, , 19F
05/02 06:31, 19F
噓
05/02 06:46, , 20F
05/02 06:46, 20F
噓
05/02 07:20, , 21F
05/02 07:20, 21F
推
05/02 07:33, , 22F
05/02 07:33, 22F
噓
05/02 07:49, , 23F
05/02 07:49, 23F
XDD已修
噓
05/02 08:24, , 24F
05/02 08:24, 24F
推
05/02 08:42, , 25F
05/02 08:42, 25F
推
05/02 08:49, , 26F
05/02 08:49, 26F
推
05/02 08:54, , 27F
05/02 08:54, 27F
推
05/02 09:19, , 28F
05/02 09:19, 28F
→
05/02 09:51, , 29F
05/02 09:51, 29F
推
05/02 09:56, , 30F
05/02 09:56, 30F
→
05/02 09:57, , 31F
05/02 09:57, 31F
我自己算出來是689啦
不過我等等用excel檢查一下好了
推
05/02 10:07, , 32F
05/02 10:07, 32F
推
05/02 10:19, , 33F
05/02 10:19, 33F
推
05/02 10:21, , 34F
05/02 10:21, 34F
噓
05/02 10:21, , 35F
05/02 10:21, 35F
推
05/02 10:22, , 36F
05/02 10:22, 36F
推
05/02 10:31, , 37F
05/02 10:31, 37F
推
05/02 10:41, , 38F
05/02 10:41, 38F
推
05/02 10:44, , 39F
05/02 10:44, 39F
→
05/02 10:47, , 40F
05/02 10:47, 40F
噓
05/02 10:48, , 41F
05/02 10:48, 41F
推
05/02 10:52, , 42F
05/02 10:52, 42F
推
05/02 10:56, , 43F
05/02 10:56, 43F
→
05/02 10:58, , 44F
05/02 10:58, 44F
推
05/02 11:17, , 45F
05/02 11:17, 45F
噓
05/02 11:22, , 46F
05/02 11:22, 46F
噓
05/02 11:50, , 47F
05/02 11:50, 47F
已改~~因為那是之前原po打錯的我沒有留意到
sorry
推
05/02 11:57, , 48F
05/02 11:57, 48F
推
05/02 12:09, , 49F
05/02 12:09, 49F
推
05/02 12:28, , 50F
05/02 12:28, 50F
※ 編輯: aibltjv (140.112.245.138), 05/02/2014 12:42:13
推
05/02 12:41, , 51F
05/02 12:41, 51F
推
05/02 12:51, , 52F
05/02 12:51, 52F
噓
05/02 12:55, , 53F
05/02 12:55, 53F
噓
05/02 13:27, , 54F
05/02 13:27, 54F
噓
05/02 15:06, , 55F
05/02 15:06, 55F
推
05/02 16:56, , 56F
05/02 16:56, 56F
推
05/02 17:38, , 57F
05/02 17:38, 57F
推
05/02 17:43, , 58F
05/02 17:43, 58F
推
05/02 19:08, , 59F
05/02 19:08, 59F
推
05/02 19:11, , 60F
05/02 19:11, 60F
→
05/02 19:11, , 61F
05/02 19:11, 61F
推
05/02 21:34, , 62F
05/02 21:34, 62F
推
05/02 21:35, , 63F
05/02 21:35, 63F
推
05/02 22:34, , 64F
05/02 22:34, 64F
噓
05/02 22:54, , 65F
05/02 22:54, 65F
推
05/02 23:19, , 66F
05/02 23:19, 66F
推
05/02 23:20, , 67F
05/02 23:20, 67F
推
05/03 00:39, , 68F
05/03 00:39, 68F
推
05/03 01:39, , 69F
05/03 01:39, 69F
推
05/03 01:40, , 70F
05/03 01:40, 70F
推
05/03 02:33, , 71F
05/03 02:33, 71F
推
05/03 04:17, , 72F
05/03 04:17, 72F
推
05/03 11:18, , 73F
05/03 11:18, 73F
推
05/03 11:48, , 74F
05/03 11:48, 74F
推
05/03 12:01, , 75F
05/03 12:01, 75F
推
05/03 13:20, , 76F
05/03 13:20, 76F
推
05/04 01:26, , 77F
05/04 01:26, 77F
推
05/06 16:22, , 78F
05/06 16:22, 78F
→
05/06 23:26, , 79F
05/06 23:26, 79F
推
05/08 23:49, , 80F
05/08 23:49, 80F
討論串 (同標題文章)
完整討論串 (本文為第 2 之 2 篇):
討論
35
87