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)