یک روش تشخیص شکست و ایزوله سازی مکانیکی با استفاده از فیلترینگ کالمن برای بهبود امنیت سیستمهای فیزیکی سایبری

---

شرح پروژه

با سلام یک پروژه کلاسی دارم که شبیه سازی آن انجام شده است. اما در ران شدن خطا می دهد. تقاضا دارم که شیه سازی آن به نحوی اصلاح شوذ که خروجی آن شامل تمام نمودارهای موجود در مقاله ای که استادم به من داده باشد. یک مقاله بیس دارم (به زبان انگلیسی). ترجمه آن را نیز دارم. قسمت سیمولینک مقاله مورد نیاز نیست. خروجی ها باید کاملا مطابق با شکل 1 تا 21 مقاله باشند. کلیه اطلاعات، فرمول ها، داده های ورودی نیز در متن مقاله آمده است.

فایل متلب را هم دارم که متاسفانه نمی توان آن را ضمیمه کرد ولی همه کدهارا در اینجا  کپی و ارسال می کنم. لازم به ذکر است که مقاله 7 سناریو دارد که همگی در کد نوشته و اجرا شده اند. در صورت انجام و تایید توسط استاد راهنما، پایان نامه ام را نیز به شما می دهم. با تشکر.

clc; clear all; close all; warning off;

ARW dhRW/ dt= QPUB +Qpermeate +Qrecycle ?QRWout;

RW(t)= 1/ ARW [1 1 1-1] {QPUB, Qpermeate, Qrecycle, QRWout},

Initial_tolerance=0.001;

x(k+1)=Ax(k)+Bu(k));

y(k)=Cx(k)+0.u(k);

x = h;

u = [QPUBQpermeateQrecycleQRWout]T;

ARW= 1.4957;

A = 1;

B=[0.6686 0.6686 0.6686 ?0.6686], C = 1,

hRW,o= 0.6,

QPUB= Qpermeate= Qrecycle= QRWout= 0,

APre-UF (dhPre-UF/dt)=QRWout +QNacl +QHcl +QNaocl- QPre-UFout,

Pre-UF (t)= 1/APre-UF [1 1 1 1-1]{QRWout, QNacl,QHcl,QNaocl,QPre-UFout};

x(k+1)=Ax(k)+Bu(k);

y(k)=Cx(k)+0.u(k);

x = h;

u = [QRWoutQNaclQHclQNaoclQPre-UFout]T;

Output_Bus=zeros(Bus_Number);

APre-UF= 1.4957;

A = 1, B =[0.6686 0.6686 0.6686 0.6686 ?0.6686], C = 1;

hPre-UF,o= 0.9;

QRWout= QNacl= QHcl= QNaocl= QPre-UFout= 0;

xk+1 =Axk +Buk;

yk =Cxk;

zk =xsp+a(k-t)+sk;

zk =xsp,new +sk;

xk=Axk-1 +Buk;

P-k=APk-1AT +Q;

rk =zk-Cx-k;

rcum,k =rcum,k-1 +rk;

u1,k =GMV101 ·c1,k;

io1,k =u1 ,k +_k;

c1,k =u1,k =io1,k;

c1 ,k =io1 ,k ±3;

if zk(LIT101) > 0.8

    iok(FIT101) = 0;

   end

if zk(LIT101) < 0.5

    iok(FIT101) > 0;

end

   if zk(LIT101) < 0.25

    iok(FIT201) = 0;

   end

   K=P-k CT(CP-k CT +R)-1;

   xk =x-k+Krk;

   Pk =(I-KC)P-k;

   if zk(FIT101) <0

       end

   if zk(LIT101) <0

       end

   if zk(C1101, RW) <0

      end

   if zk(FIT1101, RW) <0

      end

   if zk(LIT101, MV-101) <0

      end

   if zk(FIT1101, LIT101) <0

       end

   if zk(FIT201, LIT301, UF) <0

    end

result=PLOTS;

end

این پروژه شامل 2 فایل مهم است، لطفا قبل از ارسال پیشنهاد حتما نسبت به بررسی این فایل اقدام فرمایید.

مهارت ها و تخصص های مورد نیاز

---

---