关于polyval和ode45的问题

如码所示，
想要用MATLAB中的ode45模拟一个变形后的流行病学SIR模型，之前没怎么用过，
里面有大量参数我直接用xlsread读取然后用polyfit和polyval拟合以后想直接代入微分方程参与计算
不知道问题在哪

function dy=NetworkSIRModel(t,y,beta,gamma,NJtoChangZ,NJtoChuZ,NJtoHA,NJtoSH,NJtoSZ,NJtoWX,NJtoYC,NJtoYZ,NJtoZJ,ChangZtoNJ,ChangZtoChuZ,ChangZtoHA,ChangZtoSH,ChangZtoSZ,ChangZtoWX,ChangZtoYC,ChangZtoYZ,ChangZtoZJ,ChuZtoNJ,ChuZtoChangZ,ChuZtoHA,ChuZtoSH,ChuZtoSZ,ChuZtoWX,ChuZtoYC,ChuZtoYZ,ChuZtoZJ,HAtoNJ,HAtoChangZ,HAtoChuZ,HAtoSH,HAtoSZ,HAtoWX,HAtoYC,HAtoYZ,HAtoZJ,SHtoNJ,SHtoChangZ,SHtoChuZ,SHtoHA,SHtoSZ,SHtoWX,SHtoYC,SHtoYZ,SHtoZJ,SZtoNJ,SZtoChangZ,SZtoChuZ,SZtoHA,SZtoSH,SZtoWX,SZtoYC,SZtoYZ,SZtoZJ,WXtoNJ,WXtoChangZ,WXtoChuZ,WXtoHA,WXtoSH,WXtoSZ,WXtoYC,WXtoYZ,WXtoZJ,YCtoNJ,YCtoChangZ,YCtoChuZ,YCtoHA,YCtoSH,YCtoSZ,YCtoWX,YCtoYZ,YCtoZJ,YZtoNJ,YZtoChangZ,YZtoChuZ,YZtoHA,YZtoSH,YZtoSZ,YZtoWX,YZtoYC,YZtoZJ,ZJtoNJ,ZJtoChangZ,ZJtoChuZ,ZJtoHA,ZJtoSH,ZJtoSZ,ZJtoWX,ZJtoYC,ZJtoYZ)
dy=zeros(30,1);
dy(1)=-beta*y(1)*y(11)/y(21)+ChangZtoNJ*y(2)/y(22)+ChuZtoNJ*y(3)/y(23)+HAtoNJ*y(4)/y(24)+SHtoNJ*y(5)/y(25)+SZtoNJ*y(6)/y(26)+WXtoNJ*y(7)/y(27)+YCtoNJ*y(8)/y(28)+YZtoNJ*y(9)/y(29)+ZJtoNJ*y(10)/y(30)-NJtoChangZ*y(1)/y(21)-NJtoChuZ*y(1)/y(21)-NJtoHA*y(1)/y(21)-NJtoSH*y(1)/y(21)-NJtoSZ*y(1)/y(21)-NJtoWX*y(1)/y(21)-NJtoYC*y(1)/y(21)-NJtoYZ*y(1)/y(21)-NJtoZJ*y(1)/y(21); %SNJ
dy(2)=-beta*y(2)*y(12)/y(22)+NJtoChangZ*y(1)/y(21)+ChuZtoChangZ*y(3)/y(23)+HAtoChangZ*y(4)/y(24)+SHtoChangZ*y(5)/y(25)+SZtoChangZ*y(6)/y(26)+WXtoChangZ*y(7)/y(27)+YCtoChangZ*y(8)/y(28)+YZtoChangZ*y(9)/y(29)+ZJtoChangZ*y(10)/y(30)-ChangZtoNJ*y(2)/y(22)-ChangZtoChuZ*y(2)/y(22)-ChangZtoHA*y(2)/y(22)-ChangZtoSH*y(2)/y(22)-ChangZtoSZ*y(2)/y(22)-ChangZtoWX*y(2)/y(22)-ChangZtoYC*y(2)/y(22)-ChangZtoYZ*y(2)/y(22)-ChangZtoZJ*y(2)/y(22); %SChangZ
dy(3)=-beta*y(3)*y(13)/y(23)+NJtoChuZ*y(1)/y(21)+ChangZtoChuZ*y(2)/y(22)++HAtoChuZ*y(4)/y(24)+SHtoChuZ*y(5)/y(25)+SZtoChuZ*y(6)/y(26)+WXtoChuZ*y(7)/y(27)+YCtoChuZ*y(8)/y(28)+YZtoChuZ*y(9)/y(29)+ZJtoChuZ*y(10)/y(30)-ChuZtoNJ*y(3)/y(23)-ChuZtoChangZ*y(3)/y(23)-ChuZtoHA*y(3)/y(23)-ChuZtoSH*y(3)/y(23)-ChuZtoSZ*y(3)/y(23)-ChuZtoWX*y(3)/y(23)-ChuZtoYC*y(3)/y(23)-ChuZtoYZ*y(3)/y(23)-ChuZtoZJ*y(3)/y(23); %SChuz
dy(4)=-beta*y(4)*y(14)/y(24)+NJtoHA*y(1)/y(21)+ChangZtoHA*y(2)/y(22)+ChuZtoHA*y(3)/y(23)+SHtoHA*y(5)/y(25)+SZtoHA*y(6)/y(26)+WXtoHA*y(7)/y(27)+YCtoHA*y(8)/y(28)+YZtoHA*y(9)/y(29)+ZJtoHA*y(10)/y(30)-HAtoNJ*y(4)/y(24)-HAtoChangZ*y(4)/y(24)-HAtoChuZ*y(4)/y(24)-HAtoSH*y(4)/y(24)-HAtoSZ*y(4)/y(24)-HAtoWX*y(4)/y(24)-HAtoYC*y(4)/y(24)-HAtoYZ*y(4)/y(24)-HAtoZJ*y(4)/y(24); %SHA
dy(5)=-beta*y(5)*y(15)/y(25)+NJtoSH*y(1)/y(21)+ChangZtoSH*y(2)/y(22)+ChuZtoSH*y(3)/y(23)+HAtoSH*y(4)/y(24)+SZtoSH*y(6)/y(26)+WXtoSH*y(7)/y(27)+YCtoSH*y(8)/y(28)+YZtoSH*y(9)/y(29)+ZJtoSH*y(10)/y(30)-SHtoNJ*y(5)/y(25)-SHtoChangZ*y(5)/y(25)-SHtoChuZ*y(5)/y(25)-SHtoHA*y(5)/y(25)-SHtoSZ*y(5)/y(25)-SHtoWX*y(5)/y(25)-SHtoYC*y(5)/y(25)-SHtoYZ*y(5)/y(25)-SHtoZJ*y(5)/y(25); %SSH
dy(6)=-beta*y(6)*y(16)/y(26)+NJtoSZ*y(1)/y(21)+ChangZtoSZ*y(2)/y(22)+ChuZtoSZ*y(3)/y(23)+HAtoSZ*y(4)/y(24)+SHtoSZ*y(5)/y(25)+WXtoSZ*y(7)/y(27)+YCtoSZ*y(8)/y(28)+YZtoSZ*y(9)/y(29)+ZJtoSZ*y(10)/y(30)-SZtoNJ*y(6)/y(26)-SZtoChangZ*y(6)/y(26)-SZtoChuZ*y(6)/y(26)-SZtoHA*y(6)/y(26)-SZtoSH*y(6)/y(26)-SZtoWX*y(6)/y(26)-SZtoYC*y(6)/y(26)-SZtoYZ*y(6)/y(26)-SZtoZJ*y(6)/y(26); %SSZ
dy(7)=-beta*y(7)*y(17)/y(27)+NJtoWX*y(1)/y(21)+ChangZtoWX*y(2)/y(22)+ChuZtoWX*y(3)/y(23)+HAtoWX*y(4)/y(24)+SHtoWX*y(5)/y(25)+SZtoWX*y(6)/y(26)+YCtoWX*y(8)/y(28)+YZtoWX*y(9)/y(29)+ZJtoWX*y(10)/y(30)-WXtoNJ*y(7)/y(27)-WXtoChangZ*y(7)/y(27)-WXtoChuZ*y(7)/y(27)-WXtoHA*y(7)/y(27)-WXtoSH*y(7)/y(27)-WXtoSZ*y(7)/y(27)-WXtoYC*y(7)/y(27)-WXtoYZ*y(7)/y(27)-WXtoZJ*y(7)/y(27); %SWX
dy(8)=-beta*y(8)*y(18)/y(28)+NJtoYC*y(1)/y(21)+ChangZtoYC*y(2)/y(22)+ChuZtoYC*y(3)/y(23)+HAtoYC*y(4)/y(24)+SHtoYC*y(5)/y(25)+SZtoYC*y(6)/y(26)+WXtoYC*Y(7)/Y(27)+YZtoYC*y(9)/y(29)+ZJtoYC*y(10)/y(30)-YCtoNJ*y(8)/y(28)-YCtoChangZ*y(8)/y(28)-YCtoChuZ*y(8)/y(28)-YCtoHA*y(8)/y(28)-YCtoSH*y(8)/y(28)-YCtoSZ*y(8)/y(28)-YCtoWX*y(8)/y(28)-YCtoYZ*y(8)/y(28)-YCtoZJ*y(8)/y(28); %SYC
dy(9)=-beta*y(9)*y(19)/y(29)+NJtoYZ*y(1)/y(21)+ChangZtoYZ*y(2)/y(22)+ChuZtoYZ*y(3)/y(23)+HAtoYZ*y(4)/y(24)+SHtoYZ*y(5)/y(25)+SZtoYZ*y(6)/y(26)+WXtoYZ*y(7)/y(27)+YCtoYZ*y(8)/y(28)+ZJtoYZ*y(10)/y(30)-YZtoNJ*y(9)/y(29)-YZtoChangZ*y(9)/y(29)-YZtoChuZ*y(9)/y(29)-YZtoHA*y(9)/y(29)-YZtoSH*y(9)/y(29)-YZtoSZ*y(9)/y(29)-YZtoWX*y(9)/y(29)-YZtoYC*y(9)/y(29)-YZtoZJ*y(9)/y(29); %SYZ
dy(10)=-beta*y(10)*y(20)/y(30)+NJtoZJ*y(1)/y(21)+ChangZtoZJ*y(2)/y(22)+ChuZtoZJ*y(3)/y(23)+HAtoZJ*y(4)/y(24)+SHtoZJ*y(5)/y(25)+SZtoZJ*y(6)/y(26)+WXtoZJ*y(7)/y(27)+YCtoZJ*y(8)/y(28)+YZtoZJ*y(9)/y(29)-ZJtoNJ*y(10)/y(30)-ZJtoChangZ*y(10)/y(30)-ZJtoChuZ*y(10)/y(30)-ZJtoHA*y(10)/y(30)-ZJtoSH*y(10)/y(30)-ZJtoSZ*y(10)/y(30)-ZJtoWX*y(10)/y(30)-ZJtoYC*y(10)/y(30)-ZJtoYZ*y(10)/y(30); %SZJ
dy(11)=beta*y(1)*y(11)/y(21)++ChangZtoNJ*y(12)/y(22)+ChuZtoNJ*y(13)/y(23)+HAtoNJ*y(14)/y(24)+SHtoNJ*y(15)/y(25)+SZtoNJ*y(16)/y(26)+WXtoNJ*y(17)/y(27)+YCtoNJ*y(18)/y(28)+YZtoNJ*y(19)/y(29)+ZJtoNJ*y(20)/y(30)-NJtoChangZ*y(11)/y(21)-NJtoChuZ*y(11)/y(21)-NJtoHA*y(11)/y(21)-NJtoSH*y(11)/y(21)-NJtoSZ*y(11)/y(21)-NJtoWX*y(11)/y(21)-NJtoYC*y(11)/y(21)-NJtoYZ*y(11)/y(21)-NJtoZJ*y(11)/y(21)-gamma*y(11);  %INJ
dy(12)=beta*y(2)*y(12)/y(22)+NJtoChangZ*y(11)/y(21)+ChuZtoChangZ*y(13)/y(23)+HAtoChangZ*y(14)/y(24)+SHtoChangZ*y(15)/y(25)+SZtoChangZ*y(16)/y(26)+WXtoChangZ*y(17)/y(27)+YCtoChangZ*y(18)/y(28)+YZtoChangZ*y(19)/y(29)+ZJtoChangZ*y(20)/y(30)-ChangZtoNJ*y(12)/y(22)-ChangZtoChuZ*y(12)/y(22)-ChangZtoHA*y(12)/y(22)-ChangZtoSH*y(12)/y(22)-ChangZtoSZ*y(12)/y(22)-ChangZtoWX*y(12)/y(22)-ChangZtoYC*y(12)/y(22)-ChangZtoYZ*y(12)/y(22)-ChangZtoZJ*y(12)/y(22)-gamma*y(12); %IChangZ
dy(13)=beta*y(3)*y(13)/y(23)+NJtoChuZ*y(11)/y(21)+ChangZtoChuZ*y(12)/y(22)++HAtoChuZ*y(14)/y(24)+SHtoChuZ*y(15)/y(25)+SZtoChuZ*y(16)/y(26)+WXtoChuZ*y(17)/y(27)+YCtoChuZ*y(18)/y(28)+YZtoChuZ*y(19)/y(29)+ZJtoChuZ*y(20)/y(30)-ChuZtoNJ*y(13)/y(23)-ChuZtoChangZ*y(13)/y(23)-ChuZtoHA*y(13)/y(23)-ChuZtoSH*y(13)/y(23)-ChuZtoSZ*y(13)/y(23)-ChuZtoWX*y(13)/y(23)-ChuZtoYC*y(13)/y(23)-ChuZtoYZ*y(13)/y(23)-ChuZtoZJ*y(13)/y(23)-gamma*y(13); %IChuZ
dy(14)=beta*y(4)*y(14)/y(24)+NJtoHA*y(11)/y(21)+ChangZtoHA*y(12)/y(22)+ChuZtoHA*y(13)/y(23)+SHtoHA*y(15)/y(25)+SZtoHA*y(16)/y(26)+WXtoHA*y(17)/y(27)+YCtoHA*y(18)/y(28)+YZtoHA*y(19)/y(29)+ZJtoHA*y(20)/y(30)-HAtoNJ*y(14)/y(24)-HAtoChangZ*y(14)/y(24)-HAtoChuZ*y(14)/y(24)-HAtoSH*y(14)/y(24)-HAtoSZ*y(14)/y(24)-HAtoWX*y(14)/y(24)-HAtoYC*y(14)/y(24)-HAtoYZ*y(14)/y(24)-HAtoZJ*y(14)/y(24)-gamma*y(14); %IHA
dy(15)=beta*y(5)*y(15)/y(25)+NJtoSH*y(11)/y(21)+ChangZtoSH*y(12)/y(22)+ChuZtoSH*y(13)/y(23)+HAtoSH*y(14)/y(24)+SZtoSH*y(16)/y(26)+WXtoSH*y(17)/y(27)+YCtoSH*y(18)/y(28)+YZtoSH*y(19)/y(29)+ZJtoSH*y(20)/y(30)-SHtoNJ*y(15)/y(25)-SHtoChangZ*y(15)/y(25)-SHtoChuZ*y(15)/y(25)-SHtoHA*y(15)/y(25)-SHtoSZ*y(15)/y(25)-SHtoWX*y(15)/y(25)-SHtoYC*y(15)/y(25)-SHtoYZ*y(15)/y(25)-SHtoZJ*y(15)/y(25)-gamma*y(15);  %ISH
dy(16)=beta*y(6)*y(16)/y(26)+NJtoSZ*y(11)/y(21)+ChangZtoSZ*y(12)/y(22)+ChuZtoSZ*y(13)/y(23)+HAtoSZ*y(14)/y(24)+SHtoSZ*y(15)/y(25)+WXtoSZ*y(17)/y(27)+YCtoSZ*y(18)/y(28)+YZtoSZ*y(19)/y(29)+ZJtoSZ*y(20)/y(30)-SZtoNJ*y(16)/y(26)-SZtoChangZ*y(16)/y(26)-SZtoChuZ*y(16)/y(26)-SZtoHA*y(16)/y(26)-SZtoSH*y(16)/y(26)-SZtoWX*y(16)/y(26)-SZtoYC*y(16)/y(26)-SZtoYZ*y(16)/y(26)-SZtoZJ*y(16)/y(26)-gamma*y(16); %ISZ
dy(17)=beta*y(7)*y(17)/y(27)+NJtoWX*y(11)/y(21)+ChangZtoWX*y(12)/y(22)+ChuZtoWX*y(13)/y(23)+HAtoWX*y(14)/y(24)+SHtoWX*y(15)/y(25)+SZtoWX*y(16)/y(26)+YCtoWX*y(18)/y(28)+YZtoWX*y(19)/y(29)+ZJtoWX*y(20)/y(30)-WXtoNJ*y(17)/y(27)-WXtoChangZ*y(17)/y(27)-WXtoChuZ*y(17)/y(27)-WXtoHA*y(17)/y(27)-WXtoSH*y(17)/y(27)-WXtoSZ*y(17)/y(27)-WXtoYC*y(17)/y(27)-WXtoYZ*y(17)/y(27)-WXtoZJ*y(17)/y(27)-gamma*y(17); %IWX
dy(18)=beta*y(8)*y(18)/y(28)+NJtoYC*y(11)/y(21)+ChangZtoYC*y(12)/y(22)+ChuZtoYC*y(13)/y(23)+HAtoYC*y(14)/y(24)+SHtoYC*y(15)/y(25)+SZtoYC*y(16)/y(26)+WXtoYC*Y(17)/Y(27)+YZtoYC*y(19)/y(29)+ZJtoYC*y(20)/y(30)-YCtoNJ*y(18)/y(28)-YCtoChangZ*y(18)/y(28)-YCtoChuZ*y(18)/y(28)-YCtoHA*y(18)/y(28)-YCtoSH*y(18)/y(28)-YCtoSZ*y(18)/y(28)-YCtoWX*y(18)/y(28)-YCtoYZ*y(18)/y(28)-YCtoZJ*y(18)/y(28)-gamma*y(18); %IYC
dy(19)=beta*y(9)*y(19)/y(29)+NJtoYZ*y(11)/y(21)+ChangZtoYZ*y(12)/y(22)+ChuZtoYZ*y(13)/y(23)+HAtoYZ*y(14)/y(24)+SHtoYZ*y(15)/y(25)+SZtoYZ*y(16)/y(26)+WXtoYZ*y(17)/y(27)+YCtoYZ*y(18)/y(28)+ZJtoYZ*y(20)/y(30)-YZtoNJ*y(19)/y(29)-YZtoChangZ*y(19)/y(29)-YZtoChuZ*y(19)/y(29)-YZtoHA*y(19)/y(29)-YZtoSH*y(19)/y(29)-YZtoSZ*y(19)/y(29)-YZtoWX*y(19)/y(29)-YZtoYC*y(19)/y(29)-YZtoZJ*y(19)/y(29)-gamma*y(19); %IYZ
dy(20)=beta*y(10)*y(20)/y(30)+NJtoZJ*y(11)/y(21)+ChangZtoZJ*y(12)/y(22)+ChuZtoZJ*y(13)/y(23)+HAtoZJ*y(14)/y(24)+SHtoZJ*y(15)/y(25)+SZtoZJ*y(16)/y(26)+WXtoZJ*y(17)/y(27)+YCtoZJ*y(18)/y(28)+YZtoZJ*y(19)/y(29)-ZJtoNJ*y(20)/y(30)-ZJtoChangZ*y(20)/y(30)-ZJtoChuZ*y(20)/y(30)-ZJtoHA*y(20)/y(30)-ZJtoSH*y(20)/y(30)-ZJtoSZ*y(20)/y(30)-ZJtoWX*y(20)/y(30)-ZJtoYC*y(20)/y(30)-ZJtoYZ*y(20)/y(30)-gamma*y(20); %IZJ
dy(21)=ChangZtoNJ+ChuZtoNJ+HAtoNJ+SHtoNJ+SZtoNJ+WXtoNJ+YCtoNJ+YZtoNJ+ZJtoNJ-NJtoChangZ-NJtoChuZ-NJtoHA-NJtoSH-NJtoSZ-NJtoWX-NJtoYC-NJtoYZ-NJtoZJ; %NNJ
dy(22)=NJtoChangZ+ChuZtoChangZ+HAtoChangZ+SHtoChangZ+SZtoChangZ+WXtoChangZ+YCtoChangZ+YZtoChangZ+ZJtoChangZ-ChangZtoNJ-ChangZtoChuZ-ChangZtoHA-ChangZtoSH-ChangZtoSZ-ChangZtoWX-ChangZtoYC-ChangZtoYZ-ChangZtoZJ; %NChangZ
dy(23)=NJtoChuZ+ChangZtoChuZ++HAtoChuZ+SHtoChuZ+SZtoChuZ+WXtoChuZ+YCtoChuZ+YZtoChuZ+ZJtoChuZ-ChuZtoNJ-ChuZtoChangZ-ChuZtoHA-ChuZtoSH-ChuZtoSZ-ChuZtoWX-ChuZtoYC-ChuZtoYZ-ChuZtoZJ; %NChuZ
dy(24)=NJtoHA+ChangZtoHA+ChuZtoHA+SHtoHA+SZtoHA+WXtoHA+YCtoHA+YZtoHA+ZJtoHA-HAtoNJ-HAtoChangZ-HAtoChuZ-HAtoSH-HAtoSZ-HAtoWX-HAtoYC-HAtoYZ-HAtoZJ; %NHA
dy(25)=NJtoSH+ChangZtoSH+ChuZtoSH+HAtoSH+SZtoSH+WXtoSH+YCtoSH+YZtoSH+ZJtoSH-SHtoNJ-SHtoChangZ-SHtoChuZ-SHtoHA-SHtoSZ-SHtoWX-SHtoYC-SHtoYZ-SHtoZJ; %NSH
dy(26)=NJtoSZ+ChangZtoSZ+ChuZtoSZ+HAtoSZ+SHtoSZ+WXtoSZ+YCtoSZ+YZtoSZ+ZJtoSZ-SZtoNJ-SZtoChangZ-SZtoChuZ-SZtoHA-SZtoSH-SZtoWX-SZtoYC-SZtoYZ-SZtoZJ; %NSZ
dy(27)=NJtoWX+ChangZtoWX+ChuZtoWX+HAtoWX+SHtoWX+SZtoWX+YCtoWX+YZtoWX+ZJtoWX-WXtoNJ-WXtoChangZ-WXtoChuZ-WXtoHA-WXtoSH-WXtoSZ-WXtoYC-WXtoYZ-WXtoZJ; %NWX
dy(28)=NJtoYC+ChangZtoYC+ChuZtoYC+HAtoYC+SHtoYC+SZtoYC+WXtoYC+YZtoYC+ZJtoYC-YCtoNJ-YCtoChangZ-YCtoChuZ-YCtoHA-YCtoSH-YCtoSZ-YCtoWX-YCtoYZ-YCtoZJ; %NYC
dy(29)=NJtoYZ+ChangZtoYZ+ChuZtoYZ+HAtoYZ+SHtoYZ+SZtoYZ+WXtoYZ+YCtoYZ+ZJtoYZ-YZtoNJ-YZtoChangZ-YZtoChuZ-YZtoHA-YZtoSH-YZtoSZ-YZtoWX-YZtoYC-YZtoZJ; %NYZ
dy(30)=NJtoZJ+ChangZtoZJ+ChuZtoZJ+HAtoZJ+SHtoZJ+SZtoZJ+WXtoZJ+YCtoZJ+YZtoZJ-ZJtoNJ-ZJtoChangZ-ZJtoChuZ-ZJtoHA-ZJtoSH-ZJtoSZ-ZJtoWX-ZJtoYC-ZJtoYZ; %NZJ
end

以上是微分方程组的代码

下面是执行部分

clear
tspan=[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59];
y0=[8085241 4067856 2466000 2992223 24870885 10419300 6178208 4301656 3238904 2550827 0 0 0 0 10 0 0 0 0 0 8085241 4067856 2466000 2992223 24870895 10419300 6178208 4301656 3238904 2550827]; % SNJ SChangZ SChuZ SHA SSH SSZ SWX SYC SYZ SZJ INJ IChangZ IChuZ IHA ISH ISZ IWX IYC IYZ IZJ NNJ NChangZ NChuZ NHA NSH NSZ NWX NYC NYZ NZJ
beta=1.0000;
gamma=0.692;

k=xlsread('常州迁入南京','Sheet1','H2:H60');
x=[k]';
a=polyfit(tspan,x,5);
ChangZtoNJ=polyval(a,tspan);

%中间全部是类似的数据代入有将近一百个我就全部省略了没有贴上来。。。

k=xlsread('扬州迁入镇江','Sheet1','H2:H60');
x=[k]';
a=polyfit(tspan,x,5);
YZtoZJ=polyval(a,tspan);

[T,Y]=ode45(@(t,y)NetworkSIRModel_v0(t,y,beta,gamma,NJtoChangZ,NJtoChuZ,NJtoHA,NJtoSH,NJtoSZ,NJtoWX,NJtoYC,NJtoYZ,NJtoZJ,ChangZtoNJ,ChangZtoChuZ,ChangZtoHA,ChangZtoSH,ChangZtoSZ,ChangZtoWX,ChangZtoYC,ChangZtoYZ,ChangZtoZJ,ChuZtoNJ,ChuZtoChangZ,ChuZtoHA,ChuZtoSH,ChuZtoSZ,ChuZtoWX,ChuZtoYC,ChuZtoYZ,ChuZtoZJ,HAtoNJ,HAtoChangZ,HAtoChuZ,HAtoSH,HAtoSZ,HAtoWX,HAtoYC,HAtoYZ,HAtoZJ,SHtoNJ,SHtoChangZ,SHtoChuZ,SHtoHA,SHtoSZ,SHtoWX,SHtoYC,SHtoYZ,SHtoZJ,SZtoNJ,SZtoChangZ,SZtoChuZ,SZtoHA,SZtoSH,SZtoWX,SZtoYC,SZtoYZ,SZtoZJ,WXtoNJ,WXtoChangZ,WXtoChuZ,WXtoHA,WXtoSH,WXtoSZ,WXtoYC,WXtoYZ,WXtoZJ,YCtoNJ,YCtoChangZ,YCtoChuZ,YCtoHA,YCtoSH,YCtoSZ,YCtoWX,YCtoYZ,YCtoZJ,YZtoNJ,YZtoChangZ,YZtoChuZ,YZtoHA,YZtoSH,YZtoSZ,YZtoWX,YZtoYC,YZtoZJ,ZJtoNJ,ZJtoChangZ,ZJtoChuZ,ZJtoHA,ZJtoSH,ZJtoSZ,ZJtoWX,ZJtoYC,ZJtoYZ),tspan,y0);
plot(T,Y(:,11),T,Y(:,12),'--',T,Y(:,13),T,Y(:,14),'--',T,Y(:,15),T,Y(:,16),'--',T,Y(:,17),T,Y(:,18),'--',T,Y(:,19),'--',T,Y(:,20),'--');
legend('NJ','ChangZ','ChuZ','HA','SH','SZ','WX','YC','YZ','ZJ')
xlabel('Time(Days)');
ylabel('Number of infectous') 

结果显示无法赋值，左右侧元素不同
无法运行！可以帮看看问题出在哪里了嘛！

[t,y] = ode45(odefun,tspan,y0)
功能为：求微分方程组 y′=f(t,y) 从 t0 到 tf 的积分（y可以是个列向量）
这里求的是数值解，y没有显示公式表达，tspan = [t0 tf]为求解时间，

检查一下你的返回值