第 1 章：你的第一个数据分析¶

% 输入实验数据
F = [0, 1, 2, 3, 4, 5];        % 拉力，单位：牛顿
x = [0, 1.2, 2.3, 3.5, 4.6, 5.8];  % 伸长量，单位：厘米

% 计算弹性系数 k = F / x（胡克定律）
k = F ./ x;                     % 逐元素除法

% 查看结果
disp('弹性系数 k (N/cm)：');
disp(k);

a = 42;              % 标量（单个数值）
name = 'MATLAB';     % 字符串（用单引号）
isReady = true;      % 逻辑值

% 分号 ; 抑制输出，不加分号则显示结果
b = 3.14             % 不加分号，会在命令窗口显示 b = 3.1400

% 行向量
row = [1, 2, 3, 4, 5];       % 逗号分隔
row = [1 2 3 4 5];            % 空格分隔也行

% 列向量
col = [1; 2; 3; 4; 5];       % 分号表示换行

% 矩阵
A = [1, 2, 3;
     4, 5, 6;
     7, 8, 9];

% 快速生成序列
x = 1:10;                     % 1 到 10，步长为 1
x = 1:2:10;                   % 1, 3, 5, 7, 9（步长为 2）
x = linspace(0, 1, 100);      % 0 到 1 之间均匀取 100 个点

A = [10, 20, 30;
     40, 50, 60;
     70, 80, 90];

A(2, 3)        % 第 2 行第 3 列 → 60
A(1, :)        % 第 1 行所有列 → [10, 20, 30]
A(:, 2)        % 所有行第 2 列 → [20; 50; 80]
A(1:2, 2:3)    % 第 1~2 行，第 2~3 列 → [20, 30; 50, 60]

A = [1, 2; 3, 4];
B = [5, 6; 7, 8];

A * B           % 矩阵乘法 → [19, 22; 43, 50]
A .* B          % 逐元素乘法 → [5, 12; 21, 32]（注意点号 .）

A / B           % 矩阵右除（A * inv(B)）
A ./ B          % 逐元素除法 → [0.2, 0.333; 0.429, 0.5]

A ^ 2           % 矩阵幂（A * A）
A .^ 2          % 逐元素幂 → [1, 4; 9, 16]

x = [0, pi/2, pi];

sin(x)          % 正弦 → [0, 1, 1.2246e-16]
cos(x)          % 余弦 → [1, 6.1232e-17, -1]
exp(x)          % e^x
log(x)          % 自然对数
sqrt(x)         % 平方根
abs(x)          % 绝对值
round(x)        % 四舍五入

% spring_analysis.m —— 弹簧弹性系数分析
clear; clc; close all;  % 清理工作环境

% 1. 输入实验数据
F = [0, 1, 2, 3, 4, 5];
x = [0, 1.2, 2.3, 3.5, 4.6, 5.8];

% 2. 计算弹性系数（排除第一个 0/0 点）
k = F(2:end) ./ x(2:end);

% 3. 统计分析
k_mean = mean(k);
k_std = std(k);

% 4. 输出结果
fprintf('弹性系数均值：%.4f N/cm\n', k_mean);
fprintf('弹性系数标准差：%.4f N/cm\n', k_std);