 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 946 | 回復(fù): 4 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
421562339新蟲 (小有名氣)
|
[求助]
有關(guān)MATLAB數(shù)學(xué)計(jì)算編程問題的求助 已有1人參與
|
|
|
各位蟲友，能不能用MATLAB幫我計(jì)算一下，這些公式是根據(jù)Kliein & Theilacker的傾斜面輻射計(jì)算公式（見圖片1），我們編寫的程序內(nèi)循環(huán)為計(jì)算傾斜面1到12月的輻射，然后求年總輻射；外循環(huán)為傾斜面角度，0到90度，程序編寫的目的是為了求傾斜面的最佳傾角，即傾角為多少時(shí)傾斜面接受的年輻射最大，我們編寫的程序循環(huán)計(jì)算有問題，每次運(yùn)行出來的結(jié)果與我們輸入的緯度及月平均水平面總輻射無關(guān)，最佳傾角總是為90度減去0度，如果外循環(huán)傾角的取值范圍發(fā)生變化，計(jì)算結(jié)果總為末值減去初始值，不是我們想要的結(jié)果，應(yīng)該是程序出問題了，以下是我們編寫的程序，望各位蟲友指教。 %傾斜面月平均太陽輻照量計(jì)算(klien&Theilacker) %r---a對(duì)于朝向赤道的傾斜面上，在北半球朝南的傾斜面上，其月平均太陽總輻照量與水平面上月平均總輻照量之比 %輸入月平均輻射總h(i)，散射輻射d(i),和直接輻射b(i) %輸入當(dāng)?shù)鼐暥萢 %c(i)----第i月的太陽赤緯角 %c=23.45*sin[360/365*(284+n)](n--為一年中叢元旦算起的天數(shù)) a=input('輸入當(dāng)?shù)鼐暥?度):') j=0; sum_i=0; for t=0:0.5:90 j=j+1 for i=1:1:12 h=[237.2 310.3 429.5 541.5 628.2 621.1 625.7 578.3 436.0 352.0 263.9 212]; d=[149.0 178.8 258.1 304.7 324.1 299.7 271.6 240.2 213.6 181.8 148.9 133.7]; b=h-d n=[17,47,75,105,135,162,198,228,258,288,318,344]; c(i)=23.45*sind(360/365*(284+n(i))); s=acosd(-tand(a)*tand(c(i))); ss=min(s,acosd(-tand(a-t)*tand(c(i)))); sss=acosd(-tand(a-t)*tand(c(i))); dd=sind(s)-pi/180*s*cosd(s); aa=0.409+0.5016*sind(s-60); bb=0.6609-0.4767*sind(s-60); r=cosd(a-t)/dd*cosd(a)*((aa-d(i)/h(i))*(sind(ss)-pi/180*ss*cosd(sss))+bb/2*(pi/180*ss+sind(ss)*(cosd(ss)-2*cosd(sss))))+d(i)/2*h(i)*(1+cosd(t))+0.1*(1-cosd(t)); ht=r*h(i); sum_i=ht+sum_i; end sum(j)=sum_i; fprintf('sum(j)=%d\n',sum(j)) end hy=max(sum(j)); find(j==max(sum(j))); tt=(j-1)*0.5; fprintf('hy=%d\n,tt=%d\n',hy,tt) {P39N7KG@Q$KPWM$H0MIX1V.png |
回復(fù)此樓» 猜你喜歡
材料科學(xué)與工程求調(diào)劑
已經(jīng)有3人回復(fù)
-
289求調(diào)劑
已經(jīng)有5人回復(fù)
-
329求調(diào)劑，一志愿西北工業(yè)大學(xué)，材料工程（085601）
已經(jīng)有6人回復(fù)
-
化學(xué)工程085602 305分求調(diào)劑
已經(jīng)有21人回復(fù)
-
26申博自薦
已經(jīng)有5人回復(fù)
-
296求調(diào)劑
已經(jīng)有5人回復(fù)
-
總分293求調(diào)劑
已經(jīng)有10人回復(fù)
-
材料專碩調(diào)劑
已經(jīng)有3人回復(fù)
-
調(diào)劑310
已經(jīng)有12人回復(fù)
-
化學(xué)0703 調(diào)劑 306分 一志愿211
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 怎樣編寫別人論文里面的有關(guān)非線性PCA的MATLAB程序
已經(jīng)有0人回復(fù)
- MATLAB軟件與基礎(chǔ)數(shù)學(xué)實(shí)驗(yàn)
已經(jīng)有10人回復(fù)
- matlab中文幫助.CHM
已經(jīng)有52人回復(fù)
- matlab 非線性微分方程求解
已經(jīng)有3人回復(fù)
- 分享——與MATLAB有關(guān)的工程數(shù)學(xué)（第二版，英文版）
已經(jīng)有120人回復(fù)
- 《機(jī)械工程控制基礎(chǔ)》【轉(zhuǎn)載】
已經(jīng)有128人回復(fù)
- 普通高等教育十二五規(guī)劃教材—《現(xiàn)代機(jī)械設(shè)計(jì)方法(第2版)》【轉(zhuǎn)載】
已經(jīng)有184人回復(fù)
- 《MATLAB--小波變換教程》
已經(jīng)有481人回復(fù)
- MathWorks 宣布推出 2011b 版 MATLAB 和 Simulink
已經(jīng)有5人回復(fù)
- 急用�。�！新蟲求助有關(guān)階梯軸的優(yōu)化設(shè)計(jì)用matlab編程
已經(jīng)有1人回復(fù)
- 【分享】計(jì)算數(shù)學(xué)電子書下載
已經(jīng)有23人回復(fù)
- 【已解決】數(shù)學(xué)建模的資源
已經(jīng)有11人回復(fù)
- 【資源】COMSOL Multiphysics 有關(guān)于化學(xué)工程和生物醫(yī)學(xué)工程方面的論文
已經(jīng)有22人回復(fù)
happy123xp
銀蟲 (小有名氣)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 581.6
- 紅花: 2
- 帖子: 63
- 在線: 49.5小時(shí)
- 蟲號(hào): 3693097
- 注冊: 2015-02-21
- 性別: GG
- 專業(yè): 流體力學(xué)
【答案】應(yīng)助回帖
|
回去跑了一下樓主的程序，并沒有解決樓主的疑問，但是卻發(fā)現(xiàn)了幾個(gè)問題。 1.當(dāng)輸入緯度為20時(shí)，計(jì)算錯(cuò)誤，出現(xiàn)了sss是復(fù)數(shù)的情況，是sss的反三角函數(shù)里面的計(jì)算值超過了1，懷疑c(i)計(jì)算式有問題。 2.無論怎么變化緯度，居然總的最大輻射量沒怎變，查后發(fā)現(xiàn)r計(jì)算式有問題，應(yīng)把dd*cosd(a)都括起來，不然計(jì)算順序有問題。 3.最后幾句話也有問題，max函數(shù)里面的矩陣不要j，應(yīng)為sum，你想想，你的j最后值肯定是181，max(sum(j)得到的永遠(yuǎn)是sum里面最后一個(gè)，find沒起到作用，應(yīng)改為 hy=max(sum); [rows,cols]=find(hy==sum);tt=(cols-1)*0.5;。 4.羅里吧嗦一大堆，最后還是發(fā)現(xiàn)問題還是沒解決，程序應(yīng)該是沒有問題了，但是得到的sum矩陣?yán)锩娴脑乜偸沁f增的，所以導(dǎo)致，最后一個(gè)總是最大的，所以找到的cols總是181，總是90度。 最后，我不了解樓主研究領(lǐng)域，但是總是覺得樓主某個(gè)理論環(huán)節(jié)有問題，才會(huì)出現(xiàn)這樣的問題，樓主可以參考參考我說的，謝謝。 |
4215623394樓2015-09-01 16:56:51
happy123xp
銀蟲 (小有名氣)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 581.6
- 紅花: 2
- 帖子: 63
- 在線: 49.5小時(shí)
- 蟲號(hào): 3693097
- 注冊: 2015-02-21
- 性別: GG
- 專業(yè): 流體力學(xué)
2樓2015-08-31 13:44:10
421562339
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1122.3
- 紅花: 3
- 帖子: 94
- 在線: 27.5小時(shí)
- 蟲號(hào): 1324594
- 注冊: 2011-06-16
- 性別: GG
- 專業(yè): 光學(xué)
3樓2015-09-01 10:07:59
421562339
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1122.3
- 紅花: 3
- 帖子: 94
- 在線: 27.5小時(shí)
- 蟲號(hào): 1324594
- 注冊: 2011-06-16
- 性別: GG
- 專業(yè): 光學(xué)
5樓2015-09-02 10:14:05
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300化學(xué)354求調(diào)劑 +6 | 101次希望 2026-03-28 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研]
 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
+10
|
nnnnnnn5 2026-03-25 | 10/500 |
|
|
[考研] 0856求調(diào)劑 +7 | 楒桉 2026-03-28 | 7/350 |
|
|
|
[考研] 調(diào)劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
|
[考研] 289求調(diào)劑 +13 | 新時(shí)代材料 2026-03-27 | 13/650 |
|
|
|
[考研] 一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑 +7 | Micky11223 2026-03-25 | 7/350 |
|
|
|
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +5 | niko- 2026-03-27 | 5/250 |
|
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分，前三科223分 +5 | dongfang59 2026-03-27 | 5/250 |
|
|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩，求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 12/600 |
|
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
|
[考研] 各位老師您好：本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|
|
|
[考研] 284求調(diào)劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|