www.gusucode.com > 司守奎《Lingo软件及应用习题解答》的源程序 > 司守奎《Lingo软件及应用习题解答》的源程序/《Lingo软件及应用习题解答》程序和数据/12第12章/Lt12_3_6.m

    clc, clear
r1=-log(1+0.018), r2=-log(1+2*0.01944)
r3=-log(1+3*0.0216), r5=-log(1+5*0.02304)
a=zeros(11); %邻接矩阵初始化;
b=[r1,r2,r3,0,r5,zeros(1,5)]
for i=1:10
    a(i,[i+1:11])=b(1:11-i);
end
a %显示邻接矩阵
b=sparse(a); %转换成MATLAB工具箱需要的稀疏矩阵
[d,path]=graphshortestpath(b,1,[2:11],'Method','Bellman-Ford')
celldisp(path) %显示细胞数组的全部元素
str=strcat('v',int2str([1:11]'));
view(biograph(b,str,'ShowWeights','on','LayoutType','equilibrium','EdgeType','segmented'))
e=exp(-d) %显示单位资金从初始到各年度的本息和
xlswrite('Ltdata12_3_1.xlsx',e)