Здравствуйте, я хотел бы знать, есть ли способ определить параметры, если я указал входные данные для модели и желаемый результат. Чтобы быть более конкретным, я моделирую модель Windkessel. Это электрическая модель из двух резисторов, одного индуктора и одного конденсатора. Требуемый результат - пульсовая волна.
Модель:
И это желаемый выход модели (напряжение на резисторе R):
Кто-нибудь может мне помочь? Спасибо
ПРИМЕЧАНИЕ.
Matlabsolutions.com предоставляет последнюю Помощь по домашним заданиям MatLab, Помощь по заданию MatLab для студентов, инженеров и исследователей в различных отраслях, таких как ECE, EEE, CSE, Mechanical, Civil со 100% выходом. Код Matlab для BE, B.Tech ,ME,M.Tech, к.т.н. Ученые со 100% конфиденциальностью гарантированы. Получите проекты MATLAB с исходным кодом для обучения и исследований.
Если вы предоставляете ввод и вывод в виде файла «.mat» и вашей модели, должна быть возможность получить значения компонентов. (Если хотите, я могу вывести для него модель схемы.)
Вы не выполняете идентификацию системы, поскольку у вас уже есть модель для вашей системы. Вы делаете оценку параметров вашей модели.
РЕДАКТИРОВАТЬ -
Модель схемы, которую я получил для вашей модели (я позволил Symbolic Math Toolbox выполнить алгебру):
syms C L R r Pi Po P2 s t
Z1 = 1/r;
Z2 = s*C + 1/R;
Z3 = 1/(s*L);
Zt = Z1 + Z2 + Z3;
I = Pi/Zt;
Po = Z2*I;
Out = ilaplace(partfrac(Po, s));
Out = simplify(Out, 'Steps',20);
Out_fcn = matlabFunction(Out, 'Vars',{t,Pi,r,C,R,L})
Out_fcn = @(t,Pi,r,C,R,L) Pi.*dirac(t)-(Pi.*exp((t.*(R+r).*(-1.0./2.0))./(C.*R.*r)).*(cosh((1.0./sqrt(L).*t.*sqrt(L.*R.^2+L.*r.^2-C.*R.^2.*r.^2.*4.0+L.*R.*r.*2.0).*(1.0./2.0))./(C.*R.*r))-1.0./sqrt(L).*sinh((1.0./sqrt(L).*t.*sqrt(L.*R.^2+L.*r.^2-C.*R.^2.*r.^2.*4.0+L.*R.*r.*2.0).*(1.0./2.0))./(C.*R.*r)).*(L.*r+L.*R-C.*R.*r.^2.*2.0).*1.0./sqrt(L.*R.^2+L.*r.^2-C.*R.^2.*r.^2.*4.0+L.*R.*r.*2.0)))./(C.*r);
Значения «Z» — это три сосредоточенных импеданса (и «Zt» — общий импеданс), где выходной сигнал берется через «Z2».
Идея проста: текущее «I» (или кровоток здесь) через серию импедансов — это просто «Pi/Zt». Затем легко рассчитать падение напряжения (или давления) на «Z2» как «Z2*I».
Чтобы подогнать вашу функцию и оценить параметры, нам нужны три вектора данных: «t», «Pi» и «Po», время, входное давление по времени и выходное давление по времени. Функция сопоставляет выходное давление, которое она оценивает с определенными параметрами, с вектором Po.
СМОТРИТЕ ПОЛНЫЙ ОТВЕТ НАЖМИТЕ НА ССЫЛКУ
