2012年12月11日 星期二

Re: [問題] 在XY平面上的蛇行點找所有的最低值!?

作者: chungyuandye (養花種魚數月亮賞星星) 看板: Mathematica
標題: Re: [問題] 在XY平面上的蛇行點找所有的最低值!?
時間: Tue Nov 6 00:35:04 2012

※ 引述《candy88257 (阿泰斯)》之銘言:
: 標題: [問題] 在XY平面上的蛇行點找所有的最低值!?
: 時間: Mon Nov 5 22:36:00 2012
:
: 如:
:
: 給定一連串X座標為連續的座標點:(1,1),(2,3),(3,5),(4,2),(5,3),(6,7),(7,2),(8,5),(9,6)
: 畫起來像是上下起伏的震動圖
: 請問如何很快速的在裡面找到所有的最低點?
: 如上面串列的所有最低點為: (4,2),(7,2)
: PS.
: 不要用畫圖的去判斷
: 因為我目前要判斷的點都是好幾百個一組...
:
: --
: ※ 發信站: 批踢踢實業坊(ptt.cc)
: ◆ From: 175.181.151.43
: 推 jurian0101:If用法有誤,你要多看範例文件熟悉MM的語法了。 11/05 22:44
: → candy88257:意思是可以用If寫出來!? 11/05 23:00
: → candy88257:感謝指導! 11/05 23:00
: → chungyuandye:SplitBy[SortBy[data,Last],Last] 11/05 23:04
: → chungyuandye:SplitBy[SortBy[data, Last], Last][[1]] 少打[[1]] 11/05 23:04

如果要找全域極小,又不考慮第一組資料,那下列語法就可以
SplitBy[SortBy[data[[2;;-1]], Last], Last][[1]]

: → candy88257:我想找到所有的最低點...就像上面例子最後抓到2個最低 11/05 23:31
: → candy88257:也就是所有點按X順序在XY平面連起來,找到所有的最低點 11/05 23:32
: → candy88257:其實不能講最低點,應該是相對低點 11/05 23:37

如果要找相對高低點,那Split這個函數可發揮妙用

(* 產生100組,每組200個數字的隨機資料 *)
temp = Sort[#] & /@ RandomReal[{0, 50 Pi}, {100, 200}];
data = Map[{#, Sin[#]} &, temp, {2}]

(* 畫出所有資料的序列圖 *)
Graphics[Line[#] & /@ Partition[#, 2, 1], Epilog -> {Red, Point[#]},
Frame -> True, ImageSize -> 1000] & /@ data

(* 第一組資料的圖形 *)
Graphics[Line[#] & /@ Partition[data[[1]], 2, 1],
Epilog -> {Red, Point[data[[1]]]}, Frame -> True, ImageSize -> 1000]


(* 先分析第一組,找出所有的相對低點。若該組資料的y值若上一組資料y值小, *)
(* 就歸為同一組。 所以若長度為一則表示資料在遞增,故剔除。 *)

low = Select[Split[data[[1]], #1[[2]] > #2[[2]] &],
Length[#] > 1 &][[All, -1]]

(* 同樣的方法可應用在相對高點 *)
high = Select[Split[data[[1]], #1[[2]] <= #2[[2]] &],
Length[#] > 1 &][[All, -1]]

(* 輸出所有資料的相對高低點 *)
datalow =
Select[Split[data[[#]], #1[[2]] > #2[[2]] &],
Length[#] > 1 &][[All, -1]] & /@ Range[Length@data]

datahigh =
Select[Split[data[[#]], #1[[2]] <= #2[[2]] &],
Length[#] > 1 &][[All, -1]] & /@ Range[Length@data]

--
養花種魚數月亮賞星星

http://chungyuandye.twbbs.org


--
※ 發信站: 批踢踢實業坊(ptt.cc)
※ 編輯: chungyuandye 來自: 218.173.133.152 (11/06 00:48)
推 candy88257:我再研究一下,感謝!!  11/07 16:36

Re: [問題] mathmatica問題

作者: chungyuandye (養花種魚數月亮賞星星) 看板: Mathematica
標題: Re: [問題] mathmatica問題
時間: Thu Nov 29 22:48:09 2012

※ 引述《abbybao (小寶)》之銘言:
: 1. http://ppt.cc/NGy7
: 想問一下我想要畫這函數但為什麼x<0的部分不會顯示呢?
: 我程式有寫錯嗎@@

n的m次方,若m為分數且n為負數時,Mathematica會設定為複數。
因為是複數,左半邊才會沒畫出來。

Power[-1, 1/3]
Power[1,1/3]

(* (-1)^(1/3),複數 *)
Power[-1, 1/3] // N

ContourPlot[1-x^(2/3)==y,{x,-5,5},{y,-5,5}]

(* 方法1,自行做調整 *)
ContourPlot[(1-y)^3==x^2,{x,-5,5},{y,-5,5}]

(* 方法2,如果你根本不想要複數,就修改Power這個函數 *)
Unprotect[Power];
Power[x_?Negative,Rational[p_,q_?OddQ]]:=(-(-x)^(1/q))^p;
Protect[Power];

Power[-1, 1/3] // N
ContourPlot[1-x^(2/3)==y,{x,-5,5},{y,-5,5}]

: 2. mathmatica有可以找反函數的功能嗎?我打關鍵字inverse function
: 好像只能
: in: InverseFunction[f][x-7]
: out: f(-1)[x-7]
: 但我希望它是能顯示真值而不是以含數形式表示
: ex: f(x)=x-7 f(-1)(x)=x+7 <----希望顯示的值
: 不知道能不能做的到
: 拜託大家了,謝謝^^

f[x_] = x - 7;
InverseFunction[f][x]


--
養花種魚數月亮賞星星

http://chungyuandye.twbbs.org


--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.173.135.172
推 abbybao:感謝回應,不過有點疑問的地方,就是我算求雙曲線sinhx的  11/30 15:35
→ abbybao:的反函數時,差了一個正負號,不知道是甚麼原因?如圖  11/30 15:36
→ abbybao:http://tinyurl.com/cjepw9e  11/30 15:37
推 jurian0101:第九版多了一個Surd[R,n]函數,求實數n次方實根  11/30 18:04
→ abbybao:可惜還沒有第九版,過一陣子論壇應該就有破解板的了  12/02 12:39

Re: [問題] mathematica可以做以下運算嗎?

作者: chungyuandye (養花種魚數月亮賞星星) 看板: Mathematica
標題: Re: [問題] mathematica可以做以下運算嗎?
時間: Fri Nov 30 23:40:31 2012

※ 引述《abbybao (小寶)》之銘言:
: 1.http://tinyurl.com/cqp98qj
: 算極座標的面積

r=4-Cos[x];
2*Integrate[r^2*Pi/(2Pi),{x,0,Pi}]

: 2.http://tinyurl.com/cx9fuut
: 算卡氏座標的弧長

y=Log[1-x^2];
temp=Refine[
Simplify[Sqrt[Dt[x]^2+Dt[y]^2]]/Dt[x],{Dt[x]>0,0<x<1/4}]
ans1=Integrate[temp,{x,0,1/4}]
ans2=TrigToExp@ans1//Simplify


: 3.http://tinyurl.com/cfrosjw
: 算極座標的弧長
: 弧長部分我只有找到arclengthfactor這個指令,但是他範例只有參數座標的
: 直角跟極座標不知道怎麼寫?

r=1+Cos[x];
f[x_]={r*Cos[x],r*Sin[x]}
D[f[x],x]
Integrate[Sqrt[D[f[x],x].D[f[x],x]],{x,0,2Pi}]


: 4.http://tinyurl.com/c2mb8pd
: 算旋轉體的體積

Limit[Sum[Pi*(4/((2/n*i)^2+4))^2*2/n,{i,1,n}],n->Infinity]

Manipulate[
Graphics3D[
Cylinder[{{2/n*(#-1),0,0},{2/n*#,0,0}},
4/((2/n*#)^2+4)]&/@Range[n]],
{n,5,200,5}
]

: 5. 1/A = 1/9955 + 1/9956 + ...... + 1/9960
: 求 A 之值 ? (ANS:約1647)
: 這題我用程式寫出來卻有兩個解 π 跟 1/π
: 不知哪邊寫錯了?
: http://ppt.cc/4P7B
: 拜託大家了,謝謝^^


--
養花種魚數月亮賞星星

http://chungyuandye.twbbs.org


--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.173.135.172
→ abbybao:感恩,我弄出來了^^  12/01 12:55

Re: [問題] 矩陣行列變換

作者: chungyuandye (養花種魚數月亮賞星星) 看板: Mathematica
標題: Re: [問題] 矩陣行列變換
時間: Sat Nov 17 17:32:54 2012

※ 引述《xnicky2000 (愛談低調的地瓜)》之銘言:
: 我現在有個矩陣
: ┌ ┐
: │ a b c d │
: │ │
: │ │
: │ e f g h │
: │ │
: │ │
: │ i j k l │
: │ │
: │ │
: │ m n o p │
: │ │
: └ ┘
: 如果我要把 m n o p
: 往上移變成
: ┌ ┐
: │ a b c d │
: │ │
: │ m n o p │
: │ │
: │ e f g h │
: │ │
: │ i j k l │
: │ │
: └ ┘
: 要怎樣做啊?? 可以用TABLE做嗎?
: 因為如果之後可能又有別行要換 可以比較好了解
: 感謝
: 版本是7.0的

data = {{a, b, c, d}, {e, f, g, h}, {i, j, k, l}, {m, m, 0, p}};

Insert[RotateLeft[data[[2 ;; -1]], -1], data[[1]], 1]

--
養花種魚數月亮賞星星

http://chungyuandye.twbbs.org


--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.173.137.59

Re: [問題] mathematica可以做以下運算嗎?

作者: chungyuandye (養花種魚數月亮賞星星) 看板: Mathematica
標題: Re: [問題] mathematica可以做以下運算嗎?
時間: Wed Nov 21 13:29:49 2012

※ 引述《abbybao (小寶)》之銘言:
: 1.在座標中給定三個點形成一個三角形
: 要怎麼求在這個三角形內包含的所有座標?
: ex: (0,0) (3,0) (3,3)包含了(2,1) 一點(邊界的點不算)

test[x_,y_]:=Boole[{y>0&&-3+x<0&&x-y>0}][[1]]
Select[Flatten[Table[{x,y},{x,0,3},{y,0,3}],1],test@#==1&]


: 2.http://ppt.cc/nD~K
: 求兩切線

ContourPlot[2x^3+2y^3==9x*y,{x,0,10},{y,0,10},
Epilog->{Red,PointSize[0.1],Point[{6,6}]}]

m[x_,y_]=Dt[y,x]/.Solve[Dt[2x^3+2y^3==9x*y,x],Dt[y,x]][[1]]
ans={x,y}/.Solve[{2x^3+2y^3==9x*y,y-6==m[x,y](x-6)},{x,y}]

y-#[[2]]==m@@#*(x-#[[1]])&/@ans[[1;;2]]

ContourPlot[2x^3+2y^3==9x*y,{x,0,10},{y,0,10},
Epilog->{Line[{ans[[1]],{6,6}}],Line[{ans[[2]],{6,6}}],
Red,PointSize[0.1],Point[{6,6}]}]

: 3.http://ppt.cc/ydlH
: 求面積
: 我的意思是有沒有指令是可以把邊界條件打出來它會幫你算面積
: ex: 此題邊界條件為 y= 2x-x^2 , y= x-2

Integrate[Abs[2x-x^2-(x-2)],Flatten@{x,x/.Solve[2x-x^2==(x-2),x]}]

: 4.列出100的所有質數

Select[{#,PrimeQ[#]}&/@Range[100],#[[2]]==True&][[All,1]]

: 拜託大家了,謝謝^^

http://chungyuandye.twbbs.org/2010/07/mathematica_11.html
http://chungyuandye.twbbs.org/2010/07/mathematica_18.html
--
我打研究室走過 那獨坐電腦前的容顏如苦瓜的糾結
靈感不來 長壽的煙霧不散
研究室如小小的寂寞的城 恰如商管的電梯向晚

http://chungyuandye.twbbs.org

--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 118.233.128.85
※ 編輯: chungyuandye 來自: 118.233.128.85 (11/21 13:32)
推 abbybao:哇靠,太強了,真的有辦法做出來,不過第一題我執行後是  11/21 16:49
→ abbybao:空集合,不知道是哪邊出問題 http://ppt.cc/SKtY  11/21 16:50
→ abbybao:應該想說會有一點(2,1)出來  11/21 16:50
→ chungyuandye:改成test@@#  11/21 21:05

Re: [問題] FindRoot 找不到後的後續工作

作者: chungyuandye (養花種魚數月亮賞星星) 看板: Mathematica
標題: Re: [問題] FindRoot 找不到後的後續工作
時間: Thu Nov 29 23:20:21 2012

※ 引述《ntust661 (TOEFL_5!)》之銘言:
: 我想說FindRoot動作失敗後,想要取而代之改變動作
: 但是我要如何打出當FindRoot失敗後,我改變我的動作
: EX:
: A = 1 ;
: s = FindRoot[ x^2 + x + A , {x , 1} ] ;
: 以下如果找根失敗,我就要停止動作接下來改 A = 2 , 3, 4 ...
: 我要怎麼做這個動作@@
: 我看 F1 他說找根失敗,他依舊會找出奇怪的根,我就不知道判斷的依據在哪
: 有請版友幫忙了ˊˋ

f[x_,A_]:=x^2+x+A;

NestWhile[(temp=x/.Quiet@FindRoot[f[x,#],{x,1}];
Print@{#,Abs[f[temp,#]]};
If[Abs[f[temp,#]]<10^-5,#,#-1])&,10,UnsameQ,2]

不過有警告訊息不表示說是失敗,大部分都是精確度或是newton method的條件
如果不想要這些睛告訊息,FindRoot改成Quiet@FindRoot。

Mathematica的FindRoot這個指令預設是Newton method,所以要求跟的話
最好的方式就是先Plot,找到一個好個initial。




--
養花種魚數月亮賞星星

http://chungyuandye.twbbs.org


--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.173.135.172

Re: [問題] mathematica旋轉體問題

作者: chungyuandye (養花種魚數月亮賞星星) 看板: Mathematica
標題: Re: [問題] mathematica旋轉體問題
時間: Sun Dec 2 09:47:52 2012

※ 引述《abbybao (小寶)》之銘言:
: 1. http://tinyurl.com/czwfb6v
: Q1:如果我想對 y=sin(x)對x=π/2 為旋轉軸旋轉要怎麼寫呢?
: ps: 圖裡是對z軸旋轉
: Q2:如果我想要標註點(1,0,0)的位置這有辦法做的到嗎?
: ps: 圖裡是點是我自己亂標的

(* RevolutionPlot3D 的轉軸必須是有原點出發的線,可參考 RevolutionAxis設定 *)
(* 所以要對x=π/2旋轉,必須把圖形平移,之後再做調整。 *)

temp = RevolutionPlot3D[Sin[Pi/2 - x], {x, -Pi/2, 0},
BoxRatios -> {1, 1, 0.4}]

Show[ListPlot3D[temp[[1, 1]] /. {x_, y_, z_} -> {x + Pi/2, y, z},
PlotStyle -> Opacity[0.5]],
Graphics3D[{PointSize[0.025], Red, Point[{1, 0, 0}]}]]


: 2. http://ppt.cc/LO-C
: 我想要看對題目裡的拋物面跟圓柱兩者圖形取交集所產生的圖形
: 這有辦法做的到嗎?

Plot3D[(x^2 + y^2)/4, {x, -10, 10}, {y, -10, 10},
RegionFunction -> Function[{x, y}, x^2 + y^2 <= 8 y],
Filling -> Bottom]

--
養花種魚數月亮賞星星

http://chungyuandye.twbbs.org


--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.173.134.123
→ abbybao:太屌了,真的弄得出來,尤其是第二題 這對我滿有幫助的  12/02 12:40
→ abbybao:第一題我執行程式後還是沒有看到點在圖型上顯示,如圖  12/02 12:42
→ abbybao:http://ppt.cc/SrnK 我把pointsize調到30還是沒看到  12/02 12:43
→ chungyuandye:30太大了,把整個點都包住了。  12/02 15:17