هفتِ بیجار

[coronaa]

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

1465752304.jpg

خطای قبلی رو هم برطرف کردم و کد جدید بصورت زیر شد:

clc;

clear;

close all;





alpha = 0.5;

a = 0;

k = 1;

t = 10;

h = 0.1;





upper_band = round((t-a)/h);

expr = (1./(h.^alpha))*symsum(((-1).^(k))*nchoosek(alpha,k)*f(t-k*h), k , 0 ,upper_band);









f(t) = limit(expr,h,0)

و الآن خطای جدید روی f هست.

Undefined function 'f' for input arguments of type 'double'.





Error in letnikov (line 12)

expr = (1./(h.^alpha))*symsum(((-1).^(k))*nchoosek(alpha,k)*f(t-k*h), k , 0

,upper_band);

مهندس سلام.
میشه روی این معادله هم یه راهنمایی بکنید؟

[mb1372]

مهندس سلام.
میشه روی این معادله هم یه راهنمایی بکنید؟

سلام
تابع زیر ظاهرا بدون مشکل کار میکنه. فقط باید دامنه و ثابت ها رو درست انتخاب کرده باشید.

function y=fun_fun(t)

alpha = 101;

a = 0;

k = 1;

h = 0.1;

y=0;

f=@(x) sin(x)

for k=0:round((t-a)/h)

    y=y+(1./(h.^alpha))*((-1).^(k))*nchoosek(alpha,k)*f(t-k*h);

end

end

[coronaa]

سلام
تابع زیر ظاهرا بدون مشکل کار میکنه. فقط باید دامنه و ثابت ها رو درست انتخاب کرده باشید.

function y=fun_fun(t)

alpha = 101;

a = 0;

k = 1;

h = 0.1;

y=0;

f=@(x) sin(x)

for k=0:round((t-a)/h)

    y=y+(1./(h.^alpha))*((-1).^(k))*nchoosek(alpha,k)*f(t-k*h);

end

end

آلفا نمیشه اعشاری باشه؟
این معادله برای مبحث مرتبه کسری هست و آلفا معمولا یک عدد کسری و بین 0 تا 2 انتخاب میشه.

[mb1372]

آلفا نمیشه اعشاری باشه؟
این معادله برای مبحث مرتبه کسری هست و آلفا معمولا یک عدد کسری و بین 0 تا 2 انتخاب میشه.

در مورد ثابت های این تابع نمیدونم. اما میدونم که کد

nchoosek(alpha,k)

فقط زمانی جواب میده که k کوچکتر از آلفا باشه (طبق تعریف ترکیب). اون مقدار از آلفا رو صرفا برای رهایی از خطا انتخاب کردم.
مقدار ثابت ها رو از یک مرجع پیدا کنید.