我无法弄清楚为什么我的循环不在matlab中

Question

在MATLAB中，我试图编写一个程序，该程序将在图（x，y）上使用3个坐标，使用这些值来求解一个方程系统，该方程将找到多项式方程的系数，y = ax^2 + BX + C，然后我可以用来绘制抛物线。

为了测试我的代码，我认为我可以从多项式开始，绘制图形，找到抛物线的最小位置，将其直接邻居用于我的其他两个位置，然后通过我的代码运行这3个位置，这应该吐出的系数我原来的多项式。但是由于某种原因，我由此产生的抛物线是正确的转移，而我对B和C的价值是不正确的。

有人看到我的问题在哪里吗？我没有想法

clear all; close all;

x = -10:10;

%Original Polynomial
y = 6.*x.^2 + 11.*x -35;

% Find 3 Locations
[max_n, max_i] = min(y)
max_il = max_i - 1                 % left neighbor of max_ni
max_nl = y(max_il)                 % value at max_il
max_ir = max_i + 1                 % left neighbor of max_ni
max_nr = y(max_ir)                 % value at max_ir

% Solve for coefficients
syms a b c
equ = (a)*(max_i)^2 + (b)*(max_i) + (c) == (max_n);
equ_l = (a)*(max_il)^2 + (b)*(max_il) + (c) == (max_nl);
equ_r = (a)*(max_ir)^2 + (b)*(max_ir) + (c) == (max_nr);

sol = solve([equ, equ_l, equ_r],[a, b, c]);

Sola = sol.a
Solb = sol.b
Solc = sol.c

% New Polynomial
p = (sol.a).*(x).^2 + (sol.b).*(x) +(sol.c);


%Plot
plot(x,y); grid on; hold on;
plot(x, p); 
axis([-10 10 -41 40])
[max_np, max_ip] = min(p)
legend('OG', 'New')

Answer 1

您将索引混淆到您的数组y和相应的x坐标。

x = -10:10;
y = 6.*x.^2 + 11.*x -35;

[max_n, max_i] = min(y)

这里。 max_i是y数组中的索引，相应的X坐标将为x（max_i）。

我建议您找到三个数据点，以适合您的曲线如下：

[~, max_i] = min(y);
pts_x = x(max_i + (-1:1));
pts_y = y(max_i + (-1:1));

然后使用pts_x（i）和pts_y（i）作为您的x和y值：

syms a b c
equ = a * pts_x.^2 + b * pts_x + c == pts_y;
sol = solve(equ, [a, b, c]);