Kod:
[ %%cosinüs fonksiyonunun sembolik ve numerik türevini aynı %%grafik üzerinde çizdiren matlab programıdır. syms t Fs Ac%değişkenleri sembolik olarak tanımladık st=Ac*cos(2*pi*Fs*t);%st fonksiyonunu tanımladık stp=diff(st,t);%st fonksiyonunun t'ye göre türevini stp'ye atadık Fs =1;%frekansa kolaylık sağlaması için 1 değerini verdik t=[0:0.01:1];%t değişkeninde 1x101 lik bir vektör oluşturduk Ac=1;%%kolaylık olması için Ac değişkenine 1 değerini verdik st=eval(st);stp=eval(stp);%sembolik olarak tanımlanan st ve stp fonksiyonlarını eval(x) komutu ile numeriğe çevirdik %%numerik türev alımı%% sntp=(st(1:100)-st(2:101))./(t(1:100)-t(2:101));%burada türev alırken semboller yerine sayıları kullandık %ve çok küçük aralıklarla cosinüs grafiğinin eğimini hesapladık(limit gibi düşünebiliriz) sntp= [sntp -0.1966];%yukardaki işlem sonucunda 100 değer ortaya çıktı %ama bizim tanımladığımız t aralığında 101 eleman olduğundan %grafiği çizdirebilmek için bir değer ekledik bu değer eklemesini de %sntp'nin elemanlarındaki 99. elemandan 100. elemanı %çıkartıp gelen sonucu da 100. elemannın yanına aktararak %örüntüyü bozmayacak biçimde sntp'nin değerlerine +1 %eklemesi yaptık figure(1)%yeni figure penceresi açtık ama çok da önemli değil yazılamasa da kod sorunsuz çalışır plot(t,st,':or');hold on;%t aralığı üzerinde st değerlerini çizdirdik bunu yaparken grafiğin geçtiği noktalara o işareti,grafiğin rengine kırmızı,çizgisine noktalı çizgi ataması yaptık plot(t,stp,'-*b');hold on;%t aralığı üzerinde stp değerlerini çizdirdik bunu yaparken grafiğin geçtiği noktalara * işareti,grafiğin rengine mavi,çizgisine düz çizgi ataması yaptık plot(t,sntp,'+--g');%t aralığı üzerinde sntp değerlerini çizdirdik bunu yaparken grafiğin geçtiği noktalara + işareti,grafiğin rengine yeşil,çizgisine kesikli çizgi ataması yaptık xlabel('time');%x ekseni ismi...:time ylabel('st,sntp and stp function');%y ekseni ismi...: st,sntp and stp function title('cos differantial');%başlık....:cos differantial grid on;%grafik pencere ızgarasını açtık %%Created by tozgen
yukarıdaki kodda t aralığını küçültüp ve diğer ayarlamaları da ona göre değiştirdiğimizde
numerik türevin doğruluğu teorik türevden uzaklaşacaktır ve biz aralığı sonsuz yapmadığımız müddetçe de asla teorik türeve eşit olmayacaktır.