| 查看: 1379 | 回復(fù): 14 | ||
燕尾蝶12新蟲 (小有名氣)
|
[求助]
急求各位大神,看下我這個(gè)編程有什么問題,萬分感謝 已有1人參與
|
|
這段程序作出了兩個(gè)三維圖形,是k和alpha關(guān)于beta和delta的。。。![]() ![]() ![]() ![]() clc;clear c=1;q=10;p=5.5; delta=0.10:0.01:1; beta=0.55:0.005:1; for i=1:length(delta) for j=1:length(beta) t=(1-delta(i)+beta(j)*delta(i))/beta(j)/delta(i)/q; m=beta(j)*delta(i)*q; u=m*c*t-(1-delta(i))*(p-c); K(i,j)=1/3/m*3^(1/3)*(u^2*(9*m+3^(1/2)*(u^2+27*m^2)^(1/2)))^(1/3)-1/3*u^2/m*3^(2/3)/(u^2*(9*m+3^(1/2)*(u^2+27*m^2)^(1/2)))^(1/3) alpha(i,j)=((K(i,j)*m+u)/2/K(i,j)/t/m; test(i,j)=alpha(i,j)-beta(j)*delta(i); end end mesh(delta,beta,K) xlabel('\delta') ylabel('\beta') zlabel('K') figure mesh(delta,beta,alpha) xlabel('\delta') ylabel('\beta') zlabel('\alpha') |
禁蟲 (正式寫手)
|
本帖內(nèi)容被屏蔽 |
禁蟲 (正式寫手)
|
本帖內(nèi)容被屏蔽 |
新蟲 (小有名氣)
新蟲 (小有名氣)
新蟲 (小有名氣)
新蟲 (小有名氣)
|
不知道我表達(dá)的是否清楚,其實(shí)我最終的目的是除了做出三維圖,還有就是判斷test(i,j)的正負(fù),您能再幫我看下下面的程序嗎? 在這個(gè)程序中test(i,j)正負(fù)不定,所以才想著要使k和alpha在[0,1]范圍內(nèi),這樣可以篩掉不滿足條件的delta(i)和beta(j)。不知道這種思路對(duì)不對(duì)?不知道您有更好的辦法嗎?萬分感謝。。。。 ![]() ![]() ![]() ![]() ![]() clc;clear c=1;q=10;p=5.5; delta=0.10:0.01:1; beta=0.55:0.005:1; for i=1:length(delta) for j=1:length(beta) t=(1-delta(i)+beta(j)*delta(i))/beta(j)/delta(i)/q; m=beta(j)*delta(i)*q; u=m*c*t-(1-delta(i))*(p-c); K(i,j)=1/3/m*3^(1/3)*(u^2*(9*m+3^(1/2)*(u^2+27*m^2)^(1/2)))^(1/3)-1/3*u^2/m*3^(2/3)/(u^2*(9*m+3^(1/2)*(u^2+27*m^2)^(1/2)))^(1/3); alpha(i,j)=(K(i,j)*m+u)/2/K(i,j)/t/m; test(i,j)=alpha(i,j)-delta(i); end end mesh(delta,beta,K) xlabel('\delta') ylabel('\beta') zlabel('K') set(gca,'zlim',[0,1]) figure mesh(delta,beta,alpha) xlabel('\delta') ylabel('\beta') zlabel('\alpha') set(gca,'zlim',[0,1]) 再次感謝。。。。 ![]() ![]() ![]() |
鐵桿木蟲 (正式寫手)
木頭蟲子
新蟲 (小有名氣)
禁蟲 (正式寫手)
|
本帖內(nèi)容被屏蔽 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 考研生物與醫(yī)藥調(diào)劑 +3 | 鐵憨憨123425 2026-03-31 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿北化085600材料專碩275|有文章專利|求調(diào)劑 +12 | Micky11223 2026-03-25 | 13/650 |
|
|
[考研] 材料與化工調(diào)劑一志愿大連海事085600,349 +6 | 吃的不少 2026-03-30 | 6/300 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +9 | @taotao 2026-03-29 | 9/450 |
|
|
[考研] 合肥區(qū)域性重點(diǎn)一本招收調(diào)劑 +4 | 6266jl 2026-03-30 | 4/200 |
|
|
[考研] 317求調(diào)劑 +8 | 十閑wx 2026-03-24 | 8/400 |
|
|
[考研] 求收留 +6 | 1943443204 2026-03-28 | 6/300 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +9 | hanamiko 2026-03-27 | 9/450 |
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +18 | Mr. Z 2026-03-25 | 19/950 |
|
|
[考研] 328求調(diào)劑 +8 | 嗯滴的基本都 2026-03-27 | 8/400 |
|
|
[考研] 考研調(diào)劑 +5 | Sanmu-124 2026-03-26 | 5/250 |
|
|
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷 +8 | cccchenso 2026-03-29 | 8/400 |
|
|
[考研] 299求調(diào)劑 +10 | 15188958825 2026-03-25 | 10/500 |
|
|
[考研] 327求調(diào)劑 +6 | 汲亦昊 2026-03-29 | 6/300 |
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|