هفتِ بیجار

[rahnema1]

ببخشید یه سوال دیگه هم دارم:
روابطی که پایین گذاشتم، بدست آورن انتگرال خطای مطلق یا همون Integral of Absolute Error (IAE) یک سیستم برای ارزیابی نتایج هست :

IAE.jpg

من یه قسمتی از این روابط رو با کد نوشتم. نمیدونم درسته اینطوری یا نه!
حالا نمیدونم چطوری باید وزنگذاری بکنم.

t=0:0.5:200;

out1=rand(1,100);

out2=rand(1,100);

e1=.1*out1;

e2=.1*out2;

f1=int(abs(e1),0,t)+int(abs(e2),0,t);

f2=int(abs(out1(t)-out1(t-1)),0,t)+int(abs(out2(t)-out2(t-1)),0,t);

دستور int برای انتگرال سمبولیک روی یک تابع هست در صورتی که شما اینجا یک بردار دارید که کافیه از دستور sum استفاده کنید
باید اون مقاله یا مقاله های مرتیط را بهدقت مطالعه کنید تا معلوم بشه وزن ها چه طوریه

[coronaa]

دستور int برای انتگرال سمبولیک روی یک تابع هست در صورتی که شما اینجا یک بردار دارید که کافیه از دستور sum استفاده کنید
باید اون مقاله یا مقاله های مرتیط را بهدقت مطالعه کنید تا معلوم بشه وزن ها چه طوریه

خب اگه از دستور sum استفاده کنیم که فقط یک عدد اسکالر داریم!
من باید باید انتگرال رو در هر زمان جدا گانه حساب کنم تا یک بردار داشته باشم. مگه اینطوری نیست؟
چون میخوام مثل همون رگرسیون یه نمودار رسم کنم.
در مورد وزن هم چیزی صحبت نکرده. احتمالاً مانند ماتریس هست که باید برای هر درایه ارزش گذاری کنیم.
نظر شما چیه؟

[rahnema1]

خب اگه از دستور sum استفاده کنیم که فقط یک عدد اسکالر داریم!
من باید باید انتگرال رو در هر زمان جدا گانه حساب کنم تا یک بردار داشته باشم. مگه اینطوری نیست؟
چون میخوام مثل همون رگرسیون یه نمودار رسم کنم.
در مورد وزن هم چیزی صحبت نکرده. احتمالاً مانند ماتریس هست که باید برای هر درایه ارزش گذاری کنیم.
نظر شما چیه؟

مگه قرار نیست انتگرال در طول زمان گرفته بشه؟ پس یک عدد به دست میاد

[coronaa]

مگه قرار نیست انتگرال در طول زمان گرفته بشه؟ پس یک عدد به دست میاد

بله درسته. فکر کنم اشتباه کردم!
ولی باید مقداری بین صفر تا 1 باشه.
حالا این کدی که میذارم کجاش اشتباهه؟

t=0:0.5:200;

out1=rand(1,100);

out2=rand(1,100);

e1=.1*out1;

e2=.1*out2;

f1=sum(abs(e1))+sum(abs(e2));

f2=sum(abs(out1(t+1)-out1(t)))+sum(abs(out2(t+1)-out2(t)));

[coronaa]

بله درسته. فکر کنم اشتباه کردم!
ولی باید مقداری بین صفر تا 1 باشه.
حالا این کدی که میذارم کجاش اشتباهه؟

t=0:0.5:200;

out1=rand(1,100);

out2=rand(1,100);

e1=.1*out1;

e2=.1*out2;

f1=sum(abs(e1))+sum(abs(e2));

f2=sum(abs(out1(t+1)-out1(t)))+sum(abs(out2(t+1)-out2(t)));

مهندس ببخشید.
تا f1 فکر کنم درست نوشتم. الآن روی f2 خطا میده. میشه بگید مشکلش چیه؟
فعلاً به وزنها کاری ندارم.

[rahnema1]

باید اینجور باشه قاعدتا

f2 = sum(abs(diff(out1))) + sum(abs(diff(out2)));

[coronaa]

باید اینجور باشه قاعدتا

f2 = sum(abs(diff(out1))) + sum(abs(diff(out2)));

وقتی delat T داره دیگه زمان دخیل نیست؟

[rahnema1]

شاید لازم باشه اون جمع ضربدر واحد اپسیلون زمان بشه که این هم تنها یک ضریبه که نهایتا در فرمول بهینه سازی تاثیری نداره

[coronaa]

سلام.
کسی میتونه راهنمایی کنه بگه این معادله دیفرانسیل چطوری حل میشه؟



کدش رو بصورت زیر نوشتم ولی خطا میده!

clc

syms T(t)

syms R

syms U

M=1 ; Cp=2 ; I=3 ; h=4 ; A=5; Tinf=25 ;





% dsolve('DT=(1/M.Cp).*((I^2*R+I*T*U)-(h*A(T-Tinf)))','t')





o = diff(T) == (1/(M.Cp)).*(((I.^2*R)+(I*T*U))-(h*A(T-Tinf)));

T(t) = dsolve(o,'T(0)=1')

خطا :

Dot indexing is not supported for variables of this type.


Error in diffffff (line 11)
o = diff(T) == (1/(M.Cp)).*(((I.^2*R)+(I*T*U))-(h*A(T-Tinf)));

[mb1372]

سلام.
کسی میتونه راهنمایی کنه بگه این معادله دیفرانسیل چطوری حل میشه؟



کدش رو بصورت زیر نوشتم ولی خطا میده!

clc

syms T(t)

syms R

syms U

M=1 ; Cp=2 ; I=3 ; h=4 ; A=5; Tinf=25 ;





% dsolve('DT=(1/M.Cp).*((I^2*R+I*T*U)-(h*A(T-Tinf)))','t')





o = diff(T) == (1/(M.Cp)).*(((I.^2*R)+(I*T*U))-(h*A(T-Tinf)));

T(t) = dsolve(o,'T(0)=1')

خطا :

سلام

ببینید جواب درست میده؟

clc

syms T(t)

syms R

syms U

M=1 ; Cp=2 ; I=3 ; h=4 ; A=5; Tinf=25 ;

 

% dsolve('DT=(1/M.Cp).*((I^2*R+I*T*U)-(h*A(T-Tinf)))','t')

 

o = diff(T) == (1/(M.*Cp)).*(((I.^2*R)+(I*T*U))-(h*A.*(T-Tinf)));

T(t) = dsolve(o,'T(0)=1')

[coronaa]

سلام

ببینید جواب درست میده؟

clc

syms T(t)

syms R

syms U

M=1 ; Cp=2 ; I=3 ; h=4 ; A=5; Tinf=25 ;

 

% dsolve('DT=(1/M.Cp).*((I^2*R+I*T*U)-(h*A(T-Tinf)))','t')

 

o = diff(T) == (1/(M.*Cp)).*(((I.^2*R)+(I*T*U))-(h*A.*(T-Tinf)));

T(t) = dsolve(o,'T(0)=1')

سلام.
بله جواب داد :)
تشکر