为何Matlab代码无报错却未生成并保存应力应变图？

你的MATLAB代码无报错但未生成图像、无法保存文件，核心问题集中在数据读取逻辑错误、保存路径未提前创建、绘图覆盖/数据有效性不足这几点，以下是具体修复步骤和优化建议：

1. 修正数据读取逻辑 #

你使用textscan(fileID, '%f')仅读取单列数据，但columnIndices指定了3列（对应Crosshead、Load、Time），这会导致后续提取列数据时要么拿到错误结果，要么因索引超出范围隐性失败（若文件实际为多列）。需根据文件分隔符修改textscan格式符，并改用内置参数跳过表头：

% 空格/制表符分隔的3列数据，若为逗号分隔则改为 '%f,%f,%f'
data = textscan(fileID, '%f %f %f', 'HeaderLines', linesToSkip);

同时可以删除手动跳过8行的循环，HeaderLines参数更简洁可靠。

2. 确保保存目录存在 #

MATLAB的saveas不会自动创建不存在的文件夹，需提前检查并创建目标目录：

saveDir = fullfile(directory, 'Stress-strain curve');
if ~exist(saveDir, 'dir')
    mkdir(saveDir);
end

另外原代码会生成xxx.txt.png这类冗余文件名，建议截取原文件名去掉.txt后缀：

saveas(gcf, fullfile(saveDir, [filename(1:end-4) '.png']));

3. 优化绘图逻辑避免覆盖 #

每次循环使用figure(1)会覆盖之前的图像，若需保留所有文件的曲线，应创建新的绘图窗口：

% 替换原figure(1)，每次生成新窗口
fig = figure;
plot(Strain_length, Stress, 'b')
hold on 
plot(Strain_length, f, 'r')
xlabel('应变')
ylabel('应力 [MPa]')
title(['试样: ', filename])
legend('应力-应变曲线','弹性阶段拟合')
hold off
% 用窗口句柄保存更可靠
saveas(fig, fullfile(saveDir, [filename(1:end-4) '.png']));
close(fig); % 关闭窗口节省内存

4. 增加数据有效性检查 #

避免因数据长度不足（如仅1个数据点）导致无有效图像：

if length(crossheadColumn) < 2 || length(loadColumn) < 2
    warning(['文件 ', filename, ' 数据长度不足，跳过处理']);
    continue;
end

5. 优化杨氏模量计算 #

直接用全量数据拟合会包含非线性区域，建议仅取弹性阶段数据（如应变≤0.01）计算模量：

elasticIdx = Strain_length <= 0.01; % 弹性应变阈值可根据试验调整
if sum(elasticIdx) < 2
    warning(['文件 ', filename, ' 弹性阶段数据不足，无法计算模量']);
    continue;
end
p = polyfit(Strain_length(elasticIdx), Stress(elasticIdx), 1);
f = polyval(p, Strain_length);
E = p(1);
all_E = [all_E; E];

完整修复后代码示例 #

% 定义数据文件目录
directory = 'C:\Users\Biologisch';

% 获取目录下所有文件
files = dir(directory);

% 配置提取列与名称
columnIndices = [1, 2, 3];
listNames = {'Crosshead', 'Load', 'Time'};
all_E = [];

% 试样参数
W0 = 45;               % 初始宽度 [mm]
D = 4*10^-3;           % 厚度 [mm]
Area = D * W0 * 10^(-6); % 横截面积 [m²]

linesToSkip = 8; % 表头行数

% 创建图像保存目录
saveDir = fullfile(directory, 'Stress-strain curve');
if ~exist(saveDir, 'dir')
    mkdir(saveDir);
end

% 遍历处理所有TXT文件
set(0, 'DefaultFigureVisible', 'on');
for i = 1:length(files)
    filename = files(i).name;
    
    if endsWith(filename, '.txt')
        fileID = fopen(fullfile(directory, filename));
        if fileID == -1
            warning(['无法打开文件 ', filename, '，跳过']);
            continue;
        end
        
        % 读取数据并跳过表头
        data = textscan(fileID, '%f %f %f', 'HeaderLines', linesToSkip);
        fclose(fileID);
        
        % 提取目标列
        crossheadColumn = data{columnIndices(1)};
        loadColumn = data{columnIndices(2)};
        
        % 检查数据有效性
        if length(crossheadColumn) < 2 || length(loadColumn) < 2
            warning(['文件 ', filename, ' 数据无效，跳过']);
            continue;
        end
        
        % 计算应力应变
        Stress = loadColumn ./ Area;
        L0 = crossheadColumn(1);
        Strain_length = (crossheadColumn - L0) / L0;   
        
        % 弹性阶段拟合计算杨氏模量
        elasticIdx = Strain_length <= 0.01;
        if sum(elasticIdx) < 2
            warning(['文件 ', filename, ' 无足够弹性阶段数据']);
            continue;
        end
        p = polyfit(Strain_length(elasticIdx), Stress(elasticIdx), 1);
        f = polyval(p, Strain_length);
        E = p(1);
        all_E = [all_E; E];
        
        % 绘制应力应变曲线
        fig = figure;
        plot(Strain_length, Stress, 'b', 'LineWidth', 1.2)
        hold on 
        plot(Strain_length, f, 'r--', 'LineWidth', 1.2)
        xlabel('应变')
        ylabel('应力 [MPa]')
        title(['应力-应变曲线: ', filename])
        legend('试验数据','弹性拟合','Location','best')
        grid on;
        hold off
        
        % 保存图像
        saveas(fig, fullfile(saveDir, [filename(1:end-4) '.png']));
        close(fig);
    end
end

% 输出所有杨氏模量
disp('所有试样杨氏模量汇总：');
disp(all_E);

内容的提问来源于stack exchange，提问作者LaLuna Kon