好吧，我遇到了一个愚蠢的问题，我自己的工作和互联网搜索都没有真正提供解决方案，所以我想问问有没有人有任何建议。我正在编写一个 Matlab 程序，该程序最终旨在绘制利率曲线，这是通过使用三次 B 样条，然后使用优化器来获得三次 B 样条基础中的系数，并使用此信息绘制图形来完成的。现在我在尝试绘制图表时遇到了一个棘手的小问题。我将在下面提供错误消息和相关代码，但总的来说，发生的事情是这样的：
要对给定度数的 B 样条进行编程，您可以从指标函数开始，然后递归地构建它们。这个指标是关于一些输入变量的（对我来说是时间，但变量实际上是什么并不是特别相关）。当我稍后针对这个变量绘制图表时，使用像 t = (0:1:30) 这样的范围，我收到一个错误，告诉我 || 的操作数 和 && 运算符必须可转换为逻辑标量值。
我进行了一些互联网搜索，发现 && 不适用于矢量输入（我知道我会在尝试绘制图形时将其传递给它），并建议使用 & 代替。我试过了，而是收到一个错误，说我的矩阵尺寸不同意，这对我来说没有意义，因为首先从来没有涉及任何矩阵。所以在这一点上，我不确定解决这个问题的好方法是什么。在 matlab 中是否有一种向量稳健的编程指标函数方法可以进行绘图？
以下是相关代码：
startvec = 0.007*ones(16,2);
%[x,fval] = fminsearch(@(x) FittingFunction(x, Rhat, today, startdays, enddays, fixedcoupondays, floatingcoupondays),startvec,optimset('Display','iter','MaxIter',50,'MaxFunEvals',50));
f0 = 0;
l0 = 0;
t = (0:1:30);
x = startvec;
for k=-3:1:12
f0 = f0 + x(k+4,1).*Splines(3,k,t);
l0 = l0 + x(k+4,2).*Splines(3,k,t);
end
h1 = figure('Name','OIS Curve','NumberTitle','off')
plot(t,f0);
title('OIS Curve as a Function of Time');
xlabel('time (years)');
ylabel('OIS Rate');
%hold off
print(h1,'-dpdf','OISCurve.pdf')
h2 = figure('Name','LIBOR Curve','NumberTitle','off')
plot(t,l0);
title('LIBOR Curve as a Function of Time');
xlabel('time (years)');
ylabel('LIBOR Rate');
%hold off
print(h2,'-dpdf','LIBORCurve.pdf')
函数 Splines 写在这里：
function spline = Splines(n,k,t)
% The following code builds a B-spline function of any degree, using the
% recursive properties of B-splines. Can be defined from k=-3 to k=12.
knots = [-1/4 -1/6 -1/12 0 1/12 1/4 1/2 3/4 1 3/2 2 4 7 11 21 31 41 51 61 71 81 91 101 111];
if (n==0)
if (t>knots(k+4) && t<knots(k+5))
spline = 1;
else
spline = 0;
end
else
spline = ((t-knots(k+4))/(knots(k+4+n)-knots(k+4)))*Splines(n-1,k,t) + ((knots(k+5+n)-t)/(knots(k+5+n)-knots(k+5)))*Splines(n-1,k+1,t);
end
我使用 && 时收到的错误消息是：
|| 的操作数 和 && 运算符必须可转换为逻辑标量值。
如果 (t>knots(k+4) && t) 样条曲线错误（第 9 行）
样条线错误（第 15 行）样条线 = ((t-knots(k+4))/(knots(k+4+n)-knots(k+4)))*Splines(n-1,k,t) + ((结(k+5+n)-t)/(结(k+5+n)-结(k+5)))*样条(n-1,k+1,t);
样条线错误（第 15 行）样条线 = ((t-knots(k+4))/(knots(k+4+n)-knots(k+4)))*Splines(n-1,k,t) + ((结(k+5+n)-t)/(结(k+5+n)-结(k+5)))*样条(n-1,k+1,t);
样条线错误（第 15 行）样条线 = ((t-knots(k+4))/(knots(k+4+n)-knots(k+4)))*Splines(n-1,k,t) + ((结(k+5+n)-t)/(结(k+5+n)-结(k+5)))*样条(n-1,k+1,t);
InterestRatesHW1 错误（第 88 行） f0 = f0 + x(k+4,1).*Splines(3,k,t);
我只使用 & 时得到的错误是：
错误使用 * 内矩阵尺寸必须一致。
样条线错误（第 15 行）样条线 = ((t-knots(k+4))/(knots(k+4+n)-knots(k+4)))*Splines(n-1,k,t) + ((结(k+5+n)-t)/(结(k+5+n)-结(k+5)))*样条(n-1,k+1,t);
样条线错误（第 15 行）样条线 = ((t-knots(k+4))/(knots(k+4+n)-knots(k+4)))*Splines(n-1,k,t) + ((结(k+5+n)-t)/(结(k+5+n)-结(k+5)))*样条(n-1,k+1,t);
InterestRatesHW1 错误（第 88 行） f0 = f0 + x(k+4,1).*Splines(3,k,t);
我很抱歉问什么是我确定的基本问题，但这变得非常令人沮丧，因此非常感谢任何帮助。