 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 937 | 回復(fù): 4 | ||
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ù)此樓» 猜你喜歡
311求調(diào)劑
已經(jīng)有4人回復(fù)
-
304求調(diào)劑
已經(jīng)有4人回復(fù)
-
354求調(diào)劑
已經(jīng)有7人回復(fù)
-
321求調(diào)劑
已經(jīng)有3人回復(fù)
-
274求調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有15人回復(fù)
-
085600材料與化工調(diào)劑 324分
已經(jīng)有8人回復(fù)
-
能源材料化學(xué)課題組招收碩士研究生8-10名
已經(jīng)有10人回復(fù)
-
267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑
已經(jīng)有8人回復(fù)
-
328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷
已經(jīng)有10人回復(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ù)
- 急用！�。⌒孪x求助有關(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
- 注冊(cè): 2015-02-21
- 性別: GG
- 專業(yè): 流體力學(xué)
4215623392樓2015-08-31 13:44:10
421562339
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1122.3
- 紅花: 3
- 帖子: 94
- 在線: 27.5小時(shí)
- 蟲號(hào): 1324594
- 注冊(cè): 2011-06-16
- 性別: GG
- 專業(yè): 光學(xué)
3樓2015-09-01 10:07:59
happy123xp
銀蟲 (小有名氣)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 581.6
- 紅花: 2
- 帖子: 63
- 在線: 49.5小時(shí)
- 蟲號(hào): 3693097
- 注冊(cè): 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)這樣的問題，樓主可以參考參考我說的，謝謝。 |
» 本帖已獲得的紅花（最新10朵）
4樓2015-09-01 16:56:51
421562339
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1122.3
- 紅花: 3
- 帖子: 94
- 在線: 27.5小時(shí)
- 蟲號(hào): 1324594
- 注冊(cè): 2011-06-16
- 性別: GG
- 專業(yè): 光學(xué)
5樓2015-09-02 10:14:05
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 421562339 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +4 | 脫穎而出 2026-03-16 | 10/500 |
|
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +8 | SUICHILD 2026-03-12 | 8/400 |
|
|
|
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +7 | Liwangman 2026-03-15 | 7/350 |
|
|
|
[考研]
 271材料工程求調(diào)劑
+4
|
.6lL 2026-03-18 | 4/200 |
|
|
|
[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +4 | 1孫悟空 2026-03-17 | 4/200 |
|
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 9/450 |
|
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號(hào) 2026-03-16 | 16/800 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 【0856】化學(xué)工程（085602）313 分，本科學(xué)科評(píng)估A類院校化學(xué)工程與工藝，誠(chéng)求調(diào)劑 +7 | 小劉快快上岸 2026-03-11 | 8/400 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
|
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|