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

回復此樓

» 猜你喜歡

中德博士后2026項目5月30日截止- 國內博士后入站、德國兩年全職已經有0人回復
面向全球申請人的中德博士后2026項目5月30日截止- 國內博后入站德國兩年全職已經有0人回復
物理學I論文潤色/翻譯怎么收費? 已經有156人回復
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經有2人回復
重慶交大材料學院26年碩士招生復試實施細則已出。歡迎加入光子學微結構與器件課題組。已經有0人回復
北工大宋曉艷教授團隊招收2026年9月入學博士【學博和專博名額均有】！已經有1人回復
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收10人已經有0人回復
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收調劑，10人已經有0人回復
計算與實驗科學工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學換能器分會S8 已經有0人回復
0702一志愿吉大B區(qū)求調劑有論文已經有2人回復

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

用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的回帖

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的回帖

查看全部 8 個回答

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的回帖

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

[考研] 考研調劑 +5	Sanmu-124 2026-03-26	5/250	2026-03-30 13:36 by chemdavid
[考研] 本科211生物醫(yī)學工程085409求調劑339分 +3	里子木yy 2026-03-29	3/150	2026-03-30 13:29 by gyzj2026
[考研] 生物技術與工程 +7	1294608413 2026-03-25	8/400	2026-03-30 11:36 by 唐沐兒
[考研] 南京大學化學調劑 +10	景隨風 2026-03-29	15/750	2026-03-30 11:21 by limeifeng
[考研] 332求調劑 +6	Lyy930824@ 2026-03-29	6/300	2026-03-30 10:53 by 1172367218
[考研] 一志愿河北工業(yè)大學0817化工278分求調劑 +14	jhybd 2026-03-23	19/950	2026-03-30 10:34 by vgtyfty
[考研] 考研調劑 +7	小蠟新筆 2026-03-29	7/350	2026-03-29 19:00 by 學員8dgXkO
[考研] 一志愿：西北大學，英一數一408-284分求調劑 +4	12.27 2026-03-27	4/200	2026-03-29 14:40 by zhshch
[考研] 調劑求院校招收 +6	鶴鯨鴿 2026-03-28	6/300	2026-03-29 08:15 by fmesaito
[考研] 一志愿211 初試270分求調劑 +8	谷雨上岸 2026-03-23	9/450	2026-03-29 01:00 by 我是小康
[考研] 一志愿太原理工安全工程300分，求調劑 +5	0857求調劑. 2026-03-24	6/300	2026-03-28 22:04 by zhq0425
[考研] 求調劑 +7	爭取九點睡 2026-03-28	8/400	2026-03-28 21:07 by 爭取九點睡
[考研] 315求調劑 +4	akie... 2026-03-28	5/250	2026-03-28 21:05 by zhq0425
[考研] 一志愿廈門大學化學學碩307求調劑 +10	y7czhao 2026-03-26	10/500	2026-03-28 14:23 by 唐沐兒
[考研] 331環(huán)境科學與工程求調劑 +3	熠然好運氣 2026-03-27	3/150	2026-03-28 04:11 by fmesaito
[考研] 0703化學338求調劑！ +6	Zuhui0306 2026-03-26	7/350	2026-03-27 10:35 by shangxh
[考研] 材料科學與工程 317求調劑 +4	JKSOIID 2026-03-26	4/200	2026-03-26 15:58 by 不吃魚的貓
[考研] 生物學學碩求調劑 +7	小羊睡著了? 2026-03-23	10/500	2026-03-25 02:24 by 清風拂揚。 m
[考研] 材料專碩找調劑 +5	哈哈哈吼吼吼哈 2026-03-23	5/250	2026-03-24 19:07 by 了了了了。。
[考研] 292求調劑 +4	鵝鵝鵝額額額額?/a> 2026-03-24	4/200	2026-03-24 16:41 by peike