matlab做数学建模，有这么段代码，这有错么？

[复制链接]

字体大小: 正常放大

杨利霞

5273 主题	82 听众	17万积分

TA的每日心情

	开心 2021-8-11 17:59

签到天数: 17 天

[LV.4]偶尔看看III

网络挑战赛参赛者

自我介绍: 本人女，毕业于内蒙古科技大学，担任文职专业，毕业专业英语。

群组: 2018美赛大象算法课程

群组: 2018美赛护航培训课程

群组: 2019年数学中国站长建

群组: 2019年数据分析师课程

群组: 2018年大象老师国赛优

电梯直达

1^#

发表于 2019-7-22 14:44 |只看该作者 |正序浏览

|招呼Ta 关注Ta

matlab做数学建模，有这么段代码，这有错么？

结果实在是跑不对
selftimeform = [14.7:0.05:15.7];
x=[1.03650000000000,1.06990000000000,1.10380000000000,1.13830000000000,1.17320000000000,1.20870000000000,1.24480000000000,1.28150000000000,1.31890000000000,1.35680000000000,1.39550000000000,1.43490000000000,1.47510000000000,1.51600000000000,1.55770000000000,1.60030000000000,1.64380000000000,1.68820000000000,1.73370000000000,1.78010000000000,1.82770000000000]
y=[0.497300000000000,0.502900000000000,0.508500000000000,0.514200000000000,0.519800000000000,0.525500000000000,0.531100000000000,0.536800000000000,0.542600000000000,0.548300000000000,0.554100000000000,0.559800000000000,0.565700000000000,0.571500000000000,0.577400000000000,0.583300000000000,0.589200000000000,0.595200000000000,0.601300000000000,0.607400000000000,0.613500000000000]
yingchang = sqrt(x.^2+y.^2)
N=107; %N为15年4月18日与1月1日的日期之差
derta=23.45*sin(360*(284+N)/365)*pi/180; %赤纬角
omiga=15*(selftimeform-(120-111.02)/15-12)*pi/180; %计算时角
for j=1:21
for fai=1:90
fenzi = sin(derta)*(sin(fai*pi/180)*sin(fai*pi/180))+cos(derta)*cos(fai*pi/180)*sin(fai*pi/180)*cos(omiga(j))-sin(derta)
fenmu = sqrt(1-(sin(derta)*sin(fai*pi/180)+cos(derta)*cos(omiga(j))*cos(fai*pi/180))*(sin(derta)*sin(fai*pi/180)+cos(derta)*cos(omiga(j))*cos(fai*pi/180)))*cos(fai*pi/180)
f(j,fai) = fenzi/fenmu - x(j)/yingchang(j)
end
end
fai = 1:1:90
plot(fai,f(1,fai))
haha = f(1,fai)
求大神指点一二

zan