结果实在是跑不对
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)
求大神指点一二
不是我敷衍你回答，这里大部分人都是干互联网技术，你的问题是数学专业型很强的问题，即使我是学生也没用过啊
（1）f(j,fai) 函数最好写一个.m文件的脚本
（2）fai = 1:1:90 ， plot(fai,f(1,fai)) ， haha = f(1,fai) 你这个是用for循环画图吧，你在同一个图上画多个图得用hold on，然后，好像有for...end
循环代码块吧。
（3）其它的，你这matlab代码写的可读性很差（数组不能四舍五入？整个函数可以写一个.m文件吧，不用再command窗口中写命令吧），自己好好学吧，