ileri mekanik titreşimler ikinci ödevi

[muslum yaman]

arkadaşlar

tan(beta22)-1/(16*alfak*beta22*(eta)kare*(1-eta)kare)=0
dekleminin analitik çözümü varmı?Bulan varmı?

[Can Cengiz]

Selamlar. Onun icin en pratik yol grafiklerin kesistigi noktayi
bulmak. Ya da iki ayri fonksiyonun degerleri arasindaki farkin belli
bir toleranstan dusuk oldugu bir degeri while dongusu icinde
arayabilirsin. Ama analitik yol hiç pratik degil. Kolay gelsin
herkese

--
Can Cengiz

[Deniz Yazgaç]

Analitik yol zaten yok.

Can'ın dediği gibi; bir tolerans dahilinde while döngüsü olabilir, ya da her
bir değişken basamağı için grafikler çizdirilir ve değerler grafikten
bakılarak yaklaşık alınır.

[FAT OZ]

yazan biri var mı o kodu:)) buln biri bana da verse

---

From: Deniz Yazgaç <deniz...@gmail.com>
To: itu_...@googlegroups.com
Sent: Sunday, May 3, 2009 5:34:49 PM

Subject: [itu_mdta] Re: ileri mekanik titreşimler ikinci ödevi

[Deniz Yazgaç]

Merhaba Arkadaşlar;
Laura'nın makalesindeki formüllerle; Gürgöze'in makalesindeki formüllere nasıl geçildiğini anlatan bir çalışma ekte



Deniz Yazgaç
YL Öğrencisi
deniz...@gmail.com

03 Mayıs 2009 Pazar 18:16 tarihinde FAT OZ <fatih_...@yahoo.com> yazdı:

[Deniz Yazgaç]

Cihan'ın yaptığı çalışma
w11 nasıl bulunur sorusunu, çalışmasıyla cevaplandırmış.

Bu arada Tuncer Toprak'ın 200 bilmem kaç problem kitabında da , sistemin nasıl modellendiği yazıyor ancak biraz karışık.

Aşağıdan indirebilirsiniz.

http://www.zshare.net/download/5954406222cacb47/

Deniz Yazgaç
YL Öğrencisi
deniz...@gmail.com

03 Mayıs 2009 Pazar 18:23 tarihinde Deniz Yazgaç <deniz...@gmail.com> yazdı:

[Can Cengiz]

Selamlar,

beta_22'yi ben yaklaşık 1,077 buldum.

Can Cengiz

2009/5/3 Deniz Yazgaç <deniz...@gmail.com>

[Umut Murat Gök]

Can selam,

 

beta_22 alfa_k ya bağlı olduğu için alfa_k nın 8 ayrı değeri için farklı değerler alıyor. Yani beta_22 için 8 elemanlı bir vektör elde etmen gerekiyor. Ona dikkat et istersen.

04 Mayıs 2009 Pazartesi 13:33 tarihinde Can Cengiz <can.c...@gmail.com> yazdı:

--
Umut Murat Gök

[Can Cengiz]

Hmmm. Haklısın. Beta_22'nin ilk kökü diye belirtiyor ama her defasında transandental denklem değişiyor. Teşekkürler, dikkat etmek gereken bir nokta daha çıktı ortaya.

Can Cengiz

2009/5/4 Umut Murat Gök <umut.mu...@gmail.com>

[FAT OZ]

onu buldum ben de ama betanın virgulden sonra 6 basamağını bulmazsak sonuclar tutmuyo. uğraşlar devam ediyor:D

---

From: Can Cengiz <can.c...@gmail.com>
To: itu_...@googlegroups.com
Sent: Monday, May 4, 2009 2:07:51 PM

[Deniz Yazgaç]

%Deniz Yazgaç
%Değişkenleri değiştirmeyi unutmayın bire bir aynısı olmasın. ek bir takım
%işler de ekleyebilirsiniz.
%Saygılar

clc
clear

ak=0.5:0.1:10; %ak değeri
am=ak; %am değeri ak'nin aynısı
beta_22=0:0.01*pi:0.4*pi; %beta_22 vektörü, 0.4*pi yeterli denklemler her zaman 0.5 pi'den evvelde kesişmek zorunda.Bunun için tanjantın grafiğine bakmanız yeterli(herhangi bir kaynaktan)
beta_11=pi; %birinci alt sistem beta_11
n=0.5 %çubuğun ortası
m_21=0.01; %M_21
for v=1:1:length(am);  %her adımda am vektöründen bir eleman alınacak.
for i=1:1:length(ak); %her adımda ak vektöründen bir eleman alınacak.
for j=1:1:length(beta_22);%her adımda beta_22'den bir eleman alınacak.
   
        denk_1(j)=tan(beta_22(j));%frekans denkleminin bir tarafı
        denk_2(j)=1./(16.*ak(i).*beta_22(j).*(n^2).*((1-n).^2));%frekans denkleminin diğer bir tarafı
        kars=abs(denk_1-denk_2);%eleman-eleman çıkart mutlak değerini al
        [deger,sira]=min(kars);%bu mutlak değerin minimumunu al(ki bu da köktür); o minimumu deger değişkeninde sakla, sırasını da sıra değişkeninde
        B(i)=beta_22(sira);%sıra neyse beta_22'den bu değeri çek ve i değeriyle indeksle, böylece ak ile ilişkili olacak.
end
   
    w(i,v)=((beta_11)^2)./(1+(  ( ((beta_11)^4)./((B(i))^2) ).* ( m_21 ./ (48.*ak(i)) ) )+( am(v).* ((beta_11)^4)) .*(( 1./( 48.*ak(i) ) )+( (n^2).*((1-n).^2) )./3) );
%her i ve v için w1/w0'ı hesapla         
     
end
end
figure(1)
contourf(ak,am,w);%Yukarıdan bakış bunu çiz (w1/w0)
grid on
xlabel('a_k parametresi');
ylabel('a_m parametresi');
title('w_1 / w_0 değerleri (Dunkerley Yaklaşımı)');

figure(2)
mesh(ak,am,w);%yetmez bir de mesh olarak çiz (w1/w0)
grid on
xlabel('a_k parametresi');
ylabel('a_m parametresi');
title('w_1 / w_0 değerleri (Dunkerley Yaklaşımı)');

Deniz Yazgaç
YL Öğrencisi
deniz...@gmail.com

04 Mayıs 2009 Pazartesi 18:44 tarihinde FAT OZ <fatih_...@yahoo.com> yazdı:

[Deniz Yazgaç]

Matlab'de preferance'den long e görüntüsünü seç Fatih.

Deniz Yazgaç
YL Öğrencisi
deniz...@gmail.com

04 Mayıs 2009 Pazartesi 18:53 tarihinde Deniz Yazgaç <deniz...@gmail.com> yazdı:

[Deniz Yazgaç]

File-preferances-command window-Numeric format.

Deniz Yazgaç
YL Öğrencisi
deniz...@gmail.com

04 Mayıs 2009 Pazartesi 18:55 tarihinde Deniz Yazgaç <deniz...@gmail.com> yazdı:

[Deniz Yazgaç]

Formülü Cihan'ın gönderdiği yoldan düzenlediyseniz, MATLAB kodundaki formül değişecektir.

Deniz Yazgaç
YL Öğrencisi
deniz...@gmail.com

04 Mayıs 2009 Pazartesi 18:56 tarihinde Deniz Yazgaç <deniz...@gmail.com> yazdı:

[FAT OZ]

deniz sanki bu programda bir sıkıntı var. hocanın makaledeki beta_22 nin, alfa_k=0.5 için değeri 1.0751 cıkıyo. bı programdan bunu nasıl cıkartacaz, yani anlamadıgım icin soruyorum:)) ben bu değeri buldum ama virgülden sonra 6 hane gerekiyo. 

---

From: Deniz Yazgaç <deniz...@gmail.com>
To: itu_...@googlegroups.com

Sent: Monday, May 4, 2009 6:53:51 PM

[Deniz Yazgaç]

bizden beta_22'yi istemiyor ki. Gönderdiğim kod, beta_22'leri hesaplıyor ancak göstermeden kullanıyor.
İstersen workspace'den B vektörüne bakabilirsin. Beta_22 orada vektörel olarak var.

Deniz Yazgaç
YL Öğrencisi
deniz...@gmail.com

04 Mayıs 2009 Pazartesi 20:11 tarihinde FAT OZ <fatih_...@yahoo.com> yazdı:

[Deniz Yazgaç]

Kodun sonuna şunları da ekleyerek. Tablo olarak değerleri görebilirsiniz.
format long
ist_ak=[0.5 1 1.5 2 2.5 3 5 10];
ist_am=[0.5 1 1.5 2 2.5 3 5 10];
tablo=zeros((length(ist_ak)+1),length(ist_am)+1);

tablo(2:(length(ist_ak)+1),1)=ist_ak;
tablo(1,2:(length(ist_am)+1))=ist_am;

for v=1:1:length(ist_am);  %her adımda ist_am vektöründen bir eleman alınacak.
for i=1:1:length(ist_ak); %her adımda ist_ak vektöründen bir eleman alınacak.

for j=1:1:length(beta_22);%her adımda beta_22'den bir eleman alınacak.
   
        denk_1(j)=tan(beta_22(j));%frekans denkleminin bir tarafı

        denk_2(j)=1./(16.*ist_ak(i).*beta_22(j).*(n^2).*((1-n).^2));%frekans denkleminin diğer bir tarafı

        kars=abs(denk_1-denk_2);%eleman-eleman çıkart mutlak değerini al
        [deger,sira]=min(kars);%bu mutlak değerin minimumunu al(ki bu da köktür); o minimumu deger değişkeninde sakla, sırasını da sıra değişkeninde
        B(i)=beta_22(sira);%sıra neyse beta_22'den bu değeri çek ve i değeriyle indeksle, böylece ak ile ilişkili olacak.
end
   

    w_ist(i,v)=((beta_11)^2)./(1+(  ( ((beta_11)^4)./((B(i))^2) ).* ( m_21 ./ (48.*ist_ak(i)) ) )+( ist_am(v).* ((beta_11)^4)) .*( 1./( 48.*ist_ak(i) ) ).*(1+ist_ak(i)) );
       
      tablo(i+1,v+1)=w_ist(i,v);
end
end
disp(tablo)

Deniz Yazgaç
YL Öğrencisi
deniz...@gmail.com

04 Mayıs 2009 Pazartesi 21:38 tarihinde Deniz Yazgaç <deniz...@gmail.com> yazdı:

[Müslüm Yaman]

gözlerimi yaşarttınız hepinizi tebrik ediyorum :)

On 4 Mayıs, 21:40, Deniz Yazgaç <denizyaz...@gmail.com> wrote:
> Kodun sonuna şunları da ekleyerek. Tablo olarak değerleri görebilirsiniz.

> *format long

> disp(tablo)*
>
> Deniz Yazgaç
> YL Öğrencisi
> denizyaz...@gmail.com

>
> 04 Mayıs 2009 Pazartesi 21:38 tarihinde Deniz Yazgaç

> <denizyaz...@gmail.com>yazdı:

>
>
>
> > bizden beta_22'yi istemiyor ki. Gönderdiğim kod, beta_22'leri hesaplıyor
> > ancak göstermeden kullanıyor.
> > İstersen workspace'den B vektörüne bakabilirsin. Beta_22 orada vektörel
> > olarak var.
>
> > Deniz Yazgaç
> > YL Öğrencisi

> > denizyaz...@gmail.com
>
> > 04 Mayıs 2009 Pazartesi 20:11 tarihinde FAT OZ <fatih_ozba...@yahoo.com>yazdı:

>
> > deniz sanki bu programda bir sıkıntı var. hocanın makaledeki beta_22 nin,
> >> alfa_k=0.5 için değeri 1.0751 cıkıyo. bı programdan bunu nasıl cıkartacaz,
> >> yani anlamadıgım icin soruyorum:)) ben bu değeri buldum ama virgülden sonra
> >> 6 hane gerekiyo.
>

> >> ------------------------------
> >> *From:* Deniz Yazgaç <denizyaz...@gmail.com>
> >> *To:* itu_...@googlegroups.com
> >> *Sent:* Monday, May 4, 2009 6:53:51 PM
>
> >> *Subject:* [itu_mdta] Re: ileri mekanik titreşimler ikinci ödevi

> >> denizyaz...@gmail.com
>
> >> 04 Mayıs 2009 Pazartesi 18:44 tarihinde FAT OZ <fatih_ozba...@yahoo.com>yazdı:

>
> >>> onu buldum ben de ama betanın virgulden sonra 6 basamağını bulmazsak
> >>> sonuclar tutmuyo. uğraşlar devam ediyor:D
>

> >>> ------------------------------
> >>> *From:* Can Cengiz <can.cen...@gmail.com>
> >>> *To:* itu_...@googlegroups.com
> >>> *Sent:* Monday, May 4, 2009 2:07:51 PM
>
> >>> *Subject:* [itu_mdta] Re: ileri mekanik titreşimler ikinci ödevi

>
> >>> Hmmm. Haklısın. Beta_22'nin ilk kökü diye belirtiyor ama her defasında
> >>> transandental denklem değişiyor. Teşekkürler, dikkat etmek gereken bir nokta
> >>> daha çıktı ortaya.
>
> >>> Can Cengiz
>

> >>> 2009/5/4 Umut Murat Gök <umut.murat....@gmail.com>

>
> >>>> Can selam,
>
> >>>> beta_22 alfa_k ya bağlı olduğu için alfa_k nın 8 ayrı değeri için farklı
> >>>> değerler alıyor. Yani beta_22 için 8 elemanlı bir vektör elde etmen
> >>>> gerekiyor. Ona dikkat et istersen.
>
> >>>> 04 Mayıs 2009 Pazartesi 13:33 tarihinde Can Cengiz <

> >>>> can.cen...@gmail.com> yazdı:

>
> >>>>> Selamlar,
>
> >>>>> beta_22'yi ben yaklaşık 1,077 buldum.
>
> >>>>> Can Cengiz
>

> >>>>> 2009/5/3 Deniz Yazgaç <denizyaz...@gmail.com>
>
> >>>>> *Cihan'ın yaptığı çalışma

> >>>>>> w11 nasıl bulunur sorusunu, çalışmasıyla cevaplandırmış.
>
> >>>>>> Bu arada Tuncer Toprak'ın 200 bilmem kaç problem kitabında da ,
> >>>>>> sistemin nasıl modellendiği yazıyor ancak biraz karışık.

> >>>>>> *
> >>>>>> *Aşağıdan indirebilirsiniz.*
>
> >>>>>> *http://www.zshare.net/download/5954406222cacb47/*
>
> >>>>>> Deniz Yazgaç
> >>>>>> YL Öğrencisi
> >>>>>> denizyaz...@gmail.com

>
> >>>>>> 03 Mayıs 2009 Pazar 18:23 tarihinde Deniz Yazgaç <

> >>>>>> denizyaz...@gmail.com> yazdı:

>
> >>>>>> Merhaba Arkadaşlar;
> >>>>>>> Laura'nın makalesindeki formüllerle; Gürgöze'in makalesindeki
> >>>>>>> formüllere nasıl geçildiğini anlatan bir çalışma ekte
>
> >>>>>>> Deniz Yazgaç
> >>>>>>> YL Öğrencisi

> >>>>>>> denizyaz...@gmail.com
>
> >>>>>>> 03 Mayıs 2009 Pazar 18:16 tarihinde FAT OZ <fatih_ozba...@yahoo.com>yazdı:

>
> >>>>>>>> yazan biri var mı o kodu:)) buln biri bana da verse

> >>>>>>>>  ------------------------------
> >>>>>>>> *From:* Deniz Yazgaç <denizyaz...@gmail.com>
> >>>>>>>> *To:* itu_...@googlegroups.com
> >>>>>>>> *Sent:* Sunday, May 3, 2009 5:34:49 PM
> >>>>>>>> *Subject:* [itu_mdta] Re: ileri mekanik titreşimler ikinci ödevi

>
> >>>>>>>> Analitik yol zaten yok.
>
> >>>>>>>> Can'ın dediği gibi; bir tolerans dahilinde while döngüsü olabilir,
> >>>>>>>> ya da her
> >>>>>>>> bir değişken basamağı için grafikler çizdirilir ve değerler
> >>>>>>>> grafikten
> >>>>>>>> bakılarak yaklaşık alınır.
>
> >>>>>>>> -----Original Message-----
> >>>>>>>> From: itu_...@googlegroups.com [mailto:itu_...@googlegroups.com]
> >>>>>>>> On Behalf
> >>>>>>>> Of Can Cengiz
> >>>>>>>> Sent: Sunday, May 03, 2009 4:29 PM
> >>>>>>>> To: itu_...@googlegroups.com
> >>>>>>>> Subject: [itu_mdta] Re: ileri mekanik titreşimler ikinci ödevi
>
> >>>>>>>> Selamlar. Onun icin en pratik yol grafiklerin kesistigi noktayi
> >>>>>>>> bulmak. Ya da iki ayri fonksiyonun degerleri arasindaki farkin belli
> >>>>>>>> bir toleranstan dusuk oldugu bir degeri while dongusu icinde
> >>>>>>>> arayabilirsin. Ama analitik yol hiç pratik degil.  Kolay gelsin
> >>>>>>>> herkese
>

> >>>>>>>> On 5/3/09, muslum yaman <muslumya...@gmail.com> wrote:
>
> >>>>>>>> > arkadaşlar
>
> >>>>>>>> > tan(beta22)-1/(16*alfak*beta22*(eta)kare*(1-eta)kare)=0
> >>>>>>>> > dekleminin analitik çözümü varmı?Bulan varmı?
>
> >>>>>>>> --
> >>>>>>>> Can Cengiz
>
> >>>>> --

> >>>>> Umut Murat Gök- Alıntıyı gizle -
>
> - Alıntıyı göster -

[cihan orhan]

bana bakın. Herkes denizin kodu kullanacak mı? ben de kullanayım mı diye düşünüyorum da?

Saygılarımla,

Cihan ORHAN

04 Mayıs 2009 Pazartesi 23:23 tarihinde Müslüm Yaman <muslu...@gmail.com> yazdı: