最便宜可交割债券(CTD)求解

源于需求：

解决需求：

%% FindCTD
%{
测试输入：
futureput =[101 12]
bondput=[125 05 1.2131
142 15 1.3792
115 31 1.1149
144 02 1.4026];
%}
judeg = zeros(1,1);
realprice = zeros(1,1);
futureput = '请输入当前国债期货价格： ';
futurep = input(futureput);%也可以输入矩阵，注意得[a b c;a b c;a b c]这样输入
future = futurep(1,1) + futurep(1,2);
bondput = '请输入需要计算的债券矩阵： ';
bond = input(bondput);%也可以输入矩阵，注意得[a b c;a b c;a b c]这样输入

[row,col]=size(bond);  %获取矩阵的行数和列数,row为行数，col为猎术
for i=1:row
     realprice(i,1) = bond(i,1)+bond(i,2)/32;
     judge(i,1) = realprice(i,1) - future * bond(i,3);
end
CTD=min(judge(:));
[x,y]=find(judge==CTD);
X = [' CTD是第  ',num2str(x),'  支债券。'];
disp(X)

——Written in Liangxiang