p4判断一个5位数是不是回文数

%% 主程序 % 1. 数据读取 %% 第一部分：数据读取和预处理 clear; clc; close all; % 添加lvm_import函数路径 addpath('C:\Users\kaixiang\Desktop\实验数据1'); % 读取数据文件 data_struct = lvm_import('D:\Fluent Macrophone 2\FL5\O100\Methane0.9\Forward\EQ0.9\PMT_0.lvm'); data_struct1 = lvm_import('D:\Fluent Macrophone 2\FL5\O0\Methane0.9\Forward\EQ0.9\PMT_1029.lvm'); % 提取数据列 p3 = data_struct.Segment1.data(:, 4); % 压力信号 p4 = data_struct1.Segment1.data(:, 1); % 参考信号 p5 = max(abs(p4)); data = (p3 / p5) * 10; % 元素级除法 fs = 50e3; t = (0:length(data)-1)/fs; % 时间轴 % 时域波形 figure; subplot(2,1,1); plot(t, data); xlabel('时间(s)'); ylabel('幅值'); title('时域波形[^2]'); % 时域特征 rms_val = sqrt(mean(data.^2)); cf = max(abs(data))/rms_val; disp(['RMS: ', num2str(rms_val), ' CF: ', num2str(cf)]); N = length(data); f = (0:N-1)*(fs/N); % 频率轴 X = fft(data)/N; % FFT计算 % 幅度谱 subplot(2,1,2); plot(f(1:N/2), abs(X(1:N/2))); xlabel('频率(Hz)'); ylabel('幅度'); title('频域分析[^1][^3]'); % 主频检测 [~, idx] = max(abs(X(1:N/2))); disp(['主频: ', num2str(f(idx)), ' Hz']); % 4. 相空间重建 tau = AMI_Method(data, 50); % AMI求延迟时间 m = CAO_Method(data, tau, 10); % CAO求嵌入维度 Y = phase_space_recon(data, m, tau); % 5. 递归图与量化 RP = RecurrencePlot(Y, 0.1*std(data)); [RR, DET] = RQA(RP); disp(['递归率: ', num2str(RR), ' 确定性: ', num2str(DET)]); function tau = AMI_Method(data, max_tau) mi = zeros(1, max_tau); for tau = 1:max_tau x = data(1:end-tau); y = data(tau+1:end); mi(tau) = mutual_info(x, y); % 需自定义互信息计算函数 end [~, tau] = findpeaks(-mi, 'NPeaks', 1); end function m = CAO_Method(data, tau, max_m) E = zeros(1, max_m); for m = 1:max_m % 重构相空间 Y = phase_space_recon(data, m, tau); % 计算统计量E1 E(m) = compute_E1(Y); % 需自定义计算函数 end % 寻找E1变化率小于阈值的最小m m = find(diff(E) < 0.01, 1) + 1; end function RP = RecurrencePlot(Y, epsilon) N = size(Y,1); RP = zeros(N); for i = 1:N for j = 1:N RP(i,j) = norm(Y(i,:)-Y(j,:)) < epsilon; end end imag