تبلیغات
مهندسی مکانیک - سری فوریه در متلب

سری فوریه در متلب

نویسنده :مسعود شمس
تاریخ:جمعه 6 آبان 1390-08:20 ب.ظ

(این اموزش به درخواست شما دوستان در وبلاگ قرار گرفته است)
.می خواهیم سری فوریه تابع شکل زیر را به کمک متلب با محاسبه 3 ,7, 13 جمله از بدست اورده و رسم کنیم


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

را محاسبه می کنیم: bnبرابر صفر می باشد.با استفاده از نرم افزار متلب ضریب an , a0 تابع فرد می باشد پس ضرایب

>> syms  x  n

>>bn=int((-10*sin(n*pi*x)),-1,0)+int((10*sin(n*pi*x)),0,1);

   bn =

    (20*sin((pi*n)/2)^2)/(pi*n) - (10*(cos(pi*n) - 1))/(pi*n)

 

با استفاده از اسکریپت زیر تابع مورد نظر را رسم می کنیم:


t=-1:0.01:1;
f=zeros(1,length(t));
y=0;
for i=1:length(t)
    if t(i)<0
        f(i)=-10;
    else
        f(i)=10;
    end
end
f=[f f f];
t=linspace(-3,3,length(f));
plot(t,f)
axis([-3 3 -12 12])
xlabel('time')
ylabel('f(t)')

که نتیجه این کد همان شکل فوق می شود
از معادله حاصل از انتگرال گیری و همچنین برنامه نوشته شده فوق با استفاده از دستور Hold on استفاده کرده محاسبه سه جمله از سری را بدست می اوریم:
for n=1:3
    an=0;
    bn=(20*sin((pi*n)/2)^2)/(pi*n) - (10*(cos(pi*n) - 1))/(pi*n);
    y=y+an*cos(n*t)+bn*sin(n*t);
end
y=[y y y];
t=linspace(-3,3,length(y));
plot(t,y,'r')
legend('f(t)','n=3')
title('f(t) fouriers series')
xlabel('time')
ylabel('f(t)')


متغیر n در حلقه نمایان گر تعداد سری محاسبه می باشد با تغییر مقذار n به 3 ,7 , 13 نتایج زیر حاصل می شود.



اسکریپت برای نمایش تمامی جمله ها در یک نمودار :
nmax=[3 7 13];
tspan=-1:0.01:1;
f=zeros(1,length(tspan));
for i=1:length(tspan)
    if tspan(i)<0
        f(i)=-10;
    else
        f(i)=10;
    end
end
f=repmat(f,1,3);
t=linspace(-3,3,length(f));
plot(t,f)
hold on
title('function f(t)')
xlabel('time')
ylabel('f(t)')
axis([-3 3 -15 15])
con=1;
color={'r' 'b' 'g'};
for i=nmax
    fun=0;
   for n=1:i
       an=0;
       bn=(20*sin((pi*n)/2)^2)/(pi*n) - (10*(cos(pi*n) - 1))/(pi*n);
       fun=fun+an*cos(n*t)+bn*sin(n*t);
   end
   fun=repmat(fun,1,3);
   tspan=linspace(-3,3,length(fun));
   plot(tspan,fun,color{con})
   spin=sprintf('fouriers series n=%g',i);
   title(spin)
   xlabel('time')
   ylabel('f(t)')
   hold on
   grid on
   con=con+1;
end
legend('f(t)','n=3','n=7','n=13')



نوع مطلب : اموزش متلب 

داغ کن - کلوب دات کام
نظرات() 
viabiovit.com/generieke-viagra-britse-online.html
سه شنبه 23 مرداد 1397 02:32 ق.ظ

Nicely put. Kudos!
where can i buy the cheapest viagra sildenafil citrate viagra no presc buy viagra sample can you really buy viagra online viagra online buying viagra without prescription viagra online pharmacy viagra buy viagra buy viagra online no prescription buy generic viagra online canada
Cialis canada
دوشنبه 7 خرداد 1397 01:33 ق.ظ

Nicely put, Cheers.
cialis 200 dollar savings card sublingual cialis online achat cialis en suisse cialis 5 mg scheda tecnica cialis 5 mg scheda tecnica cialis dosage cialis online nederland generic cialis soft gels low dose cialis blood pressure we recommend cialis best buy
Cialis canada
پنجشنبه 20 اردیبهشت 1397 02:39 ق.ظ

Effectively voiced indeed. .
generic cialis pill online trusted tabled cialis softabs cialis 5 mg funziona cialis tablets for sale cialis for sale south africa rx cialis para comprar we recommend cheapest cialis cialis herbs cialis usa cost cialis 10mg prix pharmaci
Cialis generic
جمعه 7 اردیبهشت 1397 06:05 ب.ظ

You made the point!
safe dosage for cialis cost of cialis per pill fast cialis online viagra vs cialis vs levitra tadalafil 5mg 40 mg cialis what if i take canadian cialis price cialis wal mart pharmacy cialis online holland buy cialis sample pack
Viagra lowest price
سه شنبه 4 اردیبهشت 1397 01:28 ق.ظ

You have made the point!
rx viagra online buy viagra generic how to buy viagra without prescription can i buy viagra over the counter sildenafil buy online online viagra no prescription how old do you have to be to buy viagra where can i purchase viagra buying viagra online australia buy online viagra uk
Cialis canada
جمعه 17 فروردین 1397 05:36 ب.ظ

Many thanks. Plenty of stuff.

online cialis cialis authentique suisse order generic cialis online cialis generico online cialis lilly tadalafi cialis manufacturer coupon cialis 5 mg schweiz cialis mit grapefruitsaft overnight cialis tadalafil cialis flussig
Online cialis
جمعه 3 فروردین 1397 12:58 ب.ظ

Nicely put, Kudos!
cialis prices generic cialis 20mg tablets we use it cialis online store cialis therapie cialis for bph cialis 5 mg funziona best generic drugs cialis we choice free trial of cialis cialis purchasing cialis savings card
Cialislaf
یکشنبه 27 اسفند 1396 08:19 ب.ظ

Nicely put, Regards!
tadalafil 20mg interactions for cialis buying cialis overnight pastillas cialis y alcoho opinioni cialis generico tadalafil tablets side effects of cialis preis cialis 20mg schweiz cialis sicuro in linea cialis generico lilly
BHW
پنجشنبه 24 فروردین 1396 07:57 ب.ظ
Heya i'm for the first time here. I found this board
and I find It truly useful & it helped me out much. I hope to present something back
and help others such as you helped me.
manicure
سه شنبه 15 فروردین 1396 08:09 ق.ظ
I delight in, cause I found exactly what I was having a look
for. You have ended my 4 day long hunt! God Bless you man.
Have a nice day. Bye
 
لبخندناراحتچشمک
نیشخندبغلسوال
قلبخجالتزبان
ماچتعجبعصبانی
عینکشیطانگریه
خندهقهقههخداحافظ
سبزقهرهورا
دستگلتفکر