最小二乘法参数辨识直线电机动力学模型参数的matlab代码，要画图

好的，以下是最小二乘法参数辨识直线电机动力学模型参数的 MATLAB 代码，并包含画图功能：

%% 读取数据
data = readmatrix('data.csv'); % 读取数据文件，假设数据文件名为 data.csv
t = data(:, 1); % 时间
u = data(:, 2); % 输入电压
y = data(:, 3); % 位移

%% 建立模型
A = [u, y]; % 构建回归矩阵
B = regress(y, A); % 最小二乘法求解参数

%% 画图
subplot(2, 1, 1);
plot(t, y, 'b', t, B(1)*u + B(2), 'r');
xlabel('时间 (s)');
ylabel('位移 (mm)');
legend('实际位移', '模型预测');
title('位移随时间变化');

subplot(2, 1, 2);
plot(u, y, 'b', u, B(1)*u + B(2), 'r');
xlabel('输入电压 (V)');
ylabel('位移 (mm)');
legend('实际位移', '模型预测');
title('位移随输入电压变化');

suptitle('直线电机动力学模型参数辨识');

其中，数据文件需要自己准备，格式为三列：时间，输入电压和位移。此代码通过最小二乘法求解直线电机动力学模型的两个参数，并绘制了位移随时间和输入电压变化的图像，其中红色曲线表示模型预测值，蓝色曲线表示实际位移值。