【Matlab】基于Prandtl−Ishlinskii的迟滞模型-RLS辨识

PI模型

PI迟滞模型的输出公式:

代码记录

此代码为根据PI模型,已知输入(正弦函数幅值为3.5)、阈值以及权重值,利用matlab生成迟滞回线。

复制代码
%% The Prandtl-Ishlinskii Hysteresis Model-
%% 20241021
clc;clear;
close all;
%% PI Model
%% Define the Input-读取Excel文件进行提取数据列
filename = 'C:\Users\jia\Desktop\yadian\data\test.xlsx';
opts = detectImportOptions(filename);
opts.VariableNamingRule = 'preserve';  % 设置保留原始列名
data = readtable(filename, opts);
t = data{:, 'Time'}; 
V_val = data{:, 'Value1'};
% L_val = data{:, 'Value2'};
v = V_val'; % 输入信号
%% Thresholds and weights
r=0:0.5:3;
n = length(r); % 阈值个数
P=[0.72    0.6    0.48    0.36    0.24    0.12    0.012]; %weights
%% Play operator Calculation of Fr
Fr=zeros(n,N);
for j=1:1:n
    M=max(v(1)-r(j),min(v(1)+r(j),0));
    for k=2:1:N  
        Fr(j,1) = M; 
        F_inc=v(k);
        F_dec=v(k);
        A=F_inc-r(j);
        B=min(F_dec+r(j),Fr(j,k-1));
        Fr(j,k)=max(A,B); % Equation (1)
    end%% End sample
end %% End threshold loop

%% The output of PI model
y=P*Fr;  
% y(1) = 0;
% y_error=y-L_val';

% H1-真实输入输出与估计输出的比较
h1=figure;
set(h1,'WindowStyle','docked');
hold on
plot(1:N, v, 'b', 'LineWidth', 1.5);  
plot(1:N, y, 'r', 'LineWidth', 1.5); 
grid on; 
xlabel('t','Interpreter','latex')
ylabel('y','Interpreter','latex')
legend('V','real-y');
% H2-绘制真实输出和估计对比
h2=figure;
set(h2,'WindowStyle','docked');
hold on;
plot(v,y,'r','linewidth',1.5);
xlabel('V');ylabel('y');
legend('real-y');

运行结果

% H1-真实输入输出随步数变化曲线

% H2-绘制输入和输出迟滞回线

相关推荐
workflower5 分钟前
活动图描述场景
开发语言·软件工程·需求分析·软件需求·敏捷流程
梦想的初衷~7 分钟前
基于现代R语言【Tidyverse、Tidymodel】的机器学习方法
开发语言·机器学习·r语言
香蕉可乐荷包蛋11 分钟前
Python学习之路(十三)-常用函数的使用,及优化
开发语言·python·学习
惜.己19 分钟前
使用python的读取xml文件,简单的处理成元组数组
xml·开发语言·python·测试工具
apihz42 分钟前
域名WHOIS信息查询免费API使用指南
android·开发语言·数据库·网络协议·tcp/ip
coding随想1 小时前
掌控网页的魔法之书:JavaScript DOM的奇幻之旅
开发语言·javascript·ecmascript
爱吃烤鸡翅的酸菜鱼1 小时前
IDEA高效开发:Database Navigator插件安装与核心使用指南
java·开发语言·数据库·编辑器·intellij-idea·database
心情好的小球藻2 小时前
Python应用进阶DAY9--类型注解Type Hinting
开发语言·python
惜.己2 小时前
使用python读取json数据,简单的处理成元组数组
开发语言·python·测试工具·json
Y4090012 小时前
C语言转Java语言,相同与相异之处
java·c语言·开发语言·笔记