文件是在路径里的，但运行以后就是显示未找到文件，无法引用EXCEL里的数据，请各位大佬帮忙看一下

下面是程序(我用的是MATLAB64位2012a，电脑是windows10,EXCEL是2010版本)

clear;

clc;

[a,b,c]=xlsread('fujian1-3.xlsx','fujian1','A4:C24');

[~,~,n]=xlsread('fujian1-3.xIsx','fujian1','A2');%%根据数据拟合出曲线

X=a(:,1);

y=a(:,2);

plot(x,y);

p=polyfit(x,y,2);

yl=polyval(p,x);

plot(x,y,'*r',x,y1,'-b');

legend('题目所给影端的位置','拟合得到的影端曲线');

title('太阳影长变化曲线');

xlabel('横坐标X(米)');

ylabel('纵坐标Y(米)');%%计算出影长

l=sqrt(x.^2+y.^2);

for i=1:length(c)

T{i}=datevec(c{i,1},'HH:MM:SS');

t(i,1)=T{i}(4)+T{i}(5)/60+T{i}(6)/3600;

end

%% 日期的处理

d=datenum(n,'yy/mm/dd');

d0=datenum('2015/1/l','yy/mm/dd');

date=d-d0+l;%n为一年中日期的序号

%%寻找经纬度

m=1;

CAL=zeros(21,2);

e0=zeros(1,2);

H=zeros(1,2);

FAI=zeros(1,2);

SEITA=zeros(1,2);

for h=1:10

for fai=0:0.6:60

for seita=15:0.6:150

cal=(h*(1-(sin((469*pi*sin((2*pi*(date + 284))/365))/3600)*sin((pi* fai)/180)+cos((pi*fai)/180)*cos((pi*(15*abs(seita/15+t-8)-180))/180)*cos((469*pi*sin((2*pi*(date+284))/365))/3600))^2)^(1/2)/(sin((469*pi*sin((2*pi*(date+284))/365))/3600)*sin((pi*fai)/180)+cos((pi* fai)/180)*cos((pi*(15*abs(seita/15+t-8)-180))/180)*cos((469*pi*sin((2*pi*(date+284))/365))/3600)));

error=sum(abs(cal-l)./l)/length(l);

if error<0.002

CAL(:,m)=cal;

e0(m)=error;

H(m)=h;

FAI(m)=fai;

SEITA(m)= -seita;

m=m+1;

end

end

end

end

%%检验所求得经纬度下影长与实际影长的关系

%第一个点

T=t;

figure(2);

subplot(1,2,1)

seita=SEITA(1);

fai=FAI(1);

h=H(1);

t=sym('t');

L1=subs(L);

for i=1:21

t=T(i);

l1(i)=eval(L1);

end

stem(T,l1-l')

title('A地所求影子长度的差值');

xlabel('北京时间');

xlabel('北京时间');

ylabel('差值');

% ezplot(L1,[T(1),T(21)]);

% title('东经109.2°北纬19.2°点太阳影长与北京时间的关系');

% xlabel ('北京时间');

% ylabel ('影子的长度');

% hold on

% plot(T,l,'r');

% legend('所求得地点的影长与时间的关系曲线','题目数据拟合出的影长与时间曲线',0);

% t=T;

% xi=14.7:0.05:15.7;

% yi=interpl(t,eval(L),xi,'spline');

% alfa=atan(h./yi);

% alfa=alfa*180/pi;

% 第二个点

subplot(1,2,2);

seita=SEITA(2);

fai=FAI(2);

h=H(2);

t=sym('t');

L2=subs(L);

for i=1:21

t=T(i);

l2(i)=eval(L2);

end

stem(T,l2-l')

titl('B地所求影子长度的差值');

xlabel('北京时间');

ylabel('差值');

% ezplot(L2,[T(1),T(21)]);

% title('东经95.4，北纬25.2点太阳影长与北京时间的关系');

% xlabel('北京时间');

% ylabel('影子的长度');

% hold on;

% plot(T,l,'r');

% legend('所求得地点的影长与时间的关系曲线','题目数据拟合出的影长与时间曲线',0);

% xi=14.7:0.05:15.7;

% t=T;

% yi=interpl(t,eval(L),xi,'spline');

% alfa=atan(h./yi);

% alfa=alfa*180/pi;

拟牛顿.png

(358.26 KB, 下载次数: 0)

2019-7-29 14:47 上传