 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 938 | 回復(fù): 4 | ||
421562339新蟲 (小有名氣)
|
[求助]
有關(guān)MATLAB數(shù)學(xué)計算編程問題的求助 已有1人參與
|
|
各位蟲友，能不能用MATLAB幫我計算一下，這些公式是根據(jù)Kliein & Theilacker的傾斜面輻射計算公式（見圖片1），我們編寫的程序內(nèi)循環(huán)為計算傾斜面1到12月的輻射，然后求年總輻射；外循環(huán)為傾斜面角度，0到90度，程序編寫的目的是為了求傾斜面的最佳傾角，即傾角為多少時傾斜面接受的年輻射最大，我們編寫的程序循環(huán)計算有問題，每次運(yùn)行出來的結(jié)果與我們輸入的緯度及月平均水平面總輻射無關(guān)，最佳傾角總是為90度減去0度，如果外循環(huán)傾角的取值范圍發(fā)生變化，計算結(jié)果總為末值減去初始值，不是我們想要的結(jié)果，應(yīng)該是程序出問題了，以下是我們編寫的程序，望各位蟲友指教。 %傾斜面月平均太陽輻照量計算(klien&Theilacker) %r---a對于朝向赤道的傾斜面上，在北半球朝南的傾斜面上，其月平均太陽總輻照量與水平面上月平均總輻照量之比 %輸入月平均輻射總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ù)此樓» 猜你喜歡
271材料工程求調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有16人回復(fù)
-
304求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料工程專碩調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有3人回復(fù)
-
求材料調(diào)劑
已經(jīng)有8人回復(fù)
-
294求調(diào)劑材料與化工專碩
已經(jīng)有5人回復(fù)
-
一志愿華中科技大學(xué)，080502，354分求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 怎樣編寫別人論文里面的有關(guān)非線性PCA的MATLAB程序
已經(jīng)有0人回復(fù)
- MATLAB軟件與基礎(chǔ)數(shù)學(xué)實驗
已經(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è)計方法(第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è)計用matlab編程
已經(jīng)有1人回復(fù)
- 【分享】計算數(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小時
- 蟲號: 3693097
- 注冊: 2015-02-21
- 性別: GG
- 專業(yè): 流體力學(xué)
4215623392樓2015-08-31 13:44:10
421562339
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1122.3
- 紅花: 3
- 帖子: 94
- 在線: 27.5小時
- 蟲號: 1324594
- 注冊: 2011-06-16
- 性別: GG
- 專業(yè): 光學(xué)
3樓2015-09-01 10:07:59
happy123xp
銀蟲 (小有名氣)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 581.6
- 紅花: 2
- 帖子: 63
- 在線: 49.5小時
- 蟲號: 3693097
- 注冊: 2015-02-21
- 性別: GG
- 專業(yè): 流體力學(xué)
【答案】應(yīng)助回帖
|
回去跑了一下樓主的程序，并沒有解決樓主的疑問，但是卻發(fā)現(xiàn)了幾個問題。 1.當(dāng)輸入緯度為20時，計算錯誤，出現(xiàn)了sss是復(fù)數(shù)的情況，是sss的反三角函數(shù)里面的計算值超過了1，懷疑c(i)計算式有問題。 2.無論怎么變化緯度，居然總的最大輻射量沒怎變，查后發(fā)現(xiàn)r計算式有問題，應(yīng)把dd*cosd(a)都括起來，不然計算順序有問題。 3.最后幾句話也有問題，max函數(shù)里面的矩陣不要j，應(yīng)為sum，你想想，你的j最后值肯定是181，max(sum(j)得到的永遠(yuǎn)是sum里面最后一個，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)致，最后一個總是最大的，所以找到的cols總是181，總是90度。 最后，我不了解樓主研究領(lǐng)域，但是總是覺得樓主某個理論環(huán)節(jié)有問題，才會出現(xiàn)這樣的問題，樓主可以參考參考我說的，謝謝。 |
» 本帖已獲得的紅花（最新10朵）
4樓2015-09-01 16:56:51
421562339
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1122.3
- 紅花: 3
- 帖子: 94
- 在線: 27.5小時
- 蟲號: 1324594
- 注冊: 2011-06-16
- 性別: GG
- 專業(yè): 光學(xué)
5樓2015-09-02 10:14:05
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 421562339 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩英一數(shù)二306 +4 | z1z2z3879 2026-03-18 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
|
[考研] 290求調(diào)劑 +6 | 孔志浩 2026-03-12 | 11/550 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy！ 2026-03-13 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
|
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
|
[考研] 求調(diào)劑 +5 | 一定有學(xué)上- 2026-03-12 | 5/250 |
|
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|