版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

wcwl888

金蟲 (著名寫手)

中華人民共和國(guó)第一任航母艦長(zhǎng)

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

[求助] 菜鳥求助代碼運(yùn)用已有1人參與

親愛的大神們：
各位好，學(xué)弟近期寫論文建模需要用到編程求解，但不懂運(yùn)用何種軟件，最近看到一篇文獻(xiàn)里使用的模型和我的比較相似，苦于看不懂不知道該用什么軟件也不知道如何運(yùn)行，故請(qǐng)教大神指點(diǎn)迷津，下面是C++源代碼：（附件里是模型需要的參數(shù)數(shù)值和運(yùn)行結(jié)果，懂的大神可以參考一下）
#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<<"請(qǐng)依次輸入?yún)?shù)值:";
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;
}

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

參數(shù).png

計(jì)算結(jié)果.png

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有6人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有203人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有68人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

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

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

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

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

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

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

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

你把源碼疏導(dǎo)c里面就可以了，注意下格式，你這直接從pdf轉(zhuǎn)的肯定錯(cuò)誤多

贊一下

回復(fù)此樓

showmethemoney

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

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

wcwl888

金蟲 (著名寫手)

中華人民共和國(guó)第一任航母艦長(zhǎng)

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

引用回帖:

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

您能把這個(gè)代碼轉(zhuǎn)換成MATLAB語言嗎？只要轉(zhuǎn)換成matlab里可以運(yùn)行的就可以了，金幣全給您！

贊一下

回復(fù)此樓

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

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

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

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

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

【答案】應(yīng)助回帖

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

引用回帖:

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

...

公式只能你自己對(duì)了

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人)

回復(fù)此樓

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

wcwl888

showmethemoney

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

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

wcwl888

金蟲 (著名寫手)

中華人民共和國(guó)第一任航母艦長(zhǎng)

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

送紅花一朵

引用回帖:

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

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 ...

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

贊一下

回復(fù)此樓

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

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

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

wcwl888

金蟲 (著名寫手)

中華人民共和國(guó)第一任航母艦長(zhǎng)

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

引用回帖:

大神，我在matlab里面運(yùn)行時(shí)出現(xiàn)如下錯(cuò)誤，您看下哪里出現(xiàn)錯(cuò)誤了？拜謝了。
??? 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.

>>

贊一下

回復(fù)此樓

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

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

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

wcwl888

金蟲 (著名寫手)

中華人民共和國(guó)第一任航母艦長(zhǎng)

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

引用回帖:

6樓: Originally posted by wcwl888 at 2015-03-20 22:19:07
大神，我在matlab里面運(yùn)行時(shí)出現(xiàn)如下錯(cuò)誤，您看下哪里出現(xiàn)錯(cuò)誤了？拜謝了。
??? 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-

...

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

贊一下

回復(fù)此樓

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

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

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

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

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

引用回帖:

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

在matlab里有個(gè)file-new-script,然后把代碼復(fù)制進(jìn)去，保存后運(yùn)行（按F5運(yùn)行）

贊一下

回復(fù)此樓

showmethemoney

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

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 wcwl888 的主題更新

返回列表

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

[考研] 336求調(diào)劑 +3	rmc8866 2026-03-21	3/150	2026-03-21 11:16 by 求調(diào)劑zz
[考研] 一志愿華中科技大學(xué)071000，求調(diào)劑 +3	沿岸有貝殼6 2026-03-21	3/150	2026-03-21 10:35 by 暮云清寒
[考研] 316求調(diào)劑 +6	梁茜雯 2026-03-19	6/300	2026-03-21 06:32 by Ecowxq666！
[考研] 08工科 320總分求調(diào)劑 +6	梨花珞晚風(fēng) 2026-03-17	6/300	2026-03-21 03:40 by JourneyLucky
[考研] 301求調(diào)劑 +10	yy要上岸呀 2026-03-17	10/500	2026-03-21 03:14 by JourneyLucky
[考研] 初始318分求調(diào)劑（有工作經(jīng)驗(yàn)） +3	1911236844 2026-03-17	3/150	2026-03-21 02:33 by JourneyLucky
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3	0703武芊慧雪304 2026-03-18	3/150	2026-03-21 00:48 by JourneyLucky
[考研] 311求調(diào)劑 +5	冬十三 2026-03-18	5/250	2026-03-21 00:16 by JourneyLucky
[考研] 294求調(diào)劑材料與化工專碩 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 by JourneyLucky
[考研] 一志愿西南交通專碩材料355 本科雙非求調(diào)劑 +5	西南交通專材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 289求調(diào)劑 +6	懷瑾握瑜l 2026-03-20	6/300	2026-03-20 20:30 by 學(xué)員8dgXkO
[考研] 0703化學(xué)調(diào)劑，六級(jí)已過，有科研經(jīng)歷 +13	曦熙兮 2026-03-15	13/650	2026-03-20 19:35 by Dream007008
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3	葵梓衛(wèi)隊(duì) 2026-03-20	3/150	2026-03-20 19:28 by zhukairuo
[考研] 求調(diào)劑 +3	暗涌afhb 2026-03-16	3/150	2026-03-20 00:28 by 河南大學(xué)校友
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國(guó)級(jí)證書40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 293求調(diào)劑 +11	zjl的號(hào) 2026-03-16	16/800	2026-03-18 08:10 by zhukairuo
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu

24小時(shí)熱門版塊排行榜

[求助] 菜鳥求助代碼運(yùn)用 已有1人參與

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

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

[求助] 菜鳥求助代碼運(yùn)用已有1人參與

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助: