İlk yörüngenin yarıçapı. Bohr'un önermeleri ve çizgi spektrumunun kökenine ilişkin açıklama

). Kemometride kullanılan genel amaçlı araçlar arasında özel yer MatLab paketini kaplar. Popülerliği alışılmadık derecede yüksek. Bunun nedeni MatLab'ın çok boyutlu verileri işlemek için güçlü ve çok yönlü olmasıdır. Paketin yapısı, onu matris hesaplamaları yapmak için uygun bir araç haline getirir. MatLab kullanılarak incelenebilecek problemlerin kapsamı şunları içerir: matris analizi, sinyal ve görüntü işleme, sinir ağları ve diğerleri. MatLab bir dildir yüksek seviye sahip olmak açık kaynak deneyimli kullanıcıların programlanan algoritmaları anlamalarını sağlar. Basit bir yerleşik programlama dili, kendi algoritmalarınızı oluşturmanızı kolaylaştırır. Uzun yıllardır MatLab'ı kullanarak şunu yarattık: büyük miktar işlevler ve ToolBox (özel araç paketleri). En popüler olanı Eigenvector Research, Inc.'in PLS ToolBox paketidir.

1. Temel bilgiler

1.1. MatLab çalışma ortamı

Programı başlatmak için simgeye çift tıklayın.

Şekilde gösterilen çalışma ortamı karşınıza açılacaktır. Çalışma ortamı MatLab 6.x

Şekilde gösterilen çalışma ortamı karşınıza açılacaktır. Çalışma ortamıÖnceki sürümlerin çalışma alanından biraz farklı olarak birçok destekleyici öğeye erişim için daha kullanışlı bir arayüze sahiptir.

    aşağıdaki unsurları içerir:

    düğmeler ve açılır liste içeren araç çubuğu; Başlatma Paneli sekmelerinin bulunduğu pencere veÇalışma alanı

    çeşitli ToolBox modüllerine ve çalışma tezgahı içeriklerine erişebileceğiniz; sekmeli pencere Komut Geçmişi Ve Geçerli Dizin

    önceden girilen komutları görüntülemek ve yeniden çağırmak ve ayrıca geçerli dizini ayarlamak için tasarlanmıştır;

    "giriş" istemini ve yanıp sönen dikey imleci içeren bir komut penceresi;

durum çubuğu. Çalışma ortamı Eğer bir çalışma ortamındaysanız Şekilde gösterilen pencerelerden bazıları eksikse Görünüm menüsünden uygun öğeleri seçmelisiniz: Komut Penceresi

Komutlar komut penceresine yazılmalıdır. Komut satırı istemini belirten » sembolünün yazılmasına gerek yoktur. Çalışma alanını görüntülemek için kaydırma çubuklarını veya sola veya sağa hareket etmek için Home, End tuşlarını ve yukarı veya aşağı hareket etmek için PageUp, PageDown tuşlarını kullanmak kullanışlıdır. Komut penceresinin çalışma alanında dolaştıktan sonra aniden yanıp sönen imlecin bulunduğu komut satırı kaybolursa, Enter tuşuna basmanız yeterlidir.

MatLab'ın o komutu çalıştırabilmesi veya ifadeyi değerlendirebilmesi için herhangi bir komut veya ifade yazmanın Enter tuşuna basılmasıyla bitmesi gerektiğini unutmamak önemlidir.

1.2. Basit hesaplamalar

Komut satırına 1+2 yazıp Enter tuşuna basın. Sonuç olarak MatLab komut penceresi aşağıdakileri görüntüler:

Pirinç. 2 Temel bileşen analizinin grafiksel gösterimi

MatLab programı ne yaptı? Önce 1+2 toplamını hesapladı, sonra sonucu özel bir değişken olan ans'a yazdı ve 3'e eşit olan değerini komut penceresinde gösterdi. Yanıtın altında, MatLab'ın daha sonraki hesaplamalar için hazır olduğunu gösteren, yanıp sönen bir imlecin bulunduğu bir komut satırı bulunur. Komut satırına yeni ifadeler yazabilir ve anlamlarını bulabilirsiniz.Önceki ifadeyle çalışmaya devam etmeniz gerekiyorsa, örneğin hesapla (1+2)/4.5), o zaman en kolay yol, ans değişkeninde saklanan mevcut sonucu kullanmaktır. ans/4.5 yazın (girerken ondalık sayılar

nokta kullanılır) ve tuşuna basın

Girmek

, çıkıyor

Pirinç. 3 Temel bileşen analizinin grafiksel gösterimi

1.3. Yankı komutları

Tüm değişken değerlerini kaydetmenin en kolay yolu Dosya menüsündeki Çalışma Alanını Farklı Kaydet seçeneğini kullanmaktır.

Bu, dizini ve dosya adını belirtmeniz gereken Çalışma Alanı Değişkenlerini Kaydet iletişim kutusunu açacaktır. Varsayılan olarak, dosyanın ana MatLab dizininin çalışma alt dizinine kaydedilmesi önerilir. Program, çalışmasının sonuçlarını mat uzantılı bir dosyaya kaydedecektir. Artık MatLab'ı kapatabilirsiniz. Bir sonraki çalışma oturumunda değişkenlerin değerlerini geri yüklemek için, kaydedilen bu dosyayı Dosya menüsünün Aç alt öğesini kullanarak açmalısınız. Artık son oturumda tanımlanan tüm değişkenler tekrar kullanılabilir. Yeni girilen komutlarda kullanılabilirler.

1.5. Dergi MatLab, yürütülebilir komutları ve sonuçları bir metin dosyasına yazma (bir çalışma günlüğü tutma) yeteneğine sahiptir; bu dosya daha sonra bir metin düzenleyiciden okunabilir veya yazdırılabilir. Günlüğe kaydetmeyi başlatmak için şu komutu kullanın: günlük MatLab, yürütülebilir komutları ve sonuçları bir metin dosyasına yazma (bir çalışma günlüğü tutma) yeteneğine sahiptir; bu dosya daha sonra bir metin düzenleyiciden okunabilir veya yazdırılabilir. Günlüğe kaydetmeyi başlatmak için şu komutu kullanın:.

Komut argümanı olarak

    iş günlüğünün saklanacağı dosyanın adını belirtmelisiniz.

    Yazılan diğer komutlar ve bunların yürütülmesinin sonuçları, örneğin bir dizi komut gibi bu dosyaya yazılacaktır.

    aşağıdaki eylemleri gerçekleştirir:

    examplepl-1.txt dosyasındaki günlüğü açar;

hesaplamalar yapar;

tüm değişkenleri MAT dosyasındaki work-1.mat dosyasına kaydeder;
günlüğü, MatLab kök dizininin çalışma alt dizinindeki example-1.txt dosyasına kaydeder ve MatLab'ı kapatır;
Bazı metin düzenleyicilerde examplel-1.txt dosyasının içeriğine bakın. Dosya aşağıdaki metni içerecektir:

a1=3;
a2=2,5;

a3=a1+a2

Çalışmayı kaydet-1 çıkış yapmak 1.6. Yardım sistemi Yardım menüsünden Yardım Penceresi seçeneğini seçtikten sonra veya araç çubuğundaki soru düğmesine tıkladığınızda MatLab Yardım penceresi görünür. Aynı işlem komutu yazılarak da yapılabilir. yardım kazan.

Tek tek konulara ilişkin yardım pencerelerini görüntülemek için şunu yazın:

yardım kazanma konusu

. Yardım penceresi size yardım komutuyla aynı bilgileri sağlar, ancak pencere arayüzü diğer yardım konularına daha kullanışlı bir bağlantı sağlar. Math Works Web sayfasının adresini kullanarak şirketin sunucusuna erişebilir ve en iyi şekilde yararlanabilirsiniz.

MatLab'ın büyük harfler ve harfler arasında ayrım yaptığını unutmayın. büyük harfler, dolayısıyla p ve P farklı değişkenlerdir. Dizilere (vektörler veya matrisler) girmek için bunların elemanları köşeli parantez içine alınır.

Dolayısıyla, 1x3'lük bir satır vektörü girmek için, satır öğelerinin boşluklarla veya virgüllerle ayrıldığı aşağıdaki komutu kullanın.

Bir sütun vektörüne girerken öğeler noktalı virgülle ayrılır.

Örneğin,

Küçük matrisleri doğrudan komut satırından girmek uygundur. Girişte bir matris, her elemanı bir satır vektörü olan bir sütun vektörü olarak düşünülebilir.

veya bir matris, her elemanı bir sütun vektörü olan bir satır vektörü olarak ele alınabilir.

2.2. Öğelere erişim

Matrislerin elemanlarına erişim iki indeks kullanılarak gerçekleştirilir - parantez içindeki satır ve sütun numaraları, örneğin, B(2,3) komutu B matrisinin ikinci satırının ve üçüncü sütununun elemanını döndürür. Bir matristen bir sütun veya satır seçmek için, matrisin sütun veya satır numarasını dizinlerden biri olarak kullanın ve diğer dizini iki nokta üst üste ile değiştirin. Örneğin A matrisinin ikinci satırını z vektörüne yazalım. Ayrıca iki nokta üst üste kullanarak matris bloklarını da seçebilirsiniz. Örneğin P matrisinden renkle işaretlenmiş bir blok seçelim. .

Çalışma ortamı değişkenlerini görüntülemeniz gerekiyorsa, komutu komut satırına yazmanız gerekir.

kim

Çalışma ortamının bir skaler (p), dört matris (A, B, P, P1) ve bir satır vektörü (z) içerdiği görülmektedir.

2.3. Temel Matris İşlemleri

Matris işlemlerini kullanırken, toplama veya çıkarma yaparken matrislerin aynı boyutta olması gerektiğini ve çarpma sırasında ilk matrisin sütun sayısının ikinci matrisin satır sayısına eşit olması gerektiğini unutmayın. Matrislerin, sayıların ve vektörlerin toplanması ve çıkarılması artı ve eksi işaretleri kullanılarak gerçekleştirilir ve çarpma yıldız işareti * ile işaretlenmiştir.

3×2 boyutunda bir matris tanıtalım

Bir matrisin bir sayıyla çarpılması da yıldız işareti kullanılarak yapılır ve hem sağda hem de solda bir sayıyla çarpabilirsiniz.

Yapı kare matris

^ operatörü kullanılarak bir tam sayının kuvveti üretilir P matrisini kendisiyle çarparak sonucunuzu kontrol edin.

2.4. Özel türde matrisler oluşturma Dikdörtgen bir matrisin sıfırlarla doldurulması yerleşik işlev tarafından gerçekleştirilir

MatLab matrisleri doldurma yeteneği sağlar rastgele sayılar. Fonksiyonun sonucu rand sıfır ile bir arasında eşit olarak dağıtılan sayıların bir matrisidir ve işlevler Randn- dağıtılan sayıların matrisi normal hukuk sıfır ortalama ve birim varyans ile.

İşlev teşhis elemanları diyagonal boyunca düzenleyerek bir vektörden diyagonal bir matris oluşturur.

2.5. Matris hesaplamaları

MatLab birçok şey içeriyor çeşitli işlevler matrislerle çalışmak için. Yani, örneğin, bir matrisin yerini değiştirmek kesme işareti kullanılarak yapılır "

Bulma ters matris fonksiyonu kullanılarak gerçekleştirilir yatırım kare matrisler için

3. MatLab ve Excel'in Bütünleştirilmesi

MatLab ve Excel'in entegre edilmesi, Excel kullanıcısının veri işleme, çeşitli hesaplamalar ve sonuçların görselleştirilmesi için çok sayıda MatLab işlevine erişmesine olanak tanır. excllink.xla eklentisi bu Excel uzantısını uygular. MatLab ve Excel arasındaki iletişim için özel fonksiyonlar tanımlanmıştır.

3.1. Excel Yapılandırması

Excel'i ayarlamadan önce birlikte çalışmak MatLab ile, MatLab'ın kurulu sürümünde Excel Link'in bulunduğundan emin olmalısınız. Ana MatLab dizininin exclink alt dizininde veya toolbox alt dizininde excllink.xla eklentisini içeren bir dosya bulunmalıdır. Excel'i başlatın ve Araçlar menüsünden Eklentiler'i seçin. Mevcut olanlarla ilgili bilgileri içeren bir iletişim kutusu açılacaktır. şu andaüst yapılar. Gözat düğmesini kullanarak excllink.xla dosyasının yolunu belirtin. Satır, iletişim kutusundaki eklentiler listesinde görünür

MatLab ile kullanım için Excel Bağlantısı 2.0

bayrak seti ile.

3.2. MatLab ve Excel arasında veri alışverişi

Excel'i başlatın, gerekli tüm ayarların önceki bölümde anlatıldığı gibi yapıldığını kontrol edin (MatLab kapatılmalıdır). Excel'in gerektirdiği şekilde ondalık basamakları ayırmak için nokta kullanarak A1'den C3'e kadar olan hücrelere bir matris girin.

Sayfadaki hücre verilerini seçin ve putmatrix düğmesine tıklayın, MatLab'ın çalışmadığını belirten bir uyarı içeren bir Excel penceresi açılır. Tamam'a tıklayın, MatLab'ın açılmasını bekleyin.

Seçilen Excel hücrelerindeki verilerin dışa aktarılması gereken MatLab tezgah değişkeninin adını belirtmek için bir giriş satırı içeren bir Excel iletişim kutusu görüntülenir. Örneğin M girin ve Tamam düğmesini kullanarak pencereyi kapatın.

MatLab komut penceresine gidin ve çalışma tezgahında üçe üçlük bir dizi içeren M değişkeninin oluşturulduğundan emin olun:

MatLab'da M matrisiyle bazı işlemler yapın, örneğin onu ters çevirin. yatırım Arama

Diğer MatLab komutları gibi bir matrisi ters çevirmek için bunu doğrudan Excel'den yapabilirsiniz. Excel Bağlantı panelinde bulunan evalstring düğmesine tıklamak, giriş satırında MatLab komutunu yazmanız gereken bir iletişim kutusunun görünmesine neden olur

IM=inv(M) .

Sonuç, MatLab ortamında komut çalıştırıldığında elde edilene benzer. Excel'e dönün, A5 hücresini geçerli hücre yapın ve getmatrix düğmesini tıklayın. Excel'e aktarılacak değişkenin adını girmenizi isteyen bir giriş satırı içeren bir iletişim kutusu görüntülenir. İÇİNDE

bu durumda

böyle bir değişken IM'dir. Tamam'ı tıklayın, A5'ten A7'ye kadar olan hücrelere ters matris öğeleri girilir. Bu nedenle, bir matrisi MatLab'a aktarmak için, Excel sayfasının uygun hücrelerini seçmelisiniz ve onu içe aktarmak için yalnızca içe aktarılan dizinin sol üst öğesi olacak bir hücre belirtmeniz gerekir. uygulamalar arasında bilgi alışverişi - kaynak veriler Excel'de bulunur, ardından MatLab'a aktarılır, orada bir şekilde işlenir ve sonuç Excel'e aktarılır. Kullanıcı, Excel Bağlantısı araç çubuğu düğmelerini kullanarak verileri aktarır.

Bilgi bir matris şeklinde sunulabilir, yani. çalışma sayfasının dikdörtgen alanı. Bir satır veya sütunda düzenlenen hücreler sırasıyla MatLab satır vektörlerine ve sütun vektörlerine aktarılır.

Satır vektörlerinin ve sütun vektörlerinin Excel'e aktarılması benzer şekilde gerçekleşir.

4. Programlama MatLab, yürütülebilir komutları ve sonuçları bir metin dosyasına yazma (bir çalışma günlüğü tutma) yeteneğine sahiptir; bu dosya daha sonra bir metin düzenleyiciden okunabilir veya yazdırılabilir. Günlüğe kaydetmeyi başlatmak için şu komutu kullanın: 4.1. M dosyaları Çok sayıda komut girmeniz ve bunları sık sık değiştirmeniz gerekiyorsa MatLab komut satırından çalışmak zorlaşır. Komut kullanarak günlük tutmak

ve çalışma ortamının sürdürülmesi çalışmayı biraz daha kolaylaştırır. En çok

uygun bir şekilde MatLab komut gruplarını yürütmek, komutları yazabileceğiniz, hepsini bir kerede veya parçalar halinde yürütebileceğiniz, bunları bir dosyaya kaydedebileceğiniz ve daha sonra kullanabileceğiniz M dosyalarının kullanılmasıdır. M dosyası düzenleyicisi M dosyalarıyla çalışacak şekilde tasarlanmıştır. Onun yardımıyla, kendi işlevlerinizi oluşturabilir ve bunları komut penceresinden de dahil olmak üzere çağırabilirsiniz. Ana MatLab penceresinin Dosya menüsünü genişletin ve Yeni öğede M dosyası alt öğesini seçin. Yeni dosya, şekilde gösterilen M dosyası düzenleyici penceresinde açılır. MatLab'da iki tür M dosyası vardır: program dosyası (

Komut Dosyası M Dosyaları

), bir dizi komut ve dosya işlevlerini içeren, (

İşlev M Dosyaları ), kullanıcı tanımlı işlevleri açıklar. 4.2. Dosya programı

Düzenleyiciye, bir grafik penceresinde iki grafiğin oluşturulmasına yol açan komutları girin

Şimdi düzenleyicinin Dosya menüsünden Farklı kaydet'i seçerek mydemo.m adlı dosyayı ana MatLab dizininin çalışma alt dizinine kaydedin. Dosyanın içerdiği tüm komutları çalıştırmak için Hata Ayıklama menüsünden Çalıştır'ı seçin. Ekranda bir grafik penceresi görünecektir ), kullanıcı tanımlı işlevleri açıklar.. Sol düğmeyi basılı tutarken fareyi kullanmayı veya tuşu basılı tutarken ok tuşlarını kullanmayı seçin Vardiya

, ilk dört komutu seçin ve bunları Metin öğesinden yürütün. Lütfen grafik penceresinde yürütülen komutlara karşılık gelen yalnızca bir grafiğin görüntülendiğini unutmayın.

Bazı komutları yürütmek için bunları seçip F9 tuşuna basmayı unutmayın.

M dosyasının bireysel blokları, yürütme sırasında atlanan ancak M dosyasıyla çalışırken kullanışlı olan yorumlarla sağlanabilir. Yorumlar yüzde işaretiyle başlar ve otomatik olarak yeşil renkte vurgulanır, örneğin:

Mevcut bir M dosyasının açılması, çalışma ortamının Dosya menüsündeki Aç öğesi veya M dosyası düzenleyicisi kullanılarak yapılır. 4.3. Dosya işlevi Yukarıda tartışılan dosya programı yalnızca bir MatLab komutları dizisidir; giriş veya çıkış argümanlarına sahip değildir. Kullanmak sayısal yöntemler ve MatLab'da kendi uygulamalarınızı programlarken, üreten dosya fonksiyonları oluşturabilmeniz gerekir.

gerekli eylemler

giriş argümanlarıyla ve eylemin sonucunu çıkış argümanlarında döndürün. Dosya işlevleriyle nasıl çalışılacağını anlamanıza yardımcı olacak birkaç basit örneğe bakalım.

Artık oluşturulan işlev yerleşik sin, cos ve diğerleriyle aynı şekilde kullanılabilir. Kendi işlevleri bir dosya programından ve başka bir dosya işlevinden çağrılabilir. Matrisleri ölçeklendirecek bir dosya işlevini kendiniz yazmaya çalışın; her sütunu o sütunun standart sapmasına bölün.

Virgülle ayrılmış bir listeye yerleştirilen çeşitli giriş bağımsız değişkenlerini içeren bir işlev dosyası yazabilirsiniz. Birden fazla değer döndüren işlevler de oluşturabilirsiniz. Bunu yapmak için, çıktı bağımsız değişkenleri, virgüllerle ayrılmış olarak çıktı bağımsız değişkenleri listesine eklenir ve listenin kendisi köşeli parantez içine alınır. İyi bir örnek saniye cinsinden belirtilen süreyi saat, dakika ve saniyeye dönüştüren bir fonksiyondur.

Çoklu çıktı argümanlarıyla dosya fonksiyonları çağrılırken sonuç uygun uzunlukta bir vektöre yazılmalıdır.

4.4 Grafik oluşturma

MatLab'ın aşağıdakiler için geniş fırsatları var: grafik görüntü vektörler ve matrislerin yanı sıra yorumlar oluşturmak ve grafikleri yazdırmak için. Birkaç önemli grafik fonksiyonunu tanımlayalım.

İşlev komplo sahip olmak çeşitli şekiller giriş parametreleriyle ilişkili; örneğin, arsa(y), y'nin öğelerinin indekslerine göre parçalı doğrusal bir grafiğini oluşturur. Eğer argüman olarak iki vektör verilirse, arsa(x,y) y'ye karşı x'in bir grafiğini üretecektir.

Örneğin, günah fonksiyonunu 0 ila 2π aralığında çizmek için aşağıdakileri yapın: ), kullanıcı tanımlı işlevleri açıklar.

Program, pencerede görüntülenen bir bağımlılık grafiği oluşturmuştur.

MatLab, her çizime otomatik olarak farklı bir renk atar (kullanıcı bunu yapmadığı sürece), böylece veri kümeleri arasında ayrım yapabilirsiniz. Takım devam etmek Mevcut bir grafiğe eğriler eklemenizi sağlar. İşlev alt olay örgüsü

birden fazla grafiği tek bir pencerede görüntülemenizi sağlar

4.5 Grafiklerin yazdırılması Dosya menüsündeki Yazdır öğesi ve komut baskı baskı MatLab grafikleri Dosya menüsündeki Yazdır öğesi ve komut. Yazdır menüsü, ortak standart yazdırma seçeneklerini seçmenizi sağlayan bir iletişim kutusu açar. Takım

çıktıda daha fazla esneklik sağlar ve M dosyalarından yazdırma üzerinde kontrole olanak tanır.

Sonuç doğrudan varsayılan yazıcıya gönderilebilir veya belirli bir dosyaya kaydedilebilir.

5.1. Merkezleme ve ölçeklendirme

Çoğu zaman analiz sırasında orijinal verileri dönüştürmek gerekir. En sık kullanılan veri dönüştürme yöntemleri, her değişkeni ortalamak ve ölçeklendirmektir. standart sapma. Matrisin merkezlenmesi için fonksiyon kodu verildi. Bu nedenle aşağıda sadece fonksiyonun kodu gösterilmektedir. terazi veri. Orijinal matrisin ortalanması gerektiğini lütfen unutmayın

fonksiyon Xs = ölçeklendirme(X)
% ölçeklendirme: çıkış matrisi X'tir
% matris X ortalanmalıdır

Xs = X * inv(diag(std(X)));

Ölçeklendirmenin sonu

5.2. SVD/PCA

Verileri sıkıştırmanın en popüler yolu çok değişkenli analiz temel bileşen analizidir (PCA). İLE matematiksel nokta PCA orijinal matrisin ayrıştırılmasıdır X, yani onu iki matrisin çarpımı olarak temsil etmek Komut Geçmişi T

X = P TP t+

e onu iki matrisin çarpımı olarak temsil etmek Matris

puanların (puanların) matrisi denir, matris bir artıklar matrisidir. onu iki matrisin çarpımı olarak temsil etmek Komut Geçmişi T Matrisleri bulmanın en basit yolu - standart aracılığıyla SVD ayrıştırmasını kullanın MatLab işlevi , isminde .

SVD

işlev = pcasvd(X)
Svd(X);
T=U*D;

P=V;

%pcasvd'nin sonu

5.3 PCA/NIPALLER X PCA hesapları ve yükleri oluşturmak için her adımda bir bileşeni hesaplayan yinelenen NIPALS algoritması kullanılır. İlk olarak orijinal matris t+ 0 , dönüştürülür (en azından merkezli; bkz.) ve bir matrise dönüşür A =0..

Sonraki başvuru 2. sonraki algoritma T Sonraki başvuru P t+t = / Sonraki başvuru P Sonraki başvuru 3. sonraki algoritma = sonraki algoritma / (sonraki algoritma P sonraki algoritma T Sonraki başvuru = t+t = sonraki algoritma / sonraki algoritma P sonraki algoritma A

) ½ 4. dönüştürülür (en azından merkezli; bkz.) ve bir matrise dönüşür 5. Yakınsamayı kontrol edin, değilse 2'ye gidin Sonraki başvurut ==Sonraki başvuru Sonrakini hesapladıktan sonra ( sonraki algoritmat ==sonraki algoritma t+ dönüştürülür (en azından merkezli; bkz.) ve bir matrise dönüşür+1 = t+t =Sonraki başvuru sonraki algoritma dönüştürülür (en azından merkezli; bkz.) ve bir matrise dönüşür-th) bileşenlerin olduğunu varsayıyoruz dönüştürülür (en azından merkezli; bkz.) ve bir matrise dönüşür+1.

Ve X.

Açık

NIPALS algoritmasının kodu bu kılavuzda okuyucuların kendileri tarafından yazılabilir; yazarlar kendi versiyonlarını sunarlar. PCA'yı hesaplarken, ana bileşenlerin sayısını (değişken sayıPC) girebilirsiniz.
Kaç bileşene ihtiyaç duyulduğunu bilmiyorsanız komut satırına = pcanipals (X) yazmalısınız ve ardından program, bileşen sayısını orijinal matrisin boyutlarından en küçüğüne eşit olarak ayarlayacaktır.

fonksiyon = pcanipals(X, sayıPC)
Bileşen sayısının yüzde olarak hesaplanması
= boyut(X); P=; T=;< X_c
Lenf(sayıPC) > 0 ise
pc = sayıPC(1);
elseif (uzunluk(sayıPC) == 0) & X_r
pc = X_r;

başka
pc = X_c;
son;

k = 1:pc için
P1 = rand(X_c, 1); T1 = X * P1; d0 = T1"*T1;
P1 = (T1" * X/(T1" * T1))"; P1 = P1/norm(P1); T1 = X * P1; d = T1" * T1;
d - d0 > 0,0001 iken;

P1 = (T1" * X/(T1" * T1)); P1 = P1/norm(P1); T1 = X * P1; d0 = T1"*T1;
d - d0 > 0,0001 iken;

Chemometrics eklentisini kullanarak PCA'nın nasıl hesaplanacağı eğitimde açıklanmaktadır.

5.4PLS1

Çok değişkenli kalibrasyon için en popüler yöntem Gizli Yapılara Projeksiyon (PLS) yöntemidir. Bu yöntem, öngörücü matrisin eşzamanlı ayrıştırılmasını içerir X ve yanıt matrisleri e:

X=P t+ t+ e=UQ t+ F onu iki matrisin çarpımı olarak temsil etmek=XW(T P W) –1

Projeksiyon, karşılık gelen vektörler arasındaki korelasyonu en üst düzeye çıkaracak şekilde tutarlı bir şekilde oluşturulmuştur. X-hesaplar Sonraki başvurut = Komut Geçmişi e-hesaplar sent =. e Veri bloğu ise birden fazla yanıt içerir (ör. k >1), başlangıç ​​verilerinin iki projeksiyonu oluşturulabilir – PLS1 ve PLS2. İlk durumda, yanıtların her biri için sen k onu iki matrisin çarpımı olarak temsil etmek (kendi projeksiyon alt uzayı inşa edilmiştir. Aynı zamanda faturalar sen T (W, ) ve yükler Q

) hangi yanıtın kullanıldığına bağlıdır. Bu yaklaşıma PLS1 denir. X Sonrakini hesapladıktan sonra ( e PLS2 yöntemi için tüm yanıtlar için ortak olan yalnızca bir projeksiyon alanı oluşturulur.

Bu kitapta PLS yönteminin ayrıntılı bir açıklaması verilmiştir. PLS1 hesapları ve yükleri oluşturmak için tekrarlayan bir algoritma kullanılır. İlk olarak orijinal matrisler
merkez

= mc(X); t+= mc(Y); ve bir matrise dönüşüyorlar 0 , dönüştürülür (en azından merkezli; bkz.) ve bir matrise dönüşür 0 ve vektör

1. F T ve bir matrise dönüşüyorlart ==0. Daha sonra bunlara aşağıdaki algoritma uygulanır. t+ dönüştürülür (en azından merkezli; bkz.) ve bir matrise dönüşür 2. F = F / (F P F w Sonraki başvuru = t+t = F 4. T = Sonraki başvuru P ve bir matrise dönüşüyorlart = / Sonraki başvuru P Sonraki başvuru 5. sen = Tve bir matrise dönüşüyorlart = / T 2 6. sonraki algoritma T Sonraki başvuru P t+t = / Sonraki başvuru P Sonraki başvuru

) ½ 4. dönüştürülür (en azından merkezli; bkz.) ve bir matrise dönüşür 5. Yakınsamayı kontrol edin, değilse 2'ye gidin Sonraki başvurut ==Sonraki başvuru Sonrakini hesapladıktan sonra ( sonraki algoritmat ==sonraki algoritma) ½ 3. t+ dönüştürülür (en azından merkezli; bkz.) ve bir matrise dönüşür+1 = t+t =Sonraki başvuru sonraki algoritma Q dönüştürülür (en azından merkezli; bkz.) ve bir matrise dönüşür-th) bileşenlerin olduğunu varsayıyoruz dönüştürülür (en azından merkezli; bkz.) ve bir matrise dönüşür+1.

. Bir sonraki bileşeni elde etmek için kalanları hesaplamanız gerekir.

t ve indeksi değiştirerek aynı algoritmayı onlara uygulayın
İşte bu algoritmanın kitaptan alınan kodu
fonksiyon = pls(x, y)
%
%PLS: PLS bileşenini hesaplar.

%Çıktı vektörleri w, t, u, q ve p'dir.

% y'den u başlangıç ​​vektörü olarak bir vektör seçin.
u = y(:, 1);

% Yakınsama kriteri çok yüksek ayarlanmış.
kri = 100;

% Buradan sona kadar olan komutlar yakınsayana kadar tekrarlanır.
iken (kri > 1e - 10)
% Her başlangıç ​​vektörü ueski olarak kaydedilir.
ueski = sen; w = (u" * x)"; w = w/norm(w);

t = x * w; q = (t" * y)"/(t" * t);
u = y * q/(q" * q);
pc = X_r;

% Yakınsama kriteri u-uold normunun u normuna bölünmesidir.
kri = norm(ueski - u)/norm(u);

% Yakınsama sonrasında p'yi hesaplayın.

p = (t" * x)"/(t" * t); % Pls'in sonuEklentiyi kullanarak PLS1'i hesaplama hakkında Kemometri

Ekle

Excel'deki Projeksiyon yöntemleri kılavuzunda açıklanmıştır. X Sonrakini hesapladıktan sonra ( e 5.5PLS2 t+ PLS2 için algoritma aşağıdaki gibidir. İlk olarak orijinal matrisler F 0 , dönüştürülür (en azından merkezli; bkz.) ve bir matrise dönüşür dönüşürler (en azından - merkez; bkz.) ve matrislere dönüşürler

0 ve sen 2. F T sen P t+ dönüştürülür (en azından merkezli; bkz.) ve bir matrise dönüşür 3. F = F / (F P F T Sonraki başvuru = t+t = F 5. T T Sonraki başvuru P Ft = / Sonraki başvuru P Sonraki başvuru 6. sen = Ft = T/ T P T=0. Daha sonra bunlara aşağıdaki algoritma uygulanır. sonraki algoritma T Sonraki başvuru P t+t = / Sonraki başvuru P Sonraki başvuru

) ½ 4. dönüştürülür (en azından merkezli; bkz.) ve bir matrise dönüşür 1. Başlangıç ​​vektörünü seçin Sonraki başvurut ==Sonraki başvuru, sonraki algoritmat ==7. Yakınsamayı kontrol edin, değilse 2 8'e gidin.t ==F, sent ==sen Sonrakini hesapladıktan sonra ( T oh) PLS2 bileşenlerinin yerleştirilmesi gerekiyor: T) ½ 3. t+ dönüştürülür (en azından merkezli; bkz.) ve bir matrise dönüşür+1 = t+t =p, w bir = Ft = +1 = F dönüştürülür (en azından merkezli; bkz.) ve bir matrise dönüşürt p t ve dönüştürülür (en azından merkezli; bkz.) ve bir matrise dönüşür-th) bileşenlerin olduğunu varsayıyoruz dönüştürülür (en azından merkezli; bkz.) ve bir matrise dönüşür+1.

İşte yine kitaptan ödünç alınan kod.

fonksiyon = plsr(x, y, a)
% PLS: PLS bileşenini hesaplar.
% Çıkış matrisleri W, T, U, Q ve P'dir.
% B regresyon katsayılarını ve SS toplamlarını içerir
% kareler için artıklar.
% a, bileşenlerin sayısıdır.
%
% Bir bileşen için: bitirmek için tüm komutları kullanın.

i=1:a için
% Karelerin toplamını hesaplayın. ss fonksiyonunu kullanın.
sx = ;
sistem = ;

% Bir bileşeni hesaplamak için pls fonksiyonunu kullanın.
= pls(x, y);

% Artıkları hesaplayın.
x = x - t * p";
y = y - t * q";

% Vektörleri matrislere kaydedin.
W = ;
T = ;
u =;
S =;
p =;
pc = X_r;

% Döngüden sonraki regresyon katsayılarını hesaplayın.
B=W*inv(P"*W)*Q";

% Nihai artık SS'yi kareler vektörlerinin toplamına ekleyin.
sx=;
sy=;

% X ve Y için ss vektörlerinden bir matris oluşturun.
SS = ;

%Kullanılan SS oranını hesaplayın.
= boyut(SS);
tt = (SS * diag(SS(1,:).^(-1)) - birler(a, b)) * (-1)

% Plsr sonu

fonksiyon = ss(x)
%SS: X matrisinin karelerinin toplamını hesaplar.
%

ss=toplam(toplam(x. * x));
%Ss sonu

Eklentiyi kullanarak PLS2'yi hesaplama hakkında % Pls'in sonuEklentiyi kullanarak PLS1'i hesaplama hakkında Kemometri

Çözüm

MatLab veri analizi için çok popüler bir araçtır. Ankete göre, tüm araştırmacıların üçte biri bu programı kullanırken, Unsrambler programı bilim adamlarının yalnızca %16'sı tarafından kullanılıyor. MatLab'ın en büyük dezavantajı yüksek fiyat

. Ayrıca MatLab rutin hesaplamalar için de iyidir. Etkileşim eksikliği, yeni, keşfedilmemiş veri kümeleri için arama ve araştırma hesaplamaları yaparken bunu sakıncalı hale getirir.

13.7. Bohr'un teorisine göre hidrojen atomu (hidrojen benzeri atom)

13.7.3. Bir atomdaki elektron yörüngeleri Buna göre ( yörünge nicemleme kuralı Somerfeld ilkesi

) bir atomdaki elektronun durağan durumlarının enerjisi, yörüngesinin yarıçapı ve bu yörüngedeki hız arasındaki ilişki formülle verilir.

mvr = nℏ,

burada m elektron kütlesidir, m = 9,11 ⋅ 10 −31 kg; v - elektron hızı; r, elektron yörüngesinin yarıçapıdır; ℏ - azaltılmış Planck sabiti, ℏ = h /2π ≈ 1,055 ⋅ 10 −34 J ⋅ s; h, Planck sabitidir, h = 6,626 ⋅ 10 −34 J ⋅ s; n baş kuantum sayısıdır.

Yörüngelerin kuantizasyonu kuralından, bir atomdaki bir elektronun durağan durumlarının yalnızca koşulun karşılandığı elektron yörüngelerine karşılık geldiği sonucu çıkar.

burada r n, n numaralı yörüngedeki elektronun yarıçapıdır; v n - n numaralı yörüngedeki elektron hızı; m elektron kütlesidir, m = 9,11 ⋅ 10 −31 kg; ℏ - azaltılmış Planck sabiti, ℏ = h /2π ≈ 1,055 ⋅ 10 −34 J ⋅ s; h, Planck sabitidir, h = 6,626 ⋅ 10 −34 J ⋅ s; n baş kuantum sayısıdır.

Sabit elektron yörüngesinin yarıçapı

r n = ℏ 2 n 2 k Z e 2 m ,

burada k = 1/4πε 0 ≈ 9 ⋅ 10 9 N ⋅ m2 /Cl 2; ε 0 - elektrik sabiti, ε 0 = 8,85 ⋅ 10 −12 F/m; Z - elemanın seri numarası; e elektron yüküdür, e = −1,6 ⋅ 10 −19 C; m elektron kütlesidir, m = 9,11 ⋅ 10 −31 kg; ℏ - azaltılmış Planck sabiti, ℏ = h /2π ≈ 1,055 ⋅ 10 −34 J ⋅ s; h, Planck sabitidir, h = 6,626 ⋅ 10 −34 J ⋅ s; n baş kuantum sayısıdır.

İlk yörünge yarıçapı hidrojen atomundaki elektron (Z = 1 ve n = 1) eşittir

r 1 = ℏ 2 k e 2 m = 0,53 ⋅ 10 − 10 m

ve denir ilk Bohr yarıçapı.

Hesaplamaları basitleştirmek için yarıçap n'inci yörüngeler Hidrojen benzeri bir atomdaki elektron için formülü kullanın

r (Å) = 0,53 ⋅ n 2 Z,

burada r (Å) angstrom cinsinden yarıçaptır (1 Å = 1,0 ⋅ 10 −10 m); Z - seri numarası kimyasal element V Periyodik tablo elementler D.I. Mendeleev; n = 1, 2, 3, … ana kuantum sayısıdır.

Sabit bir yörüngede elektronun hızı hidrojen benzeri bir atomda formülle belirlenir

v n = k Z e 2 n ℏ ,

burada k = 1/4πε 0 ≈ 9 ⋅ 10 9 N ⋅ m2 /Cl 2; ε 0 - elektriksel sabit, ε 0 = 8,85 ⋅ 10 −12 F/m; Z - elemanın seri numarası; e elektron yüküdür, e = −1,6 ⋅ 10 −19 C; ℏ - azaltılmış Planck sabiti, ℏ = = h /2π ≈ 1,055 ⋅ 10 −34 J ⋅ s; h, Planck sabitidir, h = 6,626 ⋅ 10 −34 J ⋅ s; n baş kuantum sayısıdır.

Birinci yörüngedeki elektron hızı bir hidrojen atomunda (Z = 1 ve n = 1) eşittir

v n = k e 2 ℏ = 2,2 ⋅ 10 6 m/s.

Değer hesaplamalarını basitleştirmek için elektron hızı n'inci yörünge hidrojen benzeri bir atomda formül kullanılır

v (m/s) = 2,2 ⋅ 10 6 ⋅ Zn ,

burada v (m/s) - m/s cinsinden hız modülü; Z, Periyodik Elementler Tablosu D.I'deki bir kimyasal elementin seri numarasıdır. Mendeleyev; n = 1, 2, 3, … ana kuantum sayısıdır.

Örnek 21. Helyum atomundaki bir elektron, birinci yörüngeden yarıçapı 9 kat daha büyük olan bir yörüngeye doğru hareket eder. Atomun absorbe ettiği enerjiyi bulunuz.

Çözüm . Bir helyum atomu tarafından emilen enerji, enerji farkına eşittir:

∆E = E 2 - E 1 ,

burada E1, r1 yörünge yarıçapına karşılık gelen elektron enerjisidir; E2, r2 yörünge yarıçapına karşılık gelen elektron enerjisidir.

Bir helyum atomundaki (Z = 2) elektron enerjileri aşağıdaki formüllerle belirlenir:

  • baş kuantum sayısı n 1 = 1 olan bir durumda -

E 1 (eV) = - 13,6 Z 2 n 1 2 = - 54,4 eV;

  • baş kuantum sayısı n 2 olan durum -

E 2 (eV) = − 54,4 n 2 2 .

E2 enerjisini belirlemek için karşılık gelen yörüngelerin yarıçapları için ifadeyi kullanırız:

  • temel kuantum sayısı n 1 = 1 olan bir yörünge için -

r 1 (Å) ≈ 0,53 n 1 2 Z = 0,265 Å;

  • baş kuantum sayısı n 2 olan yörüngeler -

r 2 (Å) ≈ 0,265 n 2 2.

Yarıçap oranı

r 2 (Å) r 1 (Å) = 0,265 n 2 2 0,265 = n 2 2

ikinci durumun baş kuantum sayısını belirlememizi sağlar:

n 2 = r 2 (Å) r 1 (Å) = 9 = 3,

burada r2/r1, koşulda belirtilen yörünge yarıçaplarının oranıdır, r2/r1 = 9.

Enerji oranından

E 2 E 1 = 1 n 2 2

bundan ikinci durumdaki bir helyum atomundaki bir elektronun enerjisinin olduğu sonucu çıkar:

E 2 = E 1 n 2 2 = − 54,4 eV 3 2 = − 6,04 eV.

Belirtilen geçiş sırasında atom tarafından emilen enerji farktır

∆E = E 2 − E 1 = −6,04 − (−54,4) = 48,4 eV.

Sonuç olarak, belirtilen geçiş sırasında atom 48,4 eV'ye eşit bir enerji emdi.

Örnek 1. Hidrojen atomu için ilk Bohr yörüngesinin yarıçapını ve üzerindeki elektronun hızını hesaplayın.

Çözüm. Yarıçap n'inci Bohr yörüngesi r n ve hız senüzerindeki elektronlar Bohr'un ilk önermesinin denklemiyle birbirine bağlanır:

ben n r n = n. (3.1)

Miktarlarla ilgili başka bir denklemin olması sen Ve r nÇekirdeğin Coulomb çekim kuvvetinin etkisi altında dairesel bir yörüngede hareket eden bir elektron için Newton'un ikinci yasasını yazıyoruz. Hidrojen atomunun çekirdeğinin, yükü mutlak değer olarak bir elektronun yüküne eşit olan bir proton olduğunu göz önünde bulundurarak şunu yazarız:

Nerede M– elektron kütlesi, – normal ivme. (3.1) ve (3.2)'yi birlikte çözersek şunu elde ederiz:

Onu buraya koyuyorum n=1, hesaplamaları yapalım:

; .

Örnek 2. Hidrojen atomundaki elektron dördüncü enerji seviyesinden ikinci enerji seviyesine geçmiştir. Yayılan fotonun enerjisini ve dalga boyunu belirleyin.

Çözüm. Foton enerjisini belirlemek için hidrojen benzeri iyonlara ilişkin seri formülü kullanırız:

, (3.3)

Nerede λ – foton dalga boyu; RRydberg sabiti; Z– göreceli birimlerdeki nükleer yük (at Z= 1 formülü hidrojenin seri formülüne girer); n 1– elektronun hareket ettiği yörüngenin numarası; n 2– elektronun hareket ettiği yörünge numarası ( n 1 Komut Geçmişi n 2– temel kuantum sayıları).

Foton enerjisi E aşağıdaki formülle ifade edilir

Bu nedenle eşitliğin her iki tarafı (13.3) ile çarpılır. hc foton enerjisi için bir ifade elde ederiz:

.

Çünkü Rhc iyonlaşma enerjisidir E ben o zaman hidrojen atomu

.

Eşitlik (3.4)'ten fotonun dalga boyunu ifade ediyoruz

Hesaplamaları sistem dışı birimlerde yapacağız: E ben= 13,6 eV; Z = 1; n 1 = 2; n 2 = 4:

eV = 2,55 ev.

M.

Örnek 3. Başlangıç ​​hızı ihmal edilebilecek bir elektron, hızlanan bir potansiyel farkından geçmiştir. kendi projeksiyon alt uzayı inşa edilmiştir. Aynı zamanda faturalar. İki durum için elektronun de Broglie dalga boyunu bulun: 1) U 1= 51V; 2) U 2= 510kV.

Çözüm. Bir parçacığın de Broglie dalga boyu onun momentumuna bağlıdır R ve formülle belirlenir

Nerede H– Planck sabiti.

Bir parçacığın momentumu, kinetik enerjisi biliniyorsa belirlenebilir. T. Momentum ve kinetik enerji arasındaki ilişki, göreli olmayan durum için (parçacığın kinetik enerjisinin dinlenme enerjisinden çok daha az olduğu zaman) ve göreli durum için (kinetik enerjinin parçacığın dinlenme enerjisi ile karşılaştırılabilir olduğu zaman) farklıdır.

Göreli olmayan durumda

Nerede m 0 parçacığın geri kalan kütlesidir.

Göreli durumda

, (3.7)

Nerede E 0 = m 0 sn 2– parçacığın dinlenme enerjisi.

(3.6) ve (3.7) ilişkilerini dikkate alan formül (3.5) yazılacaktır:

Göreli olmayan durumda

Göreli durumda

. (3.9)

Problem koşullarında belirtilen potansiyel farkından geçen bir elektronun kinetik enerjilerini karşılaştıralım. U 1= 51V ve U 2= 510 kV, elektronun geri kalan enerjisi ile ve buna bağlı olarak de Broglie dalga boyunu hesaplamak için (3.8) veya (3.9) formüllerinden hangisinin kullanılması gerektiğine karar vereceğiz.


Bilindiği gibi hızlanan bir potansiyel farkından geçen elektronun kinetik enerjisi kendi projeksiyon alt uzayı inşa edilmiştir. Aynı zamanda faturalar,

onu iki matrisin çarpımı olarak temsil etmek = AB.

İlk durumda T 1 = AB 1= 51 eV = 0,51 · 10 -4 MeV, elektronun dinlenme enerjisinden çok daha azdır E 0 = m 0 sn 2= 0,51 MeV. Bu nedenle bu durumda formül (3.8)'i uygulayabiliriz. Hesaplamaları basitleştirmek için şunu unutmayın: T 1 = 10 -4 m 0 c 2. Bu ifadeyi formül (3.8)'de yerine koyarak, formda yeniden yazıyoruz.

.

Compton dalga boyu olduğu düşünülürse λ , alıyoruz

Çünkü λ = 14.43, o halde

İkinci durumda kinetik enerji T2 = AB 2= 510 keV = 0,51 MeV, yani. elektronun dinlenme enerjisine eşittir. Bu durumda başvurulması gerekmektedir. göreli formül(3.9). Bunu göz önünde bulundurarak T2= 0,51 MeV = m 0 sn 2, formülü (3.9) kullanarak buluyoruz

,

λ değerini yerine koyalım ve hesaplamaları yapalım:

Örnek 4. Hidrojen atomundaki bir elektronun kinetik enerjisi büyüklük sırasına göredir. T= 10 eV. Belirsizlik ilişkisini kullanarak bir atomun minimum doğrusal boyutlarını tahmin edin.

Çözüm. Konum ve momentum için belirsizlik ilişkisi şu şekildedir:

Nerede Dx– parçacığın (bu durumda elektron) koordinatının belirsizliği; Dr x– parçacık (elektron) momentumunun belirsizliği; – Planck sabiti.

Belirsizlik ilişkisinden, bir parçacığın uzaydaki konumu ne kadar doğru belirlenirse, momentumun ve dolayısıyla parçacığın enerjisinin de o kadar belirsiz hale geldiği sonucu çıkar. Atomun doğrusal boyutları olsun ben o zaman atomun elektronu belirsizliğin olduğu bölgede bir yere yerleştirilecektir.

Ayrıklığın varlığı enerji seviyeleri atomların (aynı zamanda moleküllerin ve atom çekirdeğinin) temel bir özelliğidir.

Enerji seviyelerinin farklılığını açıklayan atomun yapısını hayal etmek için bildiğimiz fizik yasalarını uygulamaya çalışalım.

En basit atom olan hidrojen atomunu ele alalım. Periyodik element tablosundaki hidrojenin sıra sayısı bire eşit Bu nedenle bir hidrojen atomu, yükü eşit olan pozitif bir çekirdek ve bir elektrondan oluşur. Çekirdek ile elektron arasında yükler arasında bir çekim kuvveti vardır. Bu kuvvetin varlığı radyal (merkezcil) ivme sağlar; bu sayede hafif bir elektron, ağır bir çekirdeğin etrafında dairesel veya eliptik bir yörüngede döner, tıpkı bir gezegenin yerçekiminin etkisi altında Güneş'in etrafında dönmesi gibi. Dolayısıyla atomun farklı olası durumları, çekirdek etrafında dönen elektronun yörüngesinin boyutundaki (ve şeklindeki) farklılığa karşılık gelir.

Bir atomdaki elektronun enerjisi, yörüngesindeki hareketin kinetik enerjisinden ve çekirdeğin elektrik alanındaki potansiyel enerjisinden oluşur. Dairesel bir yörüngedeki bir elektronun enerjisinin ve dolayısıyla bir bütün olarak atomun enerjisinin yörüngenin yarıçapına bağlı olduğu gösterilebilir (paragrafın sonuna bakın): yörüngenin daha küçük bir yarıçapı buna karşılık gelir. atomun daha küçük bir enerjisine. Ancak § 204'te gördüğümüz gibi, bir atomun enerjisi herhangi bir değeri değil, yalnızca belirli seçilmiş değerleri alabilir. Enerji yörüngenin yarıçapı tarafından belirlendiğinden, atomun her enerji seviyesi seçilen belirli bir yarıçaptaki bir yörüngeye karşılık gelir.

Bir hidrojen atomundaki bir elektronun olası dairesel yörüngelerinin resmi Şekil 2'de gösterilmektedir. 367. Bir atomun ana enerji seviyesi, en küçük yarıçaptaki bir yörüngeye karşılık gelir.

Pirinç. 367. Hidrojen atomundaki bir elektronun olası yörüngeleri: yörüngelerin yarıçapı orantılı olarak artar, yani. ilişki içinde vb.

Normalde elektron bu yörüngededir. Yeterli miktarda enerji verildiğinde elektron başka bir enerji düzeyine geçer, yani dış yörüngelerden birine “sıçrar”. Belirtildiği gibi böyle uyarılmış bir durumda atom kararsızdır. Bir süre sonra elektron daha düşük bir seviyeye hareket eder, yani daha küçük yarıçaplı bir yörüngeye “sıçrar”. Bir elektronun uzak bir yörüngeden yakın bir yörüngeye geçişine bir ışık kuantumunun yayılması eşlik eder.

Dolayısıyla, atomun nükleer modelinden ve enerji seviyelerinin farklılığından, atomda seçilmiş, "izin verilen" elektron yörüngelerinin varlığı ortaya çıkar. Bir elektronun neden rastgele bir yarıçaptaki bir yörüngede bir çekirdeğin etrafında dönemediği sorusu ortaya çıkıyor. ne içinde fiziksel fark izin verilen ve yasa dışı yörüngeler?

Ders kitabının önceki bölümlerinden aşina olduğumuz mekanik ve elektrik yasaları (bkz. cilt I, II) bu sorulara herhangi bir yanıt vermiyor. Bu yasalar açısından bakıldığında tüm yörüngeler tamamen eşittir. Özel yörüngelerin varlığı bu yasalarla çelişmektedir.

Bildiğimiz fizik yasalarıyla eşit derecede çarpıcı bir çelişki, atomun (temel durumdaki) kararlılığıdır. İvmeyle hareket eden herhangi bir yükün elektromanyetik dalgalar yaydığını biliyoruz. Elektromanyetik radyasyon enerjinizi yanınızda götürün. Bir atomda elektron, küçük yarıçaplı bir yörüngede yüksek hızda hareket eder ve bu nedenle büyük bir merkezcil ivmeye sahiptir. Bildiğimiz yasalara göre bir elektronun, elektromanyetik dalgalar halinde yayarak enerji kaybetmesi gerekir. Ancak yukarıda da belirtildiği gibi elektron enerji kaybederse yörüngesinin yarıçapı küçülür. Sonuç olarak elektron sabit yarıçaplı bir yörüngede dönemez. Hesaplamalar, radyasyon nedeniyle yörünge yarıçapının azalması sonucunda elektronun saniyenin yüz milyonda biri kadar bir sürede çekirdeğe düşmesi gerektiğini gösteriyor. Bu sonuç, atomların kararlılığını gösteren günlük deneyimimizle keskin bir şekilde çelişmektedir.

Demek ki atomun yapısına ilişkin deneyden elde edilen veriler ile yine deneysel olarak bulunan mekanik ve elektriğin temel kanunları arasında bir çelişki vardır.

Ancak bahsedilen yasaların çok fazla sayıda elektron ve çok sayıda atom içeren cisimler üzerinde yapılan deneylerde bulunup test edildiğini de unutmamak gerekir. Bu yasaların atomdaki tek bir elektronun hareketine uygulanacağına inanmamız için hiçbir nedenimiz yok. Dahası, bir atomdaki elektronun davranışı ile yasalar arasındaki tutarsızlık klasik fizik bu yasaların atomik olaylara uygulanamayacağını gösterir (ayrıca bkz. § 210).

Yukarıda atomun sözde gezegen modelinin ana hatlarını çizdik; elektronların izin verilen yörüngelerde dönmesi fikri atom çekirdeği. Gezegen modelini gerekçelendirirken klasik fizik yasalarını kullandık. Ancak, daha önce de belirtildiği gibi ve § 210'da daha ayrıntılı olarak göreceğimiz gibi, bir atomdaki elektronun hareketi, içinde yer aldığı olgular alanına aittir. klasik mekanik uygulanamaz. Bu nedenle, "mikro dünya"nın daha derinlemesine incelenmesinin, gezegen modelinin eksikliğini ve kabaca yaklaşıklığını göstermesi şaşırtıcı değildir; atomun gerçek resmi daha karmaşıktır. Bununla birlikte, bu model atomun temel özelliklerinin çoğunu doğru bir şekilde yansıtmaktadır ve bu nedenle, yaklaşımına rağmen bazen kullanılmaktadır.

Bir hidrojen atomunun enerjisinin elektron yörüngesinin yarıçapına bağımlılığını düşünelim. Yarıçaplı bir yörünge boyunca elektron hareketinin kinetik enerjisini, merkezcil ivmenin yüklerin Coulomb çekim kuvveti tarafından sağlanması koşulundan (SI sisteminde) belirleriz. Bu kuvvetin yarattığı ivmeyi merkezcil ivmeye eşitlediğimizde elektronun kinetik enerjisinin yörünge yarıçapı ile ters orantılı olduğunu buluruz. .

Yarıçapı ve olan iki yörünge seçelim. İkinci yörüngedeki elektron dönüşünün kinetik enerjisi birinciden bir miktar daha fazladır .

Yörüngeler birbirinden çok uzak değilse, o zaman . Bu nedenle paydadaki miktar ihmal edilebilir ve kinetik enerjilerdeki fark yaklaşık olarak eşit olacaktır.

Elektronun potansiyel enerjisi ise tam tersine, ilk uzak yörüngede daha büyüktür, çünkü elektronu seriden çıkarmak için elektron ile çekirdek arasında etki eden elektriksel çekim kuvvetlerine karşı iş yapılması gerekir; bu çalışma potansiyel enerjiyi arttırmaya yöneliktir.

Bir elektronun yakın bir yörüngeden uzak bir yörüngeye radyal bir yol boyunca aktarılmasına izin verin. Yol uzunluğu . Elektrik gücü bu yol boyunca büyüklük olarak sabit değildir. Ancak yörüngeler birbirine yakın olduğundan, işin yaklaşık bir hesaplaması için, elektronun çekirdeğe ortalama mesafesindeki kuvvetin değerini kullanmak mümkündür: . Coulomb yasasına göre kuvvet ve yol üzerindeki iş, potansiyel enerjideki artışa eşit olacaktır.

Böylece bir elektron uzak bir yörüngeden yakın bir yörüngeye hareket ettiğinde potansiyel enerjisindeki azalma kinetik enerjisindeki artışın iki katına eşit olur. Bu teoremi, aralarındaki mesafenin koşulu karşıladığı yakın yörüngeler için kanıtladık. Ardışık yakın yörünge çiftleri arasındaki geçişler sırasında elektron enerjisindeki değişiklikleri toplayarak, teoremin keyfi olarak uzak yörüngeler için de geçerli olduğuna ikna olduk.

Şimdi sonsuz uzaklıktaki bir yörüngeyi ele alalım; Elektronun üzerindeki potansiyel enerjisini potansiyel enerjinin kaynağı olarak alalım, yani . Kinetik enerji sıfıra gider; bir yörüngeden yarıçaplı son bir yörüngeye geçerken bir miktar artacaktır. Potansiyel enerji miktarın iki katı kadar azalacaktır, yani.

.(206.1)

Bu nedenle elektronun toplam enerjisi şuna eşittir: ; Yörünge yarıçapı ne kadar küçükse, o kadar küçüktür (eksi işareti!).

Birkaç ay süren çalışmanın ardından Bohr, 1913'te atomun kuantum teorisini yayınladı. Bu teori üç varsayıma dayanmaktadır.

Bohr'un ilk varsayımı :

Bir atom, klasik fiziğin izin verdiği tüm durumlarda olmayabilir, ancak yalnızca her biri kendine özgü E n enerjisine sahip olan özel, kuantum (veya durağan) durumlarda olabilir. Durağan durumda bir atom enerji yaymaz veya absorbe etmez.

Saniye Bohr'un varsayımı:

Bir atom bir durağan durumdan diğerine geçtiğinde, durağan durumların enerjileri arasındaki farka eşit ћω enerjili bir ışık kuantumu yayılır veya emilir (Şekil 25.5):

ћω = |E n 2 -E n 1 |

(25.1)

E n 1 başlangıç ​​durumundaki enerjidir, E n 2 ise son durumdaki enerjidir. Bohr'un varsayımı:

Üçüncü

Sabit bir durumda, bir elektron yalnızca yarıçapı koşulu karşılayan böyle bir (“izin verilen”) yörünge boyunca hareket edebilir:

m·υ·r=n·ћ (25,2) Durağanlık koşulu elektron yörüngeleri

, burada m·υ·r elektronun açısal momentumudur, n ise kuantum durumunun sayısıdır (n =1, 2, 3, ...). Kuantum durumunun sayısını ve bu durumdaki atomun enerjisini belirleyen n tam sayısına denir. .

baş kuantum sayısı

Teorisini en basit atom olan hidrojen atomuna uygulayan Bohr, deneysel verilerle tamamen uyumlu sonuçlar elde etti. En basit atom olan hidrojen atomunu ele alalım. Çekirdeğin etrafında dairesel bir yörüngede dönen, bir proton ve bir elektron içeren bir çekirdekten oluşur. Elektron çekirdeğin yanından etkilenmektedir. Coulomb kuvveti

(25.3)

çekim, ona merkezcil ivme kazandırıyor, dolayısıyla

Bohr'un ilk varsayımının yerine getirilmesi gerektiğinden, elektron yörüngelerinin durağan durumunu kullanacağız. Buradan v hızını belirleyelim.

(25.4)

karesini alın ve (25.4) yerine koyun. Ortaya çıkan ifadeden bulduğumuz

dolayısıyla bir hidrojen atomundaki elektron yörüngelerinin yarıçapı şuna eşittir:

(25.5)

Sabitlerin değerlerini (25.5) yerine koyup n = 1 sayarak atom fiziğinde bir uzunluk birimi olan ilk Bohr yarıçapının değerini elde ederiz:

rB = 0,528-10-10 m.

§ 25.3 Hidrojen atomunun enerjisi

Bohr modeline göre atomun çekirdeği hareketsiz kabul edilir, bu nedenle toplam enerji Bir atomun E'si, elektronun dönüşünün kinetik enerjisi E k ile elektronun çekirdekle etkileşiminin potansiyel enerjisi E p'nin toplamıdır:


(25.6)

E'nin ortaya çıkan değeri negatiftir, çünkü sonsuz büyük bir mesafede bulunan iki yükün potansiyel enerjisi varsayılır. sıfıra eşit. Yükler birbirine yaklaştığında potansiyel enerji azalır.

Belirli bir atomun sahip olduğu enerjinin her değeri durağan durum, isminde enerji seviyesi . N ne kadar büyükse, elektron çekirdekten o kadar uzaktır ve enerji seviyesi de o kadar yüksek olur.

Bir atomun enerji seviyeleri genellikle yatay çizgilerle ve bir atomun bir durağan durumdan diğerine geçişleri oklarla gösterilir (Şekil 25.6).

Bir atom daha yüksek bir seviyeden daha düşük bir seviyeye hareket ettiğinde (bu, bir elektronun çekirdeğe daha yakın bir yörüngeye "sıçramasına" karşılık gelir), bir miktar ışık yayılır. Absorbsiyon sırasında ise tam tersine, bir atomun üzerine gelen bir kuantum (foton) olayı, atomu daha düşük enerjili bir durumdan daha yüksek enerjili bir duruma aktarır; fotonun kendisi kaybolur ve onu soğuran elektron, çekirdekten daha uzak bir yörüngeye yerleşir.

İLE atomun durumu сn = 1 olarak adlandırılır temel veya normal durum . Bu durumda, atomun enerjisi minimumdur ve içinde kalabilir (yokluğunda) dış etkiler) istenilen süre boyunca.

n>1 olan diğer tüm durumlara denir heyecanlı . Bir atom, çok kısa bir süre boyunca (yaklaşık 10-8 saniye) uyarılmış bir durumda kalabilir, ardından karşılık gelen kuantayı yayarak kendiliğinden (hemen veya kademeli olarak, seviye seviye) temel duruma geçer.

Temel durumda, hidrojen atomunun enerjisi E = -13,6 eV'dir. Uyarılmış durumlara geçerken enerjisi artar.

Bir elektronun birinci Bohr yörüngesinden “sonsuzluğa” çıkarılması için harcanması gereken minimum enerjiye denir iyonlaşma enerjisi W і veya bir hidrojen atomunun bağlanma enerjisi.

Bu nedenle temel durumdaki bir hidrojen atomunu iyonize etmek için ona ΔE enerjisi verilmesi gerekir. = Wi = 13,6 eV. Eğer ΔE enerjisi ona aktarılırsa < W i, o zaman ΔE = E n -E i olduğunda atom E p enerjili bir duruma girecek ve ΔE ≠ E n -E i olduğunda enerji emilimi gerçekleşmeyecek ve atom aynı durumda kalacaktır.

Enerji emiliminin bu (“sıçrayış benzeri”) doğası, herhangi bir kimyasal elementin atomları için gözlemlenmelidir. Cıva atomları için, 1913'te Alman deneysel fizikçiler D. Frank ve G. Hertz tarafından keşfedildi. Deneyleri, atomun kuantum teorisinin gelişmesinde çok önemli bir rol oynayan atomlarda ayrık enerji seviyelerinin varlığını doğruladı.



Makaleyi beğendin mi? Arkadaşlarınızla paylaşın!