MATLAB 板


LINE

謝謝celestialgod! 簡單來說,我有三個公式想要記算,分別是S1, F2, A2的算式 但各自公式成立之前都有一些不一樣的限制式需要被滿足, 如果限制式有 > 0才計算數值,不然就顯示0。 所以, 1. 在S1的前面應該要有: if......(S1成立的條件式) S1 = ...... else S1 = 0 if......(F2成立的條件式) F2 = ...... else F2 = 0 if......(A2成立的條件式) A2 = ...... else A2 = 0 2. 這樣一輪下來應該就會有S1, F2, A2的值 (有數字或是0) 接著去比較在b = 0, D = 0.1的狀況下哪個是Max並且print name出來 然後迴圈後在b = 0.1, D = 0.1的狀況下再比較一次 ...經過很多輪 然後迴圈後在b = 1, D = 1的狀況下再比較一次 3. 最後希望可以整理出在不同b跟D的組合下,什麼時候S1, F2, A2會是max 然後想要畫出一個圖http://i.imgur.com/NNTHMT9.jpg 像右邊圖,在不同的組合下面max的值是哪一個 不好意思再麻煩各位幫我解惑了 > < 謝謝! ※ 引述《celestialgod (天)》之銘言: 有些問題看得不是很懂... 有問題推文問再回你吧~"~ q=0.9; a=0.9; c0=0.2; c1=0.01; X=2^(1/2); [b, D] = meshgrid(0:0.1:1, 0.1:0.1:1); S1 = (q*(a+b(:))-2*c0).^2/(16*(q-c1)); F2 = (q^2*(2*D(:) + a - b(:)) .* (a*q - 2*c0 + b(:)*q))/(16*(c1 - q).^2); A2 = ((X*q*(q*(a + b(:))-2*c0)).^2)/((4*(q-c1)).^3); Y = [S1, F2, A2]; Y(q - b(:) > 0 & q*(q + b(:)) - 2*c0 > 0) = 0; if q - c1 <= 0 Y(:) = 0; end [Z, code] = max(Y, [], 2); names = {'S1', 'F2', 'A2'}; maxNames = arrayfun(@(x) names{x}, code, 'uni', false); 儲存的話,要同時存chr跟數字的話,就用open跟fprintf解決吧,網路上很多 存數字就用dlmwrite就好 ※ 引述《kevinjean (Kevin)》之銘言: : 各位版上前輩好, : 我想要請教幾個問題希望大家可以幫助我 > < : --- : q=0.9, a=0.9, c0=0.2, c1=0.01, X=2^(1/2) : for b=0:0.1:1 : for D=0.1:0.1:1 : if q - b > 0 & q - c1 > 0 & q*(q + b) - 2*c0 > 0 : S1 = (q*(a+b)-2*c0)^2/(16*(q-c1)) : F2 = (q^2*(2*D + a - b)*(a*q - 2*c0 + b*q))/(16*(c1 - q)^2) : A2 = ((X*q*(q*(a + b)-2*c0))^2)/((4*(q-c1))^3) : else : S1 = 0 : F2 = 0 : A2 = 0 : end : end : Y = [S1 F2 A2] : Z = max(Y) : fprintf('Max = \n',Z) : end : --- : 總而言之我希望調控b跟D分別跑出一組一組的S1, F2, A2, : 並且在每一組內找出這三個中最大的值,顯示對應的代號 : 所以應該會有10 x 9 = 90組Max的值, : 可能像是 : Max = S1 : Max = S1 : Max = S1 : Max = S1 : Max = F2 : Max = F2 : Max = A2 : Max = A2 : Max = A2 : Max = A2 : Max = A2 : Max = A2 : ... : 目前想請大家幫的忙有幾個: : 1. 我無法在每一個b - D組合得到一個Max值 : 跑完之後他會先出現很多組S1, F2, A2, : 然後在最後一組之後出現: : Y = [xxx yyy zzz] : Z = yyy (假設是最大值) : 而不是每一筆資料都是一個Max值 : 2. 我無法得到代號,而是數值 : 我知道很直覺的我的矩陣裡面放的是S1, F2, A2計算出來的值, : 但我應該要怎麼寫才可以讓Max出來的是代號呢? : 3. 我無法儲存最後得到的Max值 : 假設順利得到了90組Max值,我要怎麼儲存他? : 因為目前好像就是會在Matlab出現90行的文字, : 但我希望可以把它儲存為 : (b1, D1, Max1) : (b2, D1, Max2) : (b3, D1, Max3) : ... : 再請大家幫助我了,謝謝各位 > < -- --



※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.109.74.87
※ 文章網址: https://webptt.com/m.aspx?n=bbs/MATLAB/M.1449458146.A.187.html ※ 編輯: celestialgod (140.109.74.87), 12/07/2015 11:18:04 --



※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.113.250.42
※ 文章網址: https://webptt.com/m.aspx?n=bbs/MATLAB/M.1449461619.A.74F.html ※ 編輯: kevinjean (140.113.250.42), 12/07/2015 12:14:06
1F:→ celestialgod: 1. 計算後,我再取條件等於0,不是一樣? 12/07 12:17
2F:→ celestialgod: 如果你要個別條件,在個別加就好 12/07 12:18
3F:→ celestialgod: 2. 我不知道你在寫什麼... 12/07 12:18
4F:→ kevinjean: 例如(b,D,max)=(0,0.1,S1);(0.1,0.1,S1);(0.2,0.1,S1) 12/07 12:21
5F:→ kevinjean: (0.3,0.1,S1);...;(0.5,0.5,A2);(0.6,0.5,A2);... 12/07 12:21
6F:→ kevinjean: ...;(1,0.9,F2);(1,1,F2) 這樣的輸出 12/07 12:22
7F:→ celestialgod: 所以我前面不就說可以用fprintf做? 12/07 12:32
8F:→ celestialgod: 前面程式OUTPUT是錯的,你也沒發現... 12/07 13:07
9F:→ celestialgod: 你到底有沒有跑過我寫的程式,再來問問題? 12/07 13:07
10F:→ kevinjean: 有我有跑過了,謝謝!目前先看到的狀況就是有一連串 12/07 13:19
11F:→ kevinjean: S1, F2, A2。我還沒有先對output的內容,我會再研究 12/07 13:19
12F:→ kevinjean: 一下,謝謝celestialgod 12/07 13:20
13F:→ celestialgod: 所以output應該都給你了,我不知道你後面要幹嘛就是 12/07 13:23
14F:推 sunev: 限制式只限制b , 解出來剛好是一個開區間。 12/07 13:38
15F:→ sunev: S1 F2 A2又剛好只是b D的二次式,應該可以手算才是 12/07 13:39







like.gif 您可能會有興趣的文章
icon.png[問題/行為] 貓晚上進房間會不會有憋尿問題
icon.pngRe: [閒聊] 選了錯誤的女孩成為魔法少女 XDDDDDDDDDD
icon.png[正妹] 瑞典 一張
icon.png[心得] EMS高領長版毛衣.墨小樓MC1002
icon.png[分享] 丹龍隔熱紙GE55+33+22
icon.png[問題] 清洗洗衣機
icon.png[尋物] 窗台下的空間
icon.png[閒聊] 双極の女神1 木魔爵
icon.png[售車] 新竹 1997 march 1297cc 白色 四門
icon.png[討論] 能從照片感受到攝影者心情嗎
icon.png[狂賀] 賀賀賀賀 賀!島村卯月!總選舉NO.1
icon.png[難過] 羨慕白皮膚的女生
icon.png閱讀文章
icon.png[黑特]
icon.png[問題] SBK S1安裝於安全帽位置
icon.png[分享] 舊woo100絕版開箱!!
icon.pngRe: [無言] 關於小包衛生紙
icon.png[開箱] E5-2683V3 RX480Strix 快睿C1 簡單測試
icon.png[心得] 蒼の海賊龍 地獄 執行者16PT
icon.png[售車] 1999年Virage iO 1.8EXi
icon.png[心得] 挑戰33 LV10 獅子座pt solo
icon.png[閒聊] 手把手教你不被桶之新手主購教學
icon.png[分享] Civic Type R 量產版官方照無預警流出
icon.png[售車] Golf 4 2.0 銀色 自排
icon.png[出售] Graco提籃汽座(有底座)2000元誠可議
icon.png[問題] 請問補牙材質掉了還能再補嗎?(台中半年內
icon.png[問題] 44th 單曲 生寫竟然都給重複的啊啊!
icon.png[心得] 華南紅卡/icash 核卡
icon.png[問題] 拔牙矯正這樣正常嗎
icon.png[贈送] 老莫高業 初業 102年版
icon.png[情報] 三大行動支付 本季掀戰火
icon.png[寶寶] 博客來Amos水蠟筆5/1特價五折
icon.pngRe: [心得] 新鮮人一些面試分享
icon.png[心得] 蒼の海賊龍 地獄 麒麟25PT
icon.pngRe: [閒聊] (君の名は。雷慎入) 君名二創漫畫翻譯
icon.pngRe: [閒聊] OGN中場影片:失蹤人口局 (英文字幕)
icon.png[問題] 台灣大哥大4G訊號差
icon.png[出售] [全國]全新千尋侘草LED燈, 水草

請輸入看板名稱,例如:WOW站內搜尋

TOP