pttchungyuandye: 12月 2010

2010年12月28日星期二

[自動轉寄] Re: 你好

作者: kai125 (阿花)
標題: Re: 你好
時間: Wed Dec 29 13:32:31 2010

我待會馬上寄mail給你

謝謝你

※ 引述《chungyuandye (養花種魚數月亮賞星星)》之銘言：
: ※ 引述《kai125 (阿花)》之銘言：
: : 不好意思我不會丟水球耶
: : 你有msn嗎?
: : 可以給我嗎
: : 謝謝
: 我沒有MSN說，直接Email給我吧
: chungyuandye at gmail
: 不過我不是念電機跟物理的
: 科科

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

[自動轉寄] (kai125) Re: 你好

作者: chungyuandye (養花種魚數月亮賞星星)
標題: Re: 你好
時間: Wed Dec 29 13:25:31 2010

※ 引述《kai125 (阿花)》之銘言：
: 不好意思我不會丟水球耶
: 你有msn嗎?
: 可以給我嗎
: 謝謝

我沒有MSN說，直接Email給我吧
chungyuandye at gmail

不過我不是念電機跟物理的
科科

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

http://cydye1069.blogspot.com
http://tinyurl.com/25nedr2

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

[自動轉寄] 你好

作者: kai125 (阿花)
標題: 你好
時間: Wed Dec 29 13:24:31 2010

不好意思我不會丟水球耶

你有msn嗎?

可以給我嗎

謝謝

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

2010年12月24日星期五

[自動轉寄] (manhoodguy) Re: [徵人] 協助資料分析

作者: chungyuandye (養花種魚數月亮賞星星)
標題: Re: [徵人] 協助資料分析
時間: Fri Dec 24 21:02:41 2010

※ 引述《manhoodguy (1111)》之銘言：
: ※ 引述《chungyuandye (養花種魚數月亮賞星星)》之銘言：
: : 可以吧資料給我試試嗎？？
: : 可以的話，請聯絡我！
: : 感恩∼
: 我可以先請您試分析一個簡單問題嗎?
: 若您有興趣可以給我站外信箱?
chungyuandye@gmail.com

--
我打研究室走過那獨坐電腦前的容顏如苦瓜的糾結
靈感不來長壽的煙霧不散
研究室如小小的寂寞的城恰如商管的電梯向晚

http://chungyuandye.twbbs.org

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

[自動轉寄] Re: [徵人] 協助資料分析

作者: manhoodguy (1111)
標題: Re: [徵人] 協助資料分析
時間: Fri Dec 24 08:13:17 2010

※ 引述《chungyuandye (養花種魚數月亮賞星星)》之銘言：
: ※ 引述《manhoodguy (1111)》之銘言：
: : 徵求單位:個人
: : 工作內容:因素或主成份資料分析
: : 工作時間:皆可
: : 工作薪資:約2000~3000
: : 工作地點:電子信件
: : PS:本版不鼓勵論文統計程式全部代跑的行為
: : 但是指導使用何種工具以及跑法則再接受範圍
: 可以吧資料給我試試嗎？？
: 可以的話，請聯絡我！
: 感恩∼

我可以先請您試分析一個簡單問題嗎?
若您有興趣可以給我站外信箱?

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

2010年12月20日星期一

[自動轉寄] (peter5214) Re: [其他] mathematica 求方程式的解

作者: chungyuandye (養花種魚數月亮賞星星)
標題: Re: [其他] mathematica 求方程式的解
時間: Mon Dec 20 16:13:20 2010

※ 引述《peter5214 (阿宏)》之銘言：
: 最近在用mathematica求equation的解，
: 遇到一個很怪的事情，常常會遇到他解出來的解
: 會附加一個 0.i
: EX: 0.517293+0.i
: 既然虛部是零又何必再列出來呢?? 還是說這有什麼涵義在??
: 因為有時後解出來的解又不會跑出 0.i這樣的東西，
: 還有看到0.-1.29896i 這樣的東西，我查了一下實部是0虛部是-1.29896
: 但一般不是就直接寫成a+bi的形式嗎? 跑出這樣的東西是不是有哪出問題了??
: 有請高手代為回答一下了 (_ _) 感謝

Solve[x^3 == 1, x]
N@Solve[x^3 == 1, x]
你會發現Mathematica並沒有出包，只是他把x^3=-1的所有的根中輸出！
如果應要忽略複數，那解除Power這個函數的封印
Unprotect[Power];
Power[x_?Negative, Rational[p_, q_?OddQ]] := (-(-x)^(1/q))^p;
Protect[Power];
Solve[x^3 == -1, x]
不過這個作法不一定好就是了∼∼

詳細內容可以參考
http://cydye1069.blogspot.com/2009/11/113-bug.html

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

http://cydye1069.blogspot.com
http://tinyurl.com/25nedr2

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

2010年12月8日星期三

[自動轉寄] (af2486) Re: 您好，想請問關於mathematica的問題

作者: chungyuandye (養花種魚數月亮賞星星)
標題: Re: 您好，想請問關於mathematica的問題
時間: Wed Dec 8 21:33:29 2010

Do那邊改成這樣看看可不可以

Quiet@FindRoot[f[(y (2 \[Pi]))/#, opt[#]] == 0, {y, 1.44677}] & /@
Range[1.54, 1.55, 0.005]

0.005是是公差，如果要做大量運算的話，自己找一下ParallelMap...

※ 引述《af2486 (我喜歡旅行)》之銘言：
: 您好，之前在暑假的時候曾經請教過mathematica的問題，非常感激
: 不過現在又有些問題想請教，不知是否有空撥冗幫我看一下呢?謝謝!
: 這是您之前幫我改過的程式碼：
: xiLimit = 15.;
: Hankel1[m_, x_] :=
: If[Abs[x] < xiLimit, BesselJ[m, x] + I BesselY[m, x],
: Sqrt[2/(Pi*x)] Exp[I*(x - m/2*Pi - Pi/4)]];
: BesselJp[m_, x_] :=
: If[Abs[x] < xiLimit,
: D[BesselJ[m, x1],
: x1] /. {x1 -> x}, -Sqrt[2/(Pi*x)] Sin[(x - m/2*Pi - Pi/4)]];
: BesselYp[m_, x_] :=
: If[Abs[x] < xiLimit, D[BesselY[m, x1], x1] /. {x1 -> x},
: Sqrt[2/(Pi*x)] Cos[(x - m/2*Pi - Pi/4)]];
: Hankel1p[m_, x_] :=
: If[Abs[x] < xiLimit, D[Hankel1[m, x1], x1] /. {x1 -> x},
: I*Sqrt[2/(Pi*x)] Exp[I*(x - m/2*Pi - Pi/4)]];
: kr[kz_, k_] := Sqrt[k^2 - kz^2];
: \[Alpha]r[kz_, k_] := Sqrt[kz^2 - k^2];
: f[kz_, opt_] := Module[{sub},
: k0 = \[Omega] Sqrt[\[Epsilon]0 \[Mu]0] /. opt;
: k1 = \[Omega] Sqrt[\[Epsilon]1 \[Mu]1] /. opt;
: k2 = \[Omega] Sqrt[\[Epsilon]2 \[Mu]2] /. opt;
: r0 = r0 /. opt;
: r1 = r1 /. opt;
: m = m /. opt;
: kr0 = kr[kz, k0];
: kr1 = kr[kz, k1];
: \[Alpha]r2 = \[Alpha]r[kz, k2];
: BJ00 = BesselJ[m, kr[kz, k0] r0];
: BJ10 = BesselJ[m, kr[kz, k1] r0];
: BJ11 = BesselJ[m, kr[kz, k1] r1];
: BY10 = BesselY[m, kr[kz, k1] r0];
: BY11 = BesselY[m, kr[kz, k1] r1];
: BJp00 = BesselJp[m, kr[kz, k0] r0];
: BJp10 = BesselJp[m, kr[kz, k1] r0];
: BJp11 = BesselJp[m, kr[kz, k1] r1];
: BYp10 = BesselYp[m, kr[kz, k1] r0];
: BYp11 = BesselYp[m, kr[kz, k1] r1];
: HA21 = Hankel1[m, I \[Alpha]r[kz, k2] r1];
: HAp21 = Hankel1p[m, I \[Alpha]r[kz, k2] r1];
: M = {{BJ00, 0, -BJ10, -BY10 , 0, 0, 0, 0},
: {0, BJ00, 0, 0, -BJ10, -BY10, 0, 0},
: {(\[Omega] \[Epsilon]0)/(kr0 r0) BJp00, (m kz)/(kr0^2 r0^2)
: BJ00, -((\[Omega] \[Epsilon]1)/(kr1 r0))
: BJp10, -((\[Omega] \[Epsilon]1)/(kr1 r0)) BYp10, -((m kz)/(
: kr1^2 r0^2)) BJ10, -((m kz)/(kr1^2 r0^2)) BY10, 0, 0},
: {(m kz)/(kr0^2 r0^2) BJ00, (\[Omega] \[Mu]0)/(kr0 r0)
: BJp00, -((m kz)/(kr1^2 r0^2)) BJ10, -((m kz)/(kr1^2 r0^2))
: BY10, -((\[Omega] \[Mu]1)/(kr1 r0))
: BJp10, -((\[Omega] \[Mu]1)/(kr1 r0)) BYp10, 0, 0 },
: {0, 0, BJ11, BY11 , 0, 0, -HA21, 0},
: {0, 0, 0, 0, BJ11, BY11, 0, -HA21},
: {0, 0, (\[Omega] \[Epsilon]1)/(kr1 r1)
: BJp11, (\[Omega] \[Epsilon]1)/(kr1 r1) BYp11, (m kz)/(
: kr1^2 r1^2) BJ11, (m kz)/(kr1^2 r1^2) BY11, (
: I \[Omega] \[Epsilon]2)/(\[Alpha]r2 r1) HAp21, (
: m kz)/(\[Alpha]r2^2 r1^2) HA21},
: {0, 0, (m kz)/(kr1^2 r1^2) BJ11, (m kz)/(kr1^2 r1^2)
: BY11, (\[Omega] \[Mu]1)/(kr1 r1) BJp11, (\[Omega] \[Mu]1)/(
: kr1 r1) BYp11, (m kz)/(\[Alpha]r2^2 r1^2) HA21, (
: I \[Omega] \[Mu]2)/(\[Alpha]r2 r1) HAp21}
: } /. opt;
: equ = Det[M];
: Return[equ];
: ]
: n0[lambda_] :=
: Sqrt[1 + (0.68671749*lambda^2)/(lambda^2 -
: 0.072675189^2) + (0.43481505*lambda^2)/(lambda^2 -
: 0.11514351^2) + (0.89656582*lambda^2)/(lambda^2 - 10.002398^2)];
: n1[lambda_] :=
: Sqrt[1 + (0.6961663*lambda^2)/(lambda^2 - 0.0684043^2) + (0.4079426*
: lambda^2)/(lambda^2 - 0.1162414^2) + (0.8974794*
: lambda^2)/(lambda^2 - 9.896161^2)];
: n2[lambda_] := 1.0;
: opt[lambda_] := {m -> 1, \[Omega] -> 2 Pi/lambda, r0 -> 4.1,
: r1 -> 62.5, \[Epsilon]0 -> n0[lambda]^2, \[Mu]0 ->
: 1., \[Epsilon]1 -> n1[lambda]^2, \[Mu]1 -> 1., \[Epsilon]2 ->
: n2[lambda]^2, \[Mu]2 -> 1., kz -> y};
: temp3 = NestList[{#[[1]] + 0.001, #[[3]], Re@y} /.
: FindRoot[
: f[y*(2 Pi/(#[[1]] + 0.001)), opt[#[[1]] + 0.001]] ==
: 0, {y, #[[3]]}] &, {0.999, 1.299, 1.45038}, 301][[2 ;; -1]]
: Export["temp3.csv", temp3[[All, 3]]]
: -------------------------------------------------------------
: 其中迴圈是用牛頓法在1~1.3之間每隔0.001找一次解
: 理想上想得到的解是約略成等差級數的解
: 可是現在遇到的問題是在某些點，找到的解會和上一點一樣，這並不是所想要的
: 請問這有可能程式哪邊出問題嗎?還是跟牛頓法有關呢?
: 不好意思麻煩了...謝謝!

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

http://cydye1069.blogspot.com
http://tinyurl.com/25nedr2

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

[自動轉寄] 您好，想請問關於mathematica的問題

作者: af2486 (我喜歡旅行)
標題: 您好，想請問關於mathematica的問題
時間: Wed Dec 8 17:39:28 2010

您好，之前在暑假的時候曾經請教過mathematica的問題，非常感激

不過現在又有些問題想請教，不知是否有空撥冗幫我看一下呢?謝謝!

這是您之前幫我改過的程式碼：

xiLimit = 15.;
Hankel1[m_, x_] :=
If[Abs[x] < xiLimit, BesselJ[m, x] + I BesselY[m, x],
Sqrt[2/(Pi*x)] Exp[I*(x - m/2*Pi - Pi/4)]];
BesselJp[m_, x_] :=
If[Abs[x] < xiLimit,
D[BesselJ[m, x1],
x1] /. {x1 -> x}, -Sqrt[2/(Pi*x)] Sin[(x - m/2*Pi - Pi/4)]];
BesselYp[m_, x_] :=
If[Abs[x] < xiLimit, D[BesselY[m, x1], x1] /. {x1 -> x},
Sqrt[2/(Pi*x)] Cos[(x - m/2*Pi - Pi/4)]];
Hankel1p[m_, x_] :=
If[Abs[x] < xiLimit, D[Hankel1[m, x1], x1] /. {x1 -> x},
I*Sqrt[2/(Pi*x)] Exp[I*(x - m/2*Pi - Pi/4)]];
kr[kz_, k_] := Sqrt[k^2 - kz^2];
\[Alpha]r[kz_, k_] := Sqrt[kz^2 - k^2];

f[kz_, opt_] := Module[{sub},
k0 = \[Omega] Sqrt[\[Epsilon]0 \[Mu]0] /. opt;
k1 = \[Omega] Sqrt[\[Epsilon]1 \[Mu]1] /. opt;
k2 = \[Omega] Sqrt[\[Epsilon]2 \[Mu]2] /. opt;
r0 = r0 /. opt;
r1 = r1 /. opt;
m = m /. opt;
kr0 = kr[kz, k0];
kr1 = kr[kz, k1];
\[Alpha]r2 = \[Alpha]r[kz, k2];
BJ00 = BesselJ[m, kr[kz, k0] r0];
BJ10 = BesselJ[m, kr[kz, k1] r0];
BJ11 = BesselJ[m, kr[kz, k1] r1];
BY10 = BesselY[m, kr[kz, k1] r0];
BY11 = BesselY[m, kr[kz, k1] r1];
BJp00 = BesselJp[m, kr[kz, k0] r0];
BJp10 = BesselJp[m, kr[kz, k1] r0];
BJp11 = BesselJp[m, kr[kz, k1] r1];
BYp10 = BesselYp[m, kr[kz, k1] r0];
BYp11 = BesselYp[m, kr[kz, k1] r1];
HA21 = Hankel1[m, I \[Alpha]r[kz, k2] r1];
HAp21 = Hankel1p[m, I \[Alpha]r[kz, k2] r1];
M = {{BJ00, 0, -BJ10, -BY10 , 0, 0, 0, 0},
{0, BJ00, 0, 0, -BJ10, -BY10, 0, 0},
{(\[Omega] \[Epsilon]0)/(kr0 r0) BJp00, (m kz)/(kr0^2 r0^2)
BJ00, -((\[Omega] \[Epsilon]1)/(kr1 r0))
BJp10, -((\[Omega] \[Epsilon]1)/(kr1 r0)) BYp10, -((m kz)/(
kr1^2 r0^2)) BJ10, -((m kz)/(kr1^2 r0^2)) BY10, 0, 0},
{(m kz)/(kr0^2 r0^2) BJ00, (\[Omega] \[Mu]0)/(kr0 r0)
BJp00, -((m kz)/(kr1^2 r0^2)) BJ10, -((m kz)/(kr1^2 r0^2))
BY10, -((\[Omega] \[Mu]1)/(kr1 r0))
BJp10, -((\[Omega] \[Mu]1)/(kr1 r0)) BYp10, 0, 0 },
{0, 0, BJ11, BY11 , 0, 0, -HA21, 0},
{0, 0, 0, 0, BJ11, BY11, 0, -HA21},
{0, 0, (\[Omega] \[Epsilon]1)/(kr1 r1)
BJp11, (\[Omega] \[Epsilon]1)/(kr1 r1) BYp11, (m kz)/(
kr1^2 r1^2) BJ11, (m kz)/(kr1^2 r1^2) BY11, (
I \[Omega] \[Epsilon]2)/(\[Alpha]r2 r1) HAp21, (
m kz)/(\[Alpha]r2^2 r1^2) HA21},
{0, 0, (m kz)/(kr1^2 r1^2) BJ11, (m kz)/(kr1^2 r1^2)
BY11, (\[Omega] \[Mu]1)/(kr1 r1) BJp11, (\[Omega] \[Mu]1)/(
kr1 r1) BYp11, (m kz)/(\[Alpha]r2^2 r1^2) HA21, (
I \[Omega] \[Mu]2)/(\[Alpha]r2 r1) HAp21}
} /. opt;
equ = Det[M];
Return[equ];
]

n0[lambda_] :=
Sqrt[1 + (0.68671749*lambda^2)/(lambda^2 -
0.072675189^2) + (0.43481505*lambda^2)/(lambda^2 -
0.11514351^2) + (0.89656582*lambda^2)/(lambda^2 - 10.002398^2)];

n1[lambda_] :=
Sqrt[1 + (0.6961663*lambda^2)/(lambda^2 - 0.0684043^2) + (0.4079426*
lambda^2)/(lambda^2 - 0.1162414^2) + (0.8974794*
lambda^2)/(lambda^2 - 9.896161^2)];

n2[lambda_] := 1.0;
opt[lambda_] := {m -> 1, \[Omega] -> 2 Pi/lambda, r0 -> 4.1,
r1 -> 62.5, \[Epsilon]0 -> n0[lambda]^2, \[Mu]0 ->
1., \[Epsilon]1 -> n1[lambda]^2, \[Mu]1 -> 1., \[Epsilon]2 ->
n2[lambda]^2, \[Mu]2 -> 1., kz -> y};

temp3 = NestList[{#[[1]] + 0.001, #[[3]], Re@y} /.
FindRoot[
f[y*(2 Pi/(#[[1]] + 0.001)), opt[#[[1]] + 0.001]] ==
0, {y, #[[3]]}] &, {0.999, 1.299, 1.45038}, 301][[2 ;; -1]]

Export["temp3.csv", temp3[[All, 3]]]

-------------------------------------------------------------
其中迴圈是用牛頓法在1~1.3之間每隔0.001找一次解

理想上想得到的解是約略成等差級數的解

可是現在遇到的問題是在某些點，找到的解會和上一點一樣，這並不是所想要的

請問這有可能程式哪邊出問題嗎?還是跟牛頓法有關呢?

不好意思麻煩了...謝謝!

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

訂閱：文章 (Atom)

2010年12月28日 星期二

2010年12月24日 星期五

2010年12月20日 星期一

2010年12月8日 星期三

2010年12月28日星期二

2010年12月24日星期五

2010年12月20日星期一

2010年12月8日星期三