版塊導航: 正在加載中...

登錄注冊

應《網絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯網跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

wcwl888

金蟲 (著名寫手)

中華人民共和國第一任航母艦長

應助: 12 (小學生)
金幣: 2639.4
散金: 210
紅花: 9
沙發(fā): 2
帖子: 1217
在線: 169.2小時
蟲號: 843785
注冊: 2009-09-08
性別: GG
專業(yè): 信息系統(tǒng)與管理

[求助] 菜鳥求助代碼運用已有1人參與

親愛的大神們：
各位好，學弟近期寫論文建模需要用到編程求解，但不懂運用何種軟件，最近看到一篇文獻里使用的模型和我的比較相似，苦于看不懂不知道該用什么軟件也不知道如何運行，故請教大神指點迷津，下面是C++源代碼：（附件里是模型需要的參數數值和運行結果，懂的大神可以參考一下）
#include<iostream.h>
void main ()
{ int i,j;
double a,b,c,1,v,hd,hm,n:
double w1,w2,w3,w11,w22,q1,lr11,lr12,lr1,lr21,q2,lr22,w,q,lr2;
cout<<"請依次輸入參數值:";
cin>>a>>b>>c>>l>>v>>hd>>hm>>n;
for(i=j=1;i<=10,j<=10;i++,j++)
{
wl=0。5*a+0。25*b*1*(2*i-1)+0。5*v-0。25*hd*1+0。5*1*hm*(i-(j-l)*n);
w2=a/3+(b*1*i-hd*1)/3+2*v/3+2*l*hm*(i-(j-1)*n)/3;
w3=a i*b*l+b*l;
if(wl<=0) w11=0;
else if(wl>=a+i*b*l-b*l) w11=a+i*b*l-b*l;
else w11=w1;
q1=l*(a/(2*c)+b*l/(4*c)*(2*i-l)-w11/(2*c)-hd*l/(4*c));
lr21=w11*q1-v*q1-l*hm*(i-(j-1)*n)*q1;
lr11=((a-w11+hd*(i-l)*l)*(a-w11+hd*(i-l)*l)*l+(b-hd)*(a-w11+hd*(i-1)*l)*l*l*(2*i-1)+((b-hd)*(b-hd)*l*l*l*(3*i*i-3*i+1))/3)/(4*c);
if(w2<a+i*b*l-b*l) w22=a+i*b*l-b*l;
else if(w2>=a+i*b*l-hd*l) w22=a+i*b*l-hd*l;
e1se w22=w2;
q2=((a+b*l*i-hd*l-w22)*(a+b*l*i-hd*l-w22))/(4*c*(b-hd));
lr22=w22*q2-v*q2-l*hm*(i-(j-1)*n)*q2;
lr12=((a-w22+hd*(i-1)*l)*(a-w22+hd*(i-1)*l)*(a-w22+b*i*l-hd*l)/(b-hd)+(b-hd)*(-w22+hd*(i-1)*l)*(a-w22+b*i*l-hd*l)/(b-hd)*(w22+b*i*l+hd*l-a-2*hd*i*l)/(b-hd)+(b-hd)*(b-hd)*(a-w22+b*i*l-hd*l)/(b-hd)*(((w22-a-hd*l*(i-1))*(w22-a-hd*l*(i-l)))/((b-hd)*(b-hd))+(w22-a-hd*l*(i-1))/(b-hd)*i*l+i*i*l*l)/3)/(4*c);
if(lr21>=lr22) {w=w11;q=q1:lr2=lr21;lr1=lr11;}
else {w=w22;q=q2;lr2=lr22;lr1=lr12;}
cout<<w<<" "<<q<<" "<<lr1<<lr2<<endl;
}

不知道代碼寫的有沒有錯誤，請懂的大神看看。幫忙備注一下每行代碼都是什么意思，該如何帶入數值進行運算出結果，如果運行比較困難的是否可以將這個代碼轉換成matlab程序語言，我學過一點matlab，但是C++代碼不會運行，如果兩種都會的大神還望不吝賜教，提供幫助的大神金幣雙手全奉供！比較急，小弟在此先拜謝了！

參數.png

計算結果.png

回復此樓

» 猜你喜歡

» 本主題相關價值貼推薦，對您同樣有幫助:

用lammps做模擬時如將新修改的勢函數讓lammps讀取運行已經有9人回復
菜鳥求助！關于柱芳烴的幾個問題已經有7人回復
求助--Openmpi-1.8集群編譯問題已經有37人回復
新手求助一個quantum-espresso的使用問題！已經有13人回復
菜鳥求助，老師給了一個分子讓優(yōu)化，可總是錯誤，求大神指點已經有22人回復
求助LAMMPS兩次運算結果不一致問題！已經有11人回復
編程菜鳥求大神引導入門。已經有14人回復
菜鳥求助用matlab如何畫三維圖已經有3人回復
求助用matlab 多元非線性擬合2 已經有21人回復
abaqus菜鳥求助已經有5人回復
LDPC碼的代碼求助已經有3人回復
sensors 投稿求助已經有10人回復
吐槽下量化新手還是TURBOMOLE使用者的苦逼之處已經有19人回復
菜鳥求助~matlab怎么樣讓一個1xn的向量中的每個元素是個矩陣？已經有11人回復
Fortran中數組的使用已經有4人回復
純菜鳥，關于icem最基礎問題已經有24人回復
求助：linux 國嵌加密部分（視頻，PPT ，代碼）已經有4人回復
求助大俠幫忙，有關matlab模擬數據已經有29人回復
菜鳥請教 asian journal of chemistry 關于DOI 已經有8人回復
【求助】LAMMPS-運行報錯已經有9人回復
【求助】用fortran怎么描述下面這樣的情況，我是新手，請大家?guī)兔?/a> 已經有8人回復

佛說：說了就是沒說，沒說就是說了，施主你又何必問呢！阿彌陀佛～

1樓 2015-03-19 21:25:12

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經驗: +2
計算強帖: 2
應助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準晶和納米晶
管轄: 計算模擬

【答案】應助回帖

感謝參與，應助指數 +1

你把源碼疏導c里面就可以了，注意下格式，你這直接從pdf轉的肯定錯誤多

贊一下

回復此樓

showmethemoney

2樓2015-03-20 12:00:04

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

wcwl888

金蟲 (著名寫手)

中華人民共和國第一任航母艦長

應助: 12 (小學生)
金幣: 2639.4
散金: 210
紅花: 9
沙發(fā): 2
帖子: 1217
在線: 169.2小時
蟲號: 843785
注冊: 2009-09-08
性別: GG
專業(yè): 信息系統(tǒng)與管理

引用回帖:

2樓: Originally posted by csgt0 at 2015-03-20 12:00:04
你把源碼疏導c里面就可以了，注意下格式，你這直接從pdf轉的肯定錯誤多

您能把這個代碼轉換成MATLAB語言嗎？只要轉換成matlab里可以運行的就可以了，金幣全給您！

贊一下

回復此樓

佛說：說了就是沒說，沒說就是說了，施主你又何必問呢！阿彌陀佛～

3樓2015-03-20 12:03:52

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經驗: +2
計算強帖: 2
應助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準晶和納米晶
管轄: 計算模擬

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
wcwl888: 金幣+100, ★★★很有幫助, 你有這份心就很滿足了，還沒驗證是否可行，金幣全送您了，有什么問題再來問您，非常感謝！ 2015-03-20 20:44:42

引用回帖:

3樓: Originally posted by wcwl888 at 2015-03-20 12:03:52
您能把這個代碼轉換成MATLAB語言嗎？只要轉換成matlab里可以運行的就可以了，金幣全給您！

...

公式只能你自己對了

CODE:

clc

clear

a=40;

b=1;

c=0.25;

l=10;

v=4;

hd=0.8;

hm=0.8;

n=1;

for i=1:10

    j=i;

            w1=0.5*a+0.25*b*1*(2*i-1)+0.5*v-0.25*hd*1+0.5*1*hm*(i-(j-l)*n);

                w2=a/3+(b*1*i-hd*1)/3+2*v/3+2*l*hm*(i-(j-1)*n)/3;

                w3=a+i*b*l+b*l;

                if w1<=0

            w11=0;

        elseif (w1>=a+i*b*l-b*l) 

            w11=a+i*b*l-b*l;

        else w11=w1;

        end

                q1=l*(a/(2*c)+b*l/(4*c)*(2*i-l)-w11/(2*c)-hd*l/(4*c));

                lr21=w11*q1-v*q1-l*hm*(i-(j-1)*n)*q1;

                lr11=((a-w11+hd*(i-l)*l)*(a-w11+hd*(i-l)*l)*l+(b-hd)*(a-w11+hd*(i-1)*l)*l*l*(2*i-1)+((b-hd)*(b-hd)*l*l*l*(3*i*i-3*i+1))/3)/(4*c);

                if(w2<a+i*b*l-b*l)

                        w22=a+i*b*l-b*l;

        elseif(w2>=a+i*b*l-hd*l)

              w22=a+i*b*l-hd*l;

                else w22=w2;

        end

                q2=((a+b*l*i-hd*l-w22)*(a+b*l*i-hd*l-w22))/(4*c*(b-hd));

                lr22=w22*q2-v*q2-l*hm*(i-(j-1)*n)*q2;

                lr12=((a-w22+hd*(i-1)*l)*(a-w22+hd*(i-1)*l)*(a-w22+b*i*l-hd*l)/(b-hd)+(b-hd)*(-w22+hd*(i-1)*l)*(a-w22+b*i*l-hd*l)/(b-hd)*(w22+b*i*l+hd*l-a-2*hd*i*l)/(b-hd)+(b-hd)*(b-hd)*(a-w22+b*i*l-hd*l)/(b-hd)*(((w22-a-hd*l*(i-1))*(w22-a-hd*l*(i-l)))/((b-hd)*(b-hd))+(w22-a-hd*l*(i-1))/(b-hd)*i*l+i*i*l*l)/3)/(4*c);

                if(lr21>=lr22)

            sprintf('i=j=%3f,w=%3f, q=%3f, lr2=%3f, lr1=3%f',i,w11,q1,lr21,lr11)

        else

            sprintf('i=j=%3f,w=%3f, q=%3f, lr2=%3f, lr1=3%f',i,w11,q1,lr21,lr11)        

        end                

end

贊一下(1人)

回復此樓

» 本帖已獲得的紅花（最新10朵）

wcwl888

showmethemoney

4樓2015-03-20 13:51:43

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

wcwl888

金蟲 (著名寫手)

中華人民共和國第一任航母艦長

應助: 12 (小學生)
金幣: 2639.4
散金: 210
紅花: 9
沙發(fā): 2
帖子: 1217
在線: 169.2小時
蟲號: 843785
注冊: 2009-09-08
性別: GG
專業(yè): 信息系統(tǒng)與管理

送紅花一朵

引用回帖:

4樓: Originally posted by csgt0 at 2015-03-20 13:51:43
公式只能你自己對了

clc
clear

a=40;
b=1;
c=0.25;
l=10;
v=4;
hd=0.8;
hm=0.8;
n=1;

for i=1:10
j=i;
w1=0.5*a+0.25*b*1*(2*i-1)+0.5*v-0.25*hd*1+0.5*1*hm*(i-(j-l)*n);
w2=a/3 ...

大牛，您好，我是不是只要把這個代碼放到matlab命令窗口運行出文件再輸參數求結果啊？
這個代碼的格式是matlab里能夠識別的腳本語言嗎？跟我以前看到的不太一樣（不是懷疑，也不敢，只是不懂，向大牛請教以下，望見諒哈）

贊一下

回復此樓

佛說：說了就是沒說，沒說就是說了，施主你又何必問呢！阿彌陀佛～

5樓2015-03-20 21:01:34

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

wcwl888

金蟲 (著名寫手)

中華人民共和國第一任航母艦長

應助: 12 (小學生)
金幣: 2639.4
散金: 210
紅花: 9
沙發(fā): 2
帖子: 1217
在線: 169.2小時
蟲號: 843785
注冊: 2009-09-08
性別: GG
專業(yè): 信息系統(tǒng)與管理

引用回帖:

大神，我在matlab里面運行時出現如下錯誤，您看下哪里出現錯誤了？拜謝了。
??? lr11=((a-w11+hd*(i-1)*l)*(a-w11+hd*(i-1)*l*l+(b-hd)*(a-w11+hd*(i-1)*l)*l*l*(2*i-

|
Error: Incomplete or misformed expression or statement.

>>

贊一下

回復此樓

佛說：說了就是沒說，沒說就是說了，施主你又何必問呢！阿彌陀佛～

6樓2015-03-20 22:19:07

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

wcwl888

金蟲 (著名寫手)

中華人民共和國第一任航母艦長

應助: 12 (小學生)
金幣: 2639.4
散金: 210
紅花: 9
沙發(fā): 2
帖子: 1217
在線: 169.2小時
蟲號: 843785
注冊: 2009-09-08
性別: GG
專業(yè): 信息系統(tǒng)與管理

引用回帖:

6樓: Originally posted by wcwl888 at 2015-03-20 22:19:07
大神，我在matlab里面運行時出現如下錯誤，您看下哪里出現錯誤了？拜謝了。
??? lr11=((a-w11+hd*(i-1)*l)*(a-w11+hd*(i-1)*l*l+(b-hd)*(a-w11+hd*(i-1)*l)*l*l*(2*i-

...

大牛您好，不知道您是不是上班比較忙，小弟這幾天一直在刷版，想看您的更新，這個問題一直沒解決，不知道您是否能抽空指教以下，謝謝。

贊一下

回復此樓

佛說：說了就是沒說，沒說就是說了，施主你又何必問呢！阿彌陀佛～

7樓2015-03-21 11:06:08

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經驗: +2
計算強帖: 2
應助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準晶和納米晶
管轄: 計算模擬

引用回帖:

7樓: Originally posted by wcwl888 at 2015-03-21 11:06:08
大牛您好，不知道您是不是上班比較忙，小弟這幾天一直在刷版，想看您的更新，這個問題一直沒解決，不知道您是否能抽空指教以下，謝謝。...

在matlab里有個file-new-script,然后把代碼復制進去，保存后運行（按F5運行）

贊一下

回復此樓

showmethemoney

8樓2015-03-23 17:59:20

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關版塊跳轉我要訂閱樓主 wcwl888 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 286求調劑 +10	Faune 2026-03-21	10/500	2026-03-21 23:34 by 314126402
[考研] 一志愿北京化工大學070300 學碩336求調劑 +3	vv迷 2026-03-21	4/200	2026-03-21 23:05 by f19980501
[考研] 一志愿西安交通大學材料工程專業(yè) 282分求調劑 +10	楓橋ZL 2026-03-18	12/600	2026-03-21 22:02 by peike
[考研] 一志愿深大，0703化學，總分302，求調劑 +4	七月-七七 2026-03-21	4/200	2026-03-21 18:20 by 學員8dgXkO
[考研] 材料學碩333求調劑 +3	北道巷 2026-03-18	3/150	2026-03-21 18:17 by 學員8dgXkO
[考研] 313求調劑 +4	肆叁貳壹22 2026-03-19	4/200	2026-03-21 17:33 by ColorlessPI
[考研] 277材料科學與工程080500求調劑 +6	自由煎餅果子 2026-03-16	6/300	2026-03-21 17:21 by 學員8dgXkO
[考研] 265求調劑 +3	Jack?k?y 2026-03-17	3/150	2026-03-21 03:17 by JourneyLucky
[考研] 083200學碩321分一志愿暨南大學求調劑 +3	innocenceF 2026-03-17	3/150	2026-03-21 02:35 by JourneyLucky
[考研] 二本跨考鄭大材料306英一數二 +3	z1z2z3879 2026-03-17	3/150	2026-03-21 02:29 by JourneyLucky
[考研] 一志愿華中科技大學，080502，354分求調劑 +5	守候夕陽CF 2026-03-18	5/250	2026-03-21 01:06 by JourneyLucky
[考研] 考研調劑求學校推薦 +3	伯樂29 2026-03-18	5/250	2026-03-20 22:59 by JourneyLucky
[考研] 323求調劑 +3	洼小桶 2026-03-18	3/150	2026-03-20 22:54 by JourneyLucky
[考研] 290求調劑 +7	^O^乜 2026-03-19	7/350	2026-03-20 21:43 by JourneyLucky
[考研] 一志愿武理材料工程348求調劑 +3	￣^￣゜汗 2026-03-19	4/200	2026-03-20 21:01 by zhukairuo
[考博] 招收博士1-2人 +3	QGZDSYS 2026-03-18	3/150	2026-03-20 11:58 by 呱呱呱呱叫
[考研] 081700化工學碩調劑 +3	【1】 2026-03-16	3/150	2026-03-19 23:40 by edmund7
[考博] 26博士申請 +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考博] 26申博 +4	八6八68 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 一志愿，福州大學材料專碩339分求調劑 +3	木子momo青爭 2026-03-15	3/150	2026-03-17 07:52 by laoshidan

24小時熱門版塊排行榜

[求助] 菜鳥求助代碼運用 已有1人參與

» 猜你喜歡

» 本主題相關價值貼推薦，對您同樣有幫助:

【答案】應助回帖

【答案】應助回帖

» 本帖已獲得的紅花（最新10朵）

[求助] 菜鳥求助代碼運用已有1人參與