 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 944 | 回復(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ù)此樓» 猜你喜歡
考研調(diào)劑
已經(jīng)有3人回復(fù)
-
317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑
已經(jīng)有9人回復(fù)
-
297求調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
已經(jīng)有10人回復(fù)
-
343求調(diào)劑
已經(jīng)有5人回復(fù)
-
085601材料工程找調(diào)劑
已經(jīng)有7人回復(fù)
-
求調(diào)劑，一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分
已經(jīng)有6人回復(fù)
-
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷
已經(jīng)有5人回復(fù)
-
340求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿北化求調(diào)劑
已經(jīng)有4人回復(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ù)
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é)
2樓2015-08-31 13:44:10
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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑，一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +6 | @taotao 2026-03-29 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿北化求調(diào)劑 +4 | Jsman 2026-03-22 | 4/200 |
|
|
|
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥 求調(diào)劑 +3 | 篤行致遠(yuǎn). 2026-03-28 | 3/150 |
|
|
[考研] 291求調(diào)劑 +4 | Y-cap 2026-03-29 | 4/200 |
|
|
|
[考研] 一志愿華東師范大學(xué)有機(jī)化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑! +3 | 真名有冰 2026-03-29 | 4/200 |
|
|
|
[考研] 0703化學(xué) +11 | 妮妮ninicgb 2026-03-27 | 11/550 |
|
|
|
[考研] 生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績(jī)338 +6 | YYYYYNNNNN 2026-03-26 | 7/350 |
|
|
|
[考研] 071000生物學(xué)求調(diào)劑，初試成績(jī)343 +7 | 小小甜面團(tuán) 2026-03-25 | 7/350 |
|
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|
|
|
[考研] 291求調(diào)劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
|
[考研] 311求調(diào)劑 +3 | 希望上岸阿小楊 2026-03-23 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +4 | 零八# 2026-03-27 | 4/200 |
|
|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷 +10 | 材料學(xué)257求調(diào)劑 2026-03-23 | 11/550 |
|
|
|
[考研] 359求調(diào)劑 +4 | 王了個(gè)楠 2026-03-25 | 4/200 |
|
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4 | 小魚愛有機(jī) 2026-03-25 | 4/200 |
|
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
|
[考研] 284求調(diào)劑 +15 | Zhao anqi 2026-03-22 | 15/750 |
|
|
|
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生（材料、化學(xué)、電化學(xué)，環(huán)境） +4 | 我愛學(xué)電池 2026-03-23 | 4/200 |
|
|
|
[考研] 300分，材料，求調(diào)劑，英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|