Please tell me what you thought about this blog. Comments, suggestions and views are welcomed.
這是留言板-既然來了，就留個言再走吧！
由 養花種魚數月亮賞星星 於 2029年7月22日星期日
下午2:43 發表
Please tell me what you thought about this blog. Comments, suggestions and views are welcomed.
讀者回應 ( 66 意見 )
訂閱發佈留言 (Atom)
請問一下，我用data= {0.02449541456596422`,0.0525439559116545`,1.4439337620566677`,1.9340496407236472`, 2.515622985353725`,3.0472338607735465`,3.7955493339996313`,4.632960983675772`,6.589281610279324`}
去畫圖，ListPlot[data, PlotJoined -> True, PlotStyle -> {Red}]，問題是我想要改坐標軸的刻度，要怎麼改啊，例如：全部改成兩倍?
ListPlot[data, PlotJoined -> True,
PlotStyle -> {Red},
Ticks -> {Range[10], Range[0, 7, 0.5]}
]
Ticks ->{x的刻度陣列,y的刻度陣列}
抱歉，我從PPT連到您的留言板，想向您請教有關EXCEL或VBA的問題
因為內容比較多
我把資料放在
http://docs.google.com/Doc?id=dxp8v42_0hh82pdt9&hl=en
拜託您幫我看看好嗎?感謝
同學來給你簽個名...
andychen
呵呵～～～真罕見
五十嵐的大老板
您好
我也寫了自己的Mathematica教學，到這裡來野人獻曝一番
http://wxlin1983.blogspot.com/2009/06/mathematica.html
W. X. Lin
老師～～我是睿昇
看來你要準備當拔比嘍
系阿～～真是感恩喔！
您好，我是在網路上看到您的關於SPSS的教學而連結過來的。
因為我正在寫論文。卻又從沒學過SPSS和統計。看中文書又絕大部分處在有看沒有懂的狀態。又苦於論文截止日近在眼前。
所以可否請您指點一下我的假設可以跑什麼分析方法。
非常謝謝您。
我是要探討廣告，試用品促銷，口碑vs消費者行為。
我的假設如下：
1.廣告對消費者行為有明顯正向關係
2.試用品的促銷手法對消費者行為有明顯正向關係
3.口碑對消費者行為有明顯正向關係
4.以上三者中，試用品促銷手法影響消費者行為的效果大於廣告，口碑影響消費者行為的效果大於試用品促銷。三者中，口碑影響消費者的效果最顯著。
另外，
我會跑一個信度分析。
但是我不確定我要不要跑效度。因為我的問卷裡當初在第一題時就忘了問：請問你是否瀏覽過此一網站。來判斷問卷是否有效。
所以我也想問一下這樣還要跑效度嗎？
非常非常感謝您。
您好，
請問在spss輸入資料時。
因為問卷其中一題題目是設定有經驗的人答題,沒經驗者跳至下一題。
所以這一題不是必答的。
那麼在輸入資料進spss時,不須答題的受試者欄位,請問應該要留空白?
還是要設value為0並在data view裡輸入0?
謝謝您
spss如果資料是空白的話，那他會是為missing data。建議你將未具經驗的受訪者設定為0；反之為1，之後在搭配選擇觀察值來分析
您好，
一般在多元迴歸分析裡，是可以設置一個依變項與多個自變項。
但是我是用複選題去問的，所以目前是依變項的部份有好幾個。(因為是把消費者的不同反應設成好幾個選項)。
然後我的自變項只有一個。
這種情形下還是用複迴歸嗎？
或是用交叉分析就可以了？
謝謝您
複選題基本上就只能作交叉分析～
在SPSS中先將複選題的集合定義好，之後做交叉分析比較方便
您好:
最近我因論文所需正在摸索mathematica這套軟體
但看到您的教學講義已經無法download
不知是否可以寄到下面信箱給我
lover5314@yahoo.com.tw
謝謝您
目前只釋出第三章
http://www.pws.stu.edu.tw/cydye/Chapter-3.pdf.zip
其他的部份還在排版中，有好的建議請提供給我！
感恩喔～
大大你好
本人只是在google search A.M.>= G.M. PROOF中
無意來到本站
雖然這兒像是沒有啥人
不過還是一個好網站
CM給你以示支持 :)
"雖然這兒像是沒有啥人"。
這不就是我太太說的話嗎？
Why Nobody Clicks Your Advertisement?
因為大家都是正常人，只有鬼才看得懂你這個部落格在寫什麼。
Orz
J-最愛石榴紅茶
哈哈 我是鬼~~而且我覺得寫得很棒呢
要幫你平反一下
一個還在為了Ph.d.努力的小鬼
太認真了～～po文時間是早上05:29
現在難得有博士班學生這麼認真的了～～
您好，我是在網路上搜尋而連結過來的，有些問題想請教~^^"
在您的網站有稍微解釋常態分布曲線如何使用Excel繪製，但是X軸希望能顯示是我的資料中的長度(mm)，而不是-4~+4之類的~@@
我已有自己使用NORMDIST算出數據，若選擇Excel中的XY散布圖，只能顯示出資料中有的數據，像是在兩側可能就會"斷掉"而不是一個完美的鐘形圖......
想請教要如何才能製作出X軸為長度，且在鐘形圖中央有Y軸的圖表呢?
(因為非相關科系，所以可能要煩請您稍為多加解釋細節了~~QQ 謝謝您~)
加兩個點，在將兩點連起來！
不然擬就email給我，我寄個範例給你吧！
你好，我是mathematica的初學者> <
想請問您...我想要使用各種分配
要用到Statistics package
請問我該如何加入這個package呢?
(我在mathematica library有下載到常態分配的package，那下載下來的檔案該如何處置呢?)
麻煩您了!!謝謝你!!
Mathematica6統計的Package好像就不需要用Need的方式呼叫了就直接可以使用。在Help找Statistics看看！有需要的可以討論～
要常態分配～
NormalDistribution[\[Mu],\[Sigma]]
PDF
xpdf=PDF[NormalDistribution[0,1],x];
Plot[xpdf,{x,-4,4}]
CDF
xcdf=CDF[NormalDistribution[0,1],x];
Plot[xcdf,{x,-4,4}]
了解常態分配的一些特性
Manipulate[
Plot[PDF[NormalDistribution[mu,sigma],x],{x,-21,21},
PlotRange->{{-21,21},{0,0.4}},
PlotStyle->Thickness[0.015]],{mu,1,5},{sigma,1,5}]
產生1000個標準常態隨機變數，並以直方圖表示
mydata=RandomReal[NormalDistribution[0,1],1000];
Histogram[mydata]
產生100個卡方自由度10的隨機變數並計算樣本平均數，重複該實驗10000次
data1=Mean@RandomReal[ChiSquareDistribution[2],200]&/@
Range[10000];
將上述資料做標準化轉換，並以直方圖匯出，驗證中央極限定理
Histogram[(data1-#[[1]])/
Sqrt[#[[2]]]&/@{#@data1&/@{Mean,Variance}}]
Manipulate的部份改一下會更清楚
Manipulate[
Plot[PDF[NormalDistribution[mu,sigma],x],{x,-21,21},
PlotRange->{{-21,21},{0,0.4}},PlotStyle->Thickness[0.015],
Filling->Axis],{mu,-5,5},{sigma,1,5}]
欣賞是一種心靈的學習成長
分享是喜悅和快樂泉源
很高興有緣在此認識相遇
共同分享這裡的喜悅
彼此共勉之
寫得太玄了～～
http://everythingonmyipod.blogspot.com/2010/01/nicholas-gunn-music-of-grand-canyon.html
可否拜託您幫我回答這個MATHEMATICA 曲線切線的問題
http://tw.knowledge.yahoo.com/question/question?qid=1510031900025&mode=w&from=question&recommend=0&.crumb=FSd7ynECnBR
拜託 謝謝您
作業多練習進步比較快！
你好!想請教你一下,當我把實驗數據從Mathematica軟體裡面叫出時,我的實驗曲線會變成藍色實線我要怎麼把我的藍色實線變換成虛線,另外我再用公式以Mathematica軟體去運算後所得到的模擬曲線會是紅色虛線,那我要怎麼把所模擬運算後的曲線變成實線,我弄了好久,都弄不出來,可否煩請大哥你賜教,很急~感謝你囉~
data={{1,0.9733},{2,0.8878},{3,8.2889},{4,5.6421},{5,9.7358},{6,7.8550},{7,8.1411},{8,8.4569},{9,4.8627},{10,5.6916}}
(*以最小平方法擬合資料*)
fit=LinearModelFit[data,x,x][x]
Plot[fit,{x,0,10},
PlotRange->{{0,12},{0.8Min@data[[All,2]],1.2Max@data[[All,2]]}},
(*擬合函數，紅線，粗細0.01*)
PlotStyle->{Red,Thickness[0.01]},
(*以下為實驗資料樣本點，藍色，虛線，兩線間空隔0.05，樣本點綠色，點數大小0.025*)
Epilog->{Thickness[0.01],Blue,Dashing[0.05],Line[data],Green,
PointSize[0.025],Point[data]}]
你好
假設 p=x*y*z(2x+5)
然後 p/(x*y*z)
但是 結果 不會是 (2x+5)
而市 x*y*z(2x+5)/x*y*z
有什麼方法可以讓真的把 x*y*z除掉 就只會看到(2x+5)而已?
謝謝
應該輸出就是2x+5，如果沒有那後面加個//Simplify
p/(x*y*z)//Simplify
謝謝您~ 請問你的E-mail 我想把解方程式的兩個不同結果跟你討論 不知是否ok
謝謝
chungyuandye在gmail
已經Mail給您了 再煩請您確認 謝謝
我是從 PTT 連過來的同學，感謝您的解說，真的幫了我很大的忙，真的是太感謝了！
還好以上的留言或回覆 我統統看不懂 不然我就是"鬼" 了
老師：您好！
我是EMBA的大姊〈黃同學〉。我有很大的問題：被電腦打敗了。我不會把奧運的data放到SPSS來操作。煩請老師寄給我，好嗎？很急！Thanks a lot！Happy New Year！
淑櫻E-mail：haung@ms.dsp.ks.edu.tw
老師您好, 前兩天有寄email給您問問題,我是要對list的elements做判斷然後取出我要的elements.
你教我一個方法
a = {{1}, {3}, {5, 6}, {2}, {4, 7}};
Do[Do[If[Length[a[[i]]] < 2, a = Delete[a, {{i}}]; Break[], a], {i,
Length[a]}], {i, Length[a]}]
但是有說這樣很沒效率可以改用Select或DeleteCases
確實...但現在我的狀況變成如下
a = {{{0.9330292`, 15.62140356585095`}, {0.9330292`,
15.643149500693339`}, {0.9330292`,
15.64315074986568`}}, {{0.9330293`,
15.600339001506518`}, {0.9330293`,
15.610411125774197`}, {0.9330293`,
15.636513912311804`}, {0.9330293`,
15.642874703177677`}, {0.9330293`,
15.64287505964573`}, {0.9330293`,
15.64287529104515`}}, {{0.9330294`,
15.590924622076912`}, {0.9330294`,
15.599960978631318`}, {0.9330294`,
15.631799083935137`}, {0.9330294`,
15.633517619820962`}, {0.9330294`,
15.642600328138865`}, {0.9330294`, 15.642600400164431`}}};
Do[Do[Do[Do[
If[Length[a[[i]]] > 2, a = Delete[a, {i, j}]; Break[], a], {j, 2,
Length[a[[i]]] - 1}], {j, 2, Length[a[[i]]] - 1}], {i,
Length[a]}], {i, Length[a]}]
a
我試過這是可行的, 但是我真正要跑的list是上千個點
run過一次等了幾十分鐘我還是把它給停了...@@
不曉得有沒有比較好的方法可以用select解決 感謝^^
{#[[1]], #[[-1]]} & /@ a
在Mathematica盡量不用Do While For
兩個問題請教戴老師
1.如果我要取出#[[2]]中的最大值 最小值 中間值，該打什麼指令?
2.Prob[L=300,R=1100,h=0.1]
如果要印出h=0.1~2(0.1為一個間隔)的圖該打什麼指令
非常感謝
data = RandomReal[{0, 10}, {5, 6}]
{{7.36756, 0.94386, 3.22745, 3.55061, 4.32153, 7.37545}, {1.91357,
5.82995, 7.67345, 7.52497, 2.59294, 3.70943}, {3.14569, 4.17957,
4.59748, 6.78449, 0.293061, 7.77275}, {4.8603, 4.50177, 6.95772,
7.11461, 1.48556, 5.04942}, {6.12433, 0.313357, 9.51655, 4.41046,
0.463432, 3.17241}}
{Max@#, Min@#, Mean@#, Median@#} & /@ data
{{7.37545, 0.94386, 4.46441, 3.93607}, {7.67345, 1.91357, 4.87405,
4.76969}, {7.77275, 0.293061, 4.46217, 4.38853}, {7.11461, 1.48556,
4.9949, 4.95486}, {9.51655, 0.313357, 4.00009, 3.79143}}
請問Prob的定義是如何？？
如果寫成Prob[L_,R_.h_]:=babl...
那這個圖可以這樣寫
{#,Prob[300,1100,#}&/@Range[0,2,0.1]
Prob[L=300,R=1100,h=#]&/@Range[0,2,0.1]
現在有個list較做a
a = {{{1, 55}, {1, 56}}, {{2, 54}, {2, 52}}, {{3, 56}, {3, 58}}, {{4,
52}, {4, 59}}, {{5, 51}, {5, 60}}, {{6, 50}, {6, 61}}}
裡面有sublist都是兩個elements,每個elements又有兩個elements
我現在要求是要留下a[[j]][[1,2]]i
而a[[j]][[2,2]]>a[[i]][[2,2]], for j>i
所以我用老師之前教我的方法寫了四個Do
Do[Do[Do[Do[
If[a[[j]][[1, 2]] > a[[i]][[1, 2]] ||
a[[j]][[2, 2]] < a[[i]][[2, 2]], a = Delete[a, {{j}}]; Break[],
a], {j, i + 1, Length[a]}], {j, i + 1, Length[a]}], {i,
Length[a] - 1}], {i, Length[a] - 1}]
結果是ok的
但是跟之前一樣碰到list很大就要跑很久
現在因為判斷的標準是要a的兩個elements來比較
我想不出來不用Do loop的方法, 麻煩老師解惑^^ 謝謝
a = {{{1, 55}, {1, 56}}, {{2, 54}, {2, 52}}, {{3, 56}, {3, 58}}, {{4,
52}, {4, 59}}, {{5, 51}, {5, 60}}, {{6, 50}, {6, 61}}};
b = {1}; i = 1;
Nest[({i, j} = #;
If[a[[i, 1, 2]] > a[[j, 1, 2]] && a[[i, 2, 2]] < a[[j, 2, 2]],
b = Insert[b, j, -1]; {i, j} = {j, j + 1}, {i, j} = {i,
j + 1}]; {i, j}) &, {1, 2}, Length@a - 1];
a[[#]] &@b
{{{1, 55}, {1, 56}}, {{4, 52}, {4, 59}}, {{5, 51}, {5, 60}}, {{6,
50}, {6, 61}}}
請問
Prob[L=300,R=1100,h=#]&/@Range[0,2,0.1]
可不可以不要印很多張圖 把所有圖畫在一張圖上?
Show[Prob[L=300,R=1100,h=#]&/@Range[0,2,0.1]]
不過我沒有這個函數，否則應該先將資料產生後再ListPlot出來這樣會比較好
你的需求應該要修改一下函數，不過你的h太多圖形複雜
Prob[LL_, RR_, hh_] :=
Block[{L = LL, R = RR, h = hh, P, \[Alpha]1, K, kpeq, Jdata, JPoint,
t, p1, p2, \[Theta]},
\[Alpha]1 = ArcTan[357/1100];
K[x_] := {(R - h)*Sin[x], (R - h)*Cos[x], 0};
P[\[Phi]_] := {L*Cos[\[Phi]], 0, 0};
kpeq[a_, b_] := Block[{v, slope, xy},
v = P[a][[1 ;; 2]] - K[b][[1 ;; 2]];
slope = v[[2]]/v[[1]];
slope*(x - P[a][[1]]) + P[a][[2]]
];
\[Theta] = 0;
JPoint[a_, b_] :=
Block[{x, aa = a, bb = b}, {x, Sqrt[R^2 - x^2]} /.
Solve[Sqrt[R^2 - x^2] == kpeq[aa, bb], x][[1]]];
Jdata = {#*Pi, ArcSin[JPoint[#*Pi, \[Alpha]1][[1]]/R]} & /@
Range[0, 1, 1/360];
(*
p1=ListLinePlot[Jdata,AspectRatio->1,AxesLabel->{"\[CurlyPhi]",
"\[Alpha]"},PlotStyle->Thickness[0.01],Ticks->{Range[0,Pi,Pi/2]}];
*)
p1 = Jdata;
t[\[Alpha]_, \[Phi]_] := ((R*
Cos[\[Alpha]])/(\[Sqrt]((R*Sin[\[Alpha]]*Cos[\[Theta]] -
L*Cos[\[Phi]])^2 + (R*Sin[\[Alpha]]*Sin[\[Theta]] -
L*Sin[\[Phi]])^2 + (R*Cos[\[Alpha]])^2))*((-R*
Sin[\[Alpha]]*Cos[\[Theta]]) (L*Cos[\[Phi]] -
R*Sin[\[Alpha]]*Cos[\[Theta]]) + (-R*Sin[\[Alpha]]*
Sin[\[Theta]]) (L*Sin[\[Phi]] -
R*Sin[\[Alpha]]*Sin[\[Theta]]) + (-R*Cos[\[Alpha]]) (-R*
Cos[\[Alpha]]))/((\[Sqrt]((L*Cos[\[Phi]] -
R*Sin[\[Alpha]]*Cos[\[Theta]])^2 + (L*Sin[\[Phi]] -
R*Sin[\[Alpha]]*Sin[\[Theta]])^2 + (-R*
Cos[\[Alpha]])^2)) (1100))*1/((R*Sin[\[Alpha]]*
Cos[\[Theta]] -
L*Cos[\[Phi]])^2 + (R*Sin[\[Alpha]]*Sin[\[Theta]] -
L*Sin[\[Phi]])^2 + (R*Cos[\[Alpha]])^2));
(*
p2=ListLinePlot[{1100*#[[2]],t[#[[2]],#[[1]]]}&/@Jdata,
PlotStyle->Thickness[0.01],AspectRatio->1,AxesLabel->{"x","t"}];
Grid[{{p1,p2}}]
*)
p2 = {1100*#[[2]], t[#[[2]], #[[1]]]} & /@ Jdata;
{p1, p2}
]
aa = Prob[300, 1100, #] & /@ Range[0, 2, 0.1];
Grid[{{ListLinePlot[aa[[All, 1]], Frame -> True, Axes -> False,
ImageSize -> 250]}, {ListLinePlot[aa[[All, 2]], Frame -> True,
Axes -> False, ImageSize -> 250]}}]
建議h值不用跑太多，加個圖例，這樣整個圖形就一目了然
LineStylesFromTooltips[plot_Graphics] :=
Cases[plot,
Tooltip[{s__, l_Line}, tt_] :>
Grid[{{Graphics[
Flatten[{s, Thickness[0.15], Line[{{0, 0}, {1, 0}}]}],
ImageSize -> {24, 8}, AspectRatio -> 8/24, ImagePadding -> 0],
tt}}], Infinity]
LineLegendFromTooltips[plot_Graphics] :=
Labeled[plot,
Style[Column[LineStylesFromTooltips[plot], Left], "TR",
ShowStringCharacters -> False], {{Right, Top}}]
mydata1 =
Tooltip[bb[[#, 1]],
"h=" <> ToString[NumberForm[(# - 1)*0.5, {2, 1}]]] & /@
Range[5];
LineLegendFromTooltips[
ListLinePlot[mydata1, PlotStyle -> Thickness[0.01], Frame -> True,
Axes -> False]]
mydata2 =
Tooltip[bb[[#, 2]],
"h=" <> ToString[NumberForm[(# - 1)*0.5, {2, 1}]]] & /@
Range[5];
LineLegendFromTooltips[
ListLinePlot[mydata2, PlotStyle -> Thickness[0.01], Frame -> True,
Axes -> False]]
老師您好
想請問把 數值A(Ic) Export到一個文字檔(aaa.txt)時
我們都知道如果當我第二次再把"數值A" Export 到aaa.txt時會覆蓋掉原先記事本裡面的值
想請問可否有辦法不要覆蓋而是一直存下去
我目前的問題是，我的"數值A"是因為manipulate方程式的變數而產生的，所以每一次我都只會在我記事本裡面看到一個值
Export["C:\\ddd\\aaa.txt", Ic, "Table"]
用>>>
For[i = 1, i <= 10, i >>> aaa.txt, i++]
檔案會存到你電腦裡"我的文件"
Prob[LL_, RR_, hh_] :=
Block[{L = LL, R = RR, h = hh, P, \[Alpha]1, K, kpeq, Jdata, JPoint,
t, p1, p2, \[Theta]},
\[Alpha]1 = ArcTan[259.1/1100];
K[x_] := {(R - h)*Sin[x], (R - h)*Cos[x], 0};
P[\[Phi]_] := {L*Cos[\[Phi]], 0, 0};
kpeq[a_, b_] := Block[{v, slope, xy},
v = P[a][[1 ;; 2]] - K[b][[1 ;; 2]];
slope = v[[2]]/v[[1]];
slope*(x - P[a][[1]]) + P[a][[2]]
];
\[Theta] = 0;
JPoint[a_, b_] :=
Block[{x, aa = a, bb = b}, {x, Sqrt[R^2 - x^2]} /.
Solve[Sqrt[R^2 - x^2] == kpeq[aa, bb], x][[1]]];
Jdata = {#*Pi, ArcSin[JPoint[#*Pi, \[Alpha]1][[1]]/R]} & /@
Range[0, 1, 1/361];
t[\[Alpha]_, \[Phi]_] :=
2*((R*Cos[\[Alpha]])/(\[Sqrt]((R*Sin[\[Alpha]]*Cos[\[Theta]] -
L*Cos[\[Phi]])^2 + (R*Sin[\[Alpha]]*Sin[\[Theta]] -
L*Sin[\[Phi]])^2 + (R*Cos[\[Alpha]])^2))*((-R*
Sin[\[Alpha]]*Cos[\[Theta]]) (L*Cos[\[Phi]] -
R*Sin[\[Alpha]]*Cos[\[Theta]]) + (-R*Sin[\[Alpha]]*
Sin[\[Theta]]) (L*Sin[\[Phi]] -
R*Sin[\[Alpha]]*Sin[\[Theta]]) + (-R*Cos[\[Alpha]]) (-R*
Cos[\[Alpha]]))/((\[Sqrt]((L*Cos[\[Phi]] -
R*Sin[\[Alpha]]*Cos[\[Theta]])^2 + (L*Sin[\[Phi]] -
R*Sin[\[Alpha]]*Sin[\[Theta]])^2 + (-R*
Cos[\[Alpha]])^2)) (1100))*1/((R*Sin[\[Alpha]]*
Cos[\[Theta]] -
L*Cos[\[Phi]])^2 + (R*Sin[\[Alpha]]*Sin[\[Theta]] -
L*Sin[\[Phi]])^2 + (R*Cos[\[Alpha]])^2));
p2 = ListLinePlot[{1100*#[[2]], t[#[[2]], #[[1]]]} & /@ Jdata,
PlotStyle -> Thickness[0.01], AspectRatio -> 1,
AxesLabel -> {"x", "t"}];
Grid[{{p1, p2}}]]
Prob[300, 1100, 0]
t是一個積分，當Pi等於0時，t應該是零，可是畫的圖t的最小值卻不是零?
檢查了好久都沒發現異樣在哪裡 如果老師可以指點一下 真是太感謝了
檔案以寄出。
抱歉我沒收到耶
可以麻煩老師在寄一次嗎?
非常感謝
Prob[LL_, RR_, hh_] :=
Block[{L = LL, R = RR, h = hh, P, \[Alpha]1, K, kpeq, Jdata, JPoint,
t, p1, p2, \[Theta], ss, tt, g, p3},
\[Alpha]1 = ArcTan[357/1100];
K[x_] := {(R - h)*Sin[x], (R - h)*Cos[x], 0};
P[\[Phi]_] := {L*Cos[\[Phi]], 0, 0};
kpeq[a_, b_] := Block[{v, slope, xy},
v = P[a][[1 ;; 2]] - K[b][[1 ;; 2]];
slope = v[[2]]/v[[1]];
slope*(x - P[a][[1]]) + P[a][[2]]
];
\[Theta] = 0;
JPoint[a_, b_] :=
Block[{x, aa = a, bb = b}, {x, Sqrt[R^2 - x^2]} /.
Solve[Sqrt[R^2 - x^2] == kpeq[aa, bb], x][[1]]];
Jdata = {#*Pi, ArcSin[JPoint[#*Pi, \[Alpha]1][[1]]/R]} & /@
Range[0, 1, 1/360];
(*
p1=ListLinePlot[Jdata,AspectRatio->1,
AxesLabel->{"\[CurlyPhi]","\[Alpha]"},PlotStyle->Thickness[0.01],
Ticks->{Range[0,Pi,Pi/2]}];
*)
p1 = Jdata;
t[\[Alpha]_, \[CurlyPhi]_] := (
R^2 Cos[\[Alpha]] (R -
L Cos[\[Theta] - \[CurlyPhi]] Sin[\[Alpha]]))/(
1100 (L^2 + R^2 -
2 L R Cos[\[Theta] - \[CurlyPhi]] Sin[\[Alpha]])^2);
ss = 1100*#[[2]] & /@ Jdata;
tt = t[#[[2]], #[[1]]] & /@ Jdata;
g = Interpolation[{Jdata[[#, 1]], tt[[#]]} & /@
Range[Length@Jdata]];
(*
p2=ListLinePlot[{1100*#[[2]],t[#[[2]],#[[1]]]}&/@Jdata,
PlotStyle->Thickness[0.01],AspectRatio->1,AxesLabel->{"x","t"}];
Grid[{{p1,p2}}]
*)
p2 = {ss[[#]], Integrate[g[x], {x, 0, Jdata[[#, 1]]}]} & /@
Range[Length@Jdata];
p3 = Plot[Evaluate@Integrate[g[x], {x, 0, y}], {y, 0, Pi}][[1,
1, -1, -1, 1]];
{p1, p2, p3}
];
LineStylesFromTooltips[plot_Graphics] :=
Cases[plot,
Tooltip[{s__, l_Line}, tt_] :>
Grid[{{Graphics[
Flatten[{s, Thickness[0.15], Line[{{0, 0}, {1, 0}}]}],
ImageSize -> {24, 8}, AspectRatio -> 8/24, ImagePadding -> 0],
tt}}], Infinity]
LineLegendFromTooltips[plot_Graphics] :=
Labeled[plot,
Style[Column[LineStylesFromTooltips[plot], Left], "TR",
ShowStringCharacters -> False], {{Right, Top}}]
bb = Prob[300, 1100, #] & /@ Range[0, 2, 0.5];
mydata1 =
Tooltip[bb[[#, 1]],
"h=" <> ToString[NumberForm[(# - 1)*0.5, {2, 1}]]] & /@
Range[5];
LineLegendFromTooltips[
ListLinePlot[mydata1, PlotStyle -> Thickness[0.01], Frame -> True,
Axes -> False]]
mydata2 =
Tooltip[bb[[#, 2]],
"h=" <> ToString[NumberForm[(# - 1)*0.5, {2, 1}]]] & /@
Range[5];
LineLegendFromTooltips[
ListLinePlot[mydata2, PlotStyle -> Thickness[0.01], Frame -> True,
Axes -> False]]
mydata3 =
Tooltip[bb[[#, 3]],
"h=" <> ToString[NumberForm[(# - 1)*0.5, {2, 1}]]] & /@
Range[5];
LineLegendFromTooltips[
ListLinePlot[mydata3, PlotStyle -> Thickness[0.01], Frame -> True,
Axes -> False]]
{bb[[1, 1, 1]], bb[[1, 2, 1]]}
{{0, 0.313821}, {345.203, 0}}
請問Mathmetica的點座標(很多組)
可以直接轉成excel檔或是matlab檔嗎?
有這種指令嗎?
Export["c:\\test.csv",RandomReal[NormalDistribution[0, 1], {100, 2}]]
檔案會存在你的C下面
挖 真是太好用了!!謝謝老師 我之前還用c來轉轉半天
不過如我要輸出Jdata 的值
Jdata = {#*2 Pi, 180/Pi*ArcSin[JPoint[#*2 Pi, \[Alpha]1][[1]]/R]} & /@
Range[0, 1, 1/360];
指令打這樣
Export["c:\\test.csv", Jdata[{360, 2}]]
不行耶>"<
Jdata的定義應該是一個陣列。你要的Jdata[{360,2}]是什麼？
應該直接Export["c:\\test.csv",Jdata]
我有試過
Export["c:\\test.csv",Jdata]
可是只會跑出 "Jdata"一個字
所以我才加[{360,2}]
Jdata的陣列有360列
不過Export["c:\\test.csv",Jdata]不行耶
謝謝
因為Jdara是不是全域函數，所以你直接夏Jdata當然跑不出來。
L = 300; R = 100; \[Alpha]1 = ArcTan[357/1100];
K[x_] := {(R - h)*Sin[x], (R - h)*Cos[x], 0};
P[\[Phi]_] := {L*Cos[\[Phi]], 0, 0};
kpeq[a_, b_] := Block[{v, slope, xy},
v = P[a][[1 ;; 2]] - K[b][[1 ;; 2]];
slope = v[[2]]/v[[1]];
slope*(x - P[a][[1]]) + P[a][[2]]
];
\[Theta] = 0;
JPoint[a_, b_] :=
Block[{x, aa = a, bb = b}, {x, Sqrt[R^2 - x^2]} /.
Solve[Sqrt[R^2 - x^2] == kpeq[aa, bb], x][[1]]];
h = 0.1;
Jdata = N@{#*Pi, ArcSin[JPoint[#*Pi, \[Alpha]1][[1]]/R]} & /@
Range[0, 1, 1/360];
Export["c:\\test.csv", Jdata]
戴老師您好，我是潘家群，謝謝老師慷慨地給我高等統計學講義。為了我的上課需要，已將它改成教育心理統計學(最後兩章沒有教到，所以沒有更改)。現在回傳給戴老師 (放在 bandongo 空間: http://www.badongo.com/cn/file/26173108)，這份講義依然根懅老師您的規範來使用。若有人需要這份講義，也可以直接與我聯繫 (panjc1019@gmail.com)。目前我已畢業於交大統計所，正在等待一個博士後的機會。謝謝戴老師的幫助，謝謝!
很久沒有使用MATHEMATICA，全部都忘光光；還好有老師的網站可以查資料。
久久上來看一次，謝謝老師~~
相見恨晚
這網頁真的太棒了!!!!!!!!!!!!!!!!!!!!!!!!!!!!
發佈留言
Please leave your name and tell me what you thought about this site. Comments, suggestions and views are welcomed.
如果這篇文章對你有幫助，那請留個訊息給我～