Thursday 20 July 2017

Moving Average Window Matlab


Cara sederhana (ad hoc) adalah dengan hanya mengambil rata-rata tertimbang (merdu oleh alfa) pada setiap titik dengan tetangganya: atau beberapa variasi daripadanya. Ya, agar lebih canggih Anda bisa mengubah data Fourier terlebih dahulu, lalu memotong frekuensi tinggi. Sesuatu seperti: Ini memotong 20 frekuensi tertinggi. Hati-hati untuk memotongnya secara simetris sebaliknya transformasi terbalik tidak lagi nyata. Anda perlu hati-hati memilih frekuensi cutoff untuk tingkat smoothing yang tepat. Ini adalah jenis penyaringan yang sangat sederhana (penyaringan kotak pada domain frekuensi), sehingga Anda dapat mencoba dengan lembut menipiskan frekuensi orde tinggi jika distorsi tidak dapat diterima. Menjawab 4 Okt 09 at 9:16 FFT isnt ide yang buruk, tapi mungkin berlebihan di sini. Berjalan atau moving averages memberikan hasil yang umumnya buruk dan harus dihindari untuk segala hal selain PR akhir (dan white noise). Id menggunakan filter Savitzky-Golay (di Matlab sgolayfilt (.)). Ini akan memberi Anda hasil terbaik untuk apa yang Anda cari - beberapa smoothing lokal sambil mempertahankan bentuk kurva. Tutorial ini membahas bagaimana menggunakan MATLAB untuk pemrosesan gambar. Beberapa keakraban dengan MATLAB diasumsikan (Anda harus tahu bagaimana menggunakan matriks dan menulis file M). Akan sangat membantu jika memiliki Toolbox Pengolahan Gambar MATLAB, tapi untungnya, tidak ada kotak peralatan yang dibutuhkan untuk sebagian besar operasi. Perintah yang membutuhkan Image Toolbox ditunjukkan dengan Image Toolbox. Representasi gambar Ada lima jenis gambar dalam MATLAB. Grayscale. Gambar grayscale M piksel tinggi dan lebar N piksel direpresentasikan sebagai matriks tipe datatype ganda dengan ukuran M N. Nilai elemen (misalnya MyImage (m, n)) menunjukkan intensitas piksel grayscale dalam 0,1 dengan 0black dan 1white. RGB truecolor. Citra truecolor red-green-blue (RGB) diwakili sebagai matriks ganda M-3 tiga dimensi. Setiap piksel memiliki komponen merah, hijau, biru di sepanjang dimensi ketiga dengan nilai 0,1, misalnya komponen warna piksel (m, n) adalah MyImage (m, n, 1) merah, MyImage (m, n, 2) hijau, MyImage (m, n, 3) biru. Indeks. Gambar indeks (paletted) diwakili dengan matriks indeks dengan ukuran M N dan matriks colormap dengan ukuran K3. Colormap memegang semua warna yang digunakan pada gambar dan matriks indeks mewakili piksel dengan mengacu pada warna pada colormap. Misalnya, jika warna ke 22 adalah magenta MyColormap (22, :) 1,0,1. Maka MyImage (m, n) 22 adalah pixel. Binary berwarna magenta. Citra biner diwakili oleh matriks logika M N dimana nilai pikselnya adalah 1 (true) atau 0 (false).uint8. Jenis ini menggunakan sedikit memori dan beberapa operasi menghitung lebih cepat dibandingkan dengan tipe ganda. Untuk kesederhanaan, tutorial ini tidak membahas uint8 lebih jauh. Grayscale biasanya merupakan format yang disukai untuk pemrosesan gambar. Dalam kasus yang membutuhkan warna, gambar warna RGB dapat didekomposisi dan ditangani sebagai tiga gambar grayscale yang terpisah. Gambar terindeks harus dikonversi ke grayscale atau RGB untuk sebagian besar operasi. Berikut adalah beberapa manipulasi dan konversi yang umum. Beberapa perintah memerlukan Image Toolbox dan ditunjukkan dengan Image Toolbox. Membaca dan menulis file gambar MATLAB bisa membaca dan menulis gambar dengan perintah imft dan imwrite. Meskipun cukup banyak format file yang didukung, ada juga yang tidak. Gunakan imformats untuk melihat apa yang didukung oleh instalasi Anda: Saat membaca gambar, masalah yang tidak menguntungkan adalah imver mengembalikan data gambar pada tipe data uint8, yang harus dikonversi menjadi dua kali lipat dan dikompres ulang sebelum digunakan. Jadi alih-alih memanggil imread secara langsung, saya menggunakan fungsi M-file berikut untuk membaca dan mengubah gambar: Klik kanan dan simpan getimage. m untuk menggunakan fungsi M ini. Jika gambar baboon. png ada di direktori saat ini (atau di suatu tempat di jalur pencarian MATLAB), Anda dapat membacanya dengan MyImage getimage (baboon. png). Anda juga bisa menggunakan jalur parsial, misalnya jika gambar berada pada direktori saat ini dengan getimage (imagesbaboon. png). Untuk menulis gambar grayscale atau RGB, gunakan Berhati-hatilah bahwa MyImage adalah matriks ganda dengan elemen dalam 0,1 jika tidak bersesuaian dengan benar, file yang tersimpan mungkin akan kosong. Saat menulis file gambar, saya sangat menyarankan penggunaan format file PNG. Format ini adalah pilihan yang andal karena lossless, mendukung truecolor RGB, dan kompres dengan cukup baik. Gunakan format lain dengan hati-hati. Operasi dasar Berikut adalah beberapa operasi dasar pada gambar grayscale u. Perintah yang membutuhkan Image Toolbox ditunjukkan dengan Image Toolbox. (Catatan: Untuk array apapun, sintaks u (:) berarti membuka gulungan u ke dalam vektor kolom. Misalnya, jika u 1,50,2, maka u (:) adalah 1052.) Misalnya, kekuatan sinyal gambar digunakan pada Menghitung rasio signal-to-noise (SNR) dan rasio signal-to-noise (PSNR) puncak. Diberikan citra bersih uclean dan gambar yang terkontaminasi kebisingan u, hati-hati dengan norma. Perilaku adalah norma (v) pada vektor v menghitung sqrt (jumlah (v.2)). Tapi norma (A) pada matriks A menghitung norma matriks L 2 yang diinduksi, Jadi norma (A) tentu tidak sqrt (jumlah (A (:) 2)). Namun demikian merupakan kesalahan yang mudah untuk menggunakan norma (A) di mana seharusnya menjadi norma (A (:)). Linear filter Linear filtering adalah teknik dasar pemrosesan sinyal. Untuk memperkenalkan secara singkat, filter linier adalah operasi dimana pada setiap piksel x m, n dari sebuah gambar, fungsi linear dievaluasi pada piksel dan tetangganya untuk menghitung nilai piksel baru y m, n. Sebuah filter linear dalam dua dimensi memiliki bentuk umum dimana x adalah input, y adalah outputnya, dan h adalah respon impuls filter. Pilihan h yang berbeda mengarah ke filter yang menghaluskan, mempertajam, dan mendeteksi tepinya, untuk memberi beberapa aplikasi. Sisi kanan dari persamaan di atas dilambangkan secara ringkas seperti h x dan disebut konvolusi h dan x. Penyaringan domain spasial Penyaringan linear dua dimensi diterapkan pada MATLAB dengan konveks. Sayangnya, conv2 hanya bisa menangani penyaringan di dekat batas gambar dengan zero-padding, yang berarti hasil penyaringan biasanya tidak sesuai untuk piksel yang dekat dengan batas. Untuk mengatasi hal ini, kita dapat memberi masukan pada gambar masukan dan menggunakan opsi yang valid saat memanggil conv2. Fungsi M berikut melakukan ini. Klik kanan dan simpan conv2padded. m untuk menggunakan fungsi M ini. Berikut adalah beberapa contohnya: Filter 2D h dikatakan dapat dipisahkan jika dapat dinyatakan sebagai produk luar dari dua filter 1D h1 dan h2. Yaitu, h h1 (:) h2 (:). Lebih cepat lewat h1 dan h2 daripada h. Seperti yang dilakukan di atas untuk jendela rata-rata bergerak dan filter Gaussian. Sebenarnya, sobel filter hx dan hy juga separablewhat adalah h1 dan h2 Fourier-domain filtering Spatial-domain filtering with conv2 dengan mudah merupakan operasi yang sangat mahal. Untuk filter K K pada gambar M N, biaya tambahan dan multiplikasi V2 (MNK 2), atau O (N 4) seandainya M N K. Untuk filter besar, penyaringan di domain Fourier lebih cepat karena biaya komputasi dikurangi menjadi O (N 2 log N). Dengan menggunakan sifat perkalian konvolusi dari transformasi Fourier, konvolusi dihitung secara ekivalen oleh Hasilnya setara dengan conv2padded (x, h) kecuali di dekat batas, di mana perhitungan di atas menggunakan perpanjangan batas periodik. Penyaringan berbasis Fourier juga dapat dilakukan dengan perpanjangan batas simetris dengan memantulkan masukan ke setiap arah: (Catatan: Metode yang lebih efisien lagi adalah penyaringan tumpang tindih FFT. Toolbar Pemrosesan Sinyal menerapkan tumpang tindih FFT-satu dimensi dalam fftfilt .) Filter nonlinier Filter nonlinier adalah operasi dimana masing-masing piksel pixel disaring, n adalah fungsi nonlinear xm, n dan tetangganya. Disini kita secara singkat membahas beberapa jenis filter nonlinier. Filter statistik pesanan Jika Anda memiliki Image Toolbox, filter statistik pesanan dapat dilakukan dengan ordfilt2 dan medfilt2. Filter statistik pesanan memilah nilai piksel di atas lingkungan dan memilih nilai terbesar k. Filter min, max, dan median adalah kasus khusus. Filter morfologi Jika Anda memiliki Image Toolbox, bwmorph menerapkan berbagai operasi morfologi pada gambar biner, seperti erosi, pelebaran, terbuka, dekat, dan kerangka. Ada juga perintah yang tersedia untuk morfologi gambar grayscale: imerode. Imdilate dan imtophat. diantara yang lain. Buat filter Anda sendiri Terkadang kami ingin menggunakan filter baru yang tidak dimiliki MATLAB. Kode di bawah ini adalah template untuk menerapkan filter. (Catatan: Klaim yang sering salah arah adalah bahwa loop di MATLAB lambat dan harus dihindari. Ini pernah benar, kembali ke MATLAB 5 dan lebih awal, tapi loop dalam versi modern cukup cepat.) Misalnya, filter mean alpha-trimmed Mengabaikan d 2 terendah dan d 2 nilai tertinggi di jendela, dan rata-rata nilai sisa (2 r 1) 2 d. Filter adalah keseimbangan antara filter median dan filter rata-rata. Filter mean alpha-trimmed dapat diimplementasikan dalam template sebagai contoh lainnya, filter bilateralnya adalahGEOS 585A, Analisis Seri Terapan Waktu Telepon: (520) 621-3457 Faks: (520) 621-8229 Jam kerja Jumat, 1:00 -6: 00 PM (silahkan email ke pertemuan jadwal) Deskripsi Kursus Alat analisis dalam domain waktu dan frekuensi diperkenalkan dalam konteks seri waktu sampel. Saya menggunakan dataset dari seri waktu sampel untuk menggambarkan metode, dan mengubah dataset setiap semester kursus ditawarkan. Tahun ini dataset sampel berasal dari proyek NSF mengenai variabilitas snowpack di American River Basin of California. Dataset ini mencakup kronologi ring pohon, indeks iklim, catatan arus sungai, dan rangkaian waktu setara salju yang diukur di stasiun kursus salju. Anda akan mengumpulkan deret waktu Anda sendiri untuk digunakan dalam kursus. Ini mungkin berasal dari proyek penelitian Anda sendiri. Kembali ke Atas Halaman Ini adalah kursus pengantar, dengan penekanan pada aspek praktis dari analisis deret waktu. Metode diperkenalkan secara hierarkis - dimulai dengan grafis terminologi dan eksplorasi, beralih ke statistik deskriptif, dan diakhiri dengan prosedur pemodelan dasar. Topik meliputi detrending, filtering, autoregressive modeling, spektral analysis dan regression. Anda menghabiskan dua minggu pertama menginstal Matlab di laptop Anda, mendapatkan pengenalan dasar tentang Matlab, dan mengumpulkan dataset Anda untuk seri waktu kursus. Dua belas topik, atau pelajaran kemudian ditutup, masing-masing diberikan seminggu, atau dua periode kelas. Dua belas tugas kelas mengikuti topik. Penugasan terdiri dari penerapan metode dengan menjalankan skrip Matlab pra-tulis (program) pada deret waktu Anda dan menafsirkan hasilnya. Kursus 3 kredit untuk siswa di kampus di University of Arizona di Tucson, dan 1 kredit untuk siswa online. Setiap deret waktu dengan kenaikan waktu konstan (mis., Bulan, bulan, tahun) adalah kandidat untuk digunakan dalam kursus. Contohnya adalah pengukuran curah hujan setiap hari, aliran arus total musiman, suhu udara rata-rata musim panas, indeks pertumbuhan pohon tahunan, indeks suhu permukaan laut, dan kenaikan harian semak semak. Sebagai hasil dari mengikuti kursus, Anda harus: memahami konsep dan terminologi time series dasar dapat memilih metode time series yang sesuai dengan tujuan dapat mengevaluasi secara kritis literatur ilmiah yang menggunakan metode time series yang dibahas telah meningkatkan pemahaman tentang sifat deret waktu dari Dataset sendiri dapat ringkas merangkum hasil analisis deret waktu secara tertulis Prasyarat Kursus statistik pendahuluan Akses ke komputer laptop yang mampu menginstal Matlab di dalamnya Izin para instruktur (mahasiswa sarjana dan mahasiswa online) Persyaratan Lain Jika Anda berada di Universitas Mahasiswa Arizona (UA) di kampus di Tucson, Anda memiliki akses ke Matlab dan kotak peralatan yang dibutuhkan melalui lisensi situs UA karena tidak memerlukan perangkat lunak biaya. Tidak ada pengalaman sebelumnya dengan Matlab yang dibutuhkan, dan pemrograman komputer bukan bagian dari kursus. Jika Anda online, bukan di kampus UA, Anda akan bisa mengikuti kursus semester musim semi 2017 sebagai iCourse. Anda harus memastikan bahwa Anda memiliki akses ke Matlab dan kotak peralatan yang diperlukan (lihat di bawah) di lokasi Anda. Akses ke internet. Tidak ada pertukaran kertas dalam kursus. Catatan dan tugas ditukar secara elektronik dan selesai diserahkan secara elektronik melalui sistem University of Arizona Desire2Learn (D2L). Versi matlab Saya memperbarui skrip dan fungsi sekarang dan kemudian menggunakan rilis lisensi situs saat ini dari Matlab, dan pembaruannya mungkin menggunakan fitur Matlab yang tidak tersedia dalam rilis Matlab sebelumnya. Untuk 2017, saya menggunakan Matlab Version 9.1.0.441655 (R2016b). Jika Anda menggunakan rilis sebelumnya, pastikan itu Matlab Release 2007b atau lebih tinggi. Selain paket Matlab utama, empat toolboxes digunakan: Statistik, Pengolahan Sinyal, Identifikasi Sistem, dan Spline (Matlab Release 2010a atau sebelumnya), atau Curve Fitting (Matlab Release 2010b atau yang lebih baru) Ketersediaan Kursus ini ditawarkan di Semester Musim Semi Setiap tahun (2015, 2017, dst.). Ini terbuka untuk mahasiswa pascasarjana dan mungkin juga diambil oleh para manula senior dengan izin instruktur. Pendaftaran siswa UA tinggal ditutup pada usia 18 untuk Semester Musim Semi 2017. Sejumlah kecil siswa online juga biasanya diakomodasi dengan menawarkan kursus dengan berbagai cara. Caranya sekarang adalah tempat iCourse yang dijelaskan di atas. Kembali ke Atas Halaman Garis Besar Kursus (Pelajaran) Jadwal biasanya memungkinkan sekitar dua minggu untuk mengumpulkan data dan menjadi terbiasa dengan Matlab. Kemudian satu minggu (dua periode kelas) dikhususkan untuk masing-masing dari 12 pelajaran atau topik. Kelas bertemu pada hari Selasa dan Kamis. Topik baru diperkenalkan pada hari Selasa, dan dilanjutkan pada hari Kamis berikutnya. Kelas hari Kamis diakhiri dengan sebuah tugas dan demonstrasi menjalankan skrip pada data sampel saya. Tugasnya jatuh tempo (harus diunggah oleh Anda ke D2L) sebelum kelas pada hari Selasa berikutnya. 12 jam pertama kelas hari Selasa itu digunakan untuk penilaian diri yang dipandu dan penilaian tugas dan pengunggahan tugas dinilai (dinilai) ke D2L. Sisanya 45 menit digunakan untuk mengenalkan topik selanjutnya. Anda harus membawa laptop Anda ke kelas pada hari Selasa. 12 pelajaran atau topik yang dibahas dalam kursus tercantum dalam garis besar kelas. Siswa online diharapkan mengikuti jadwal penyerahan tugas yang sama dengan siswa yang tinggal, namun tidak memiliki akses ke ceramah. Tugas yang dikirim dari siswa online tidak dinilai sendiri, namun dinilai oleh saya. Siswa online harus memiliki akses ke D2L untuk mengirimkan tugas. Semester musim semi 2017 Kelas bertemu dua kali seminggu selama 75 menit, 9: 00-10: 15 AM TTh, di kamar 424 (Ruang Konferensi) Gedung Cincin Pohon Bryant Bannister (bangunan 45B). Hari pertama kelas adalah 12 Januari (Kam). Hari terakhir kelas adalah 2 Mei (sel). Tidak ada kelas selama minggu Spring Break (Mar 11-19). Anda menganalisis data pilihan Anda sendiri di kelas tugas. Sebagaimana tercantum dalam ikhtisar kursus. Ada banyak fleksibilitas dalam pemilihan deret waktu. Saya akan membuat katalog rangkaian waktu yang sesuai, tapi yang terbaik adalah memfokuskan kursus pada kumpulan data Anda sendiri. Tugas pertama melibatkan menjalankan skrip yang menyimpan data dan metadata yang telah Anda kumpulkan di file mat, format asli Matlab. Tugas selanjutnya menarik data dari file mat untuk analisis deret waktu. Penugasan 12 topik tersebut dibahas secara berurutan sepanjang semester, yang mencakup sekitar 15 minggu. Tentang dua minggu pertama (pertemuan kelas 4-5) digunakan untuk beberapa bahan pengantar, menentukan dan mengumpulkan deret waktu Anda, dan menyiapkan Matlab di laptop Anda. Setiap minggu setelah itu dikhususkan untuk salah satu dari 12 topik topik. Setiap tugas terdiri dari membaca bab catatan, menjalankan skrip Matlab terkait yang menerapkan metode analisis time series pilihan ke data Anda, dan menuliskan interpretasi Anda terhadap hasilnya. Tugas memerlukan pemahaman tentang topik kuliah serta kemampuan untuk menggunakan komputer dan perangkat lunak. Anda mengirimkan tugas dengan mengunggahnya ke D2L sebelum kelas Selasa saat topik berikutnya diperkenalkan. Semester pertama kelas Selasa itu digunakan untuk penilaian diri yang dipandu oleh penugasan, termasuk mengunggah PDF dengan self-grade ke D2L. Saya memeriksa satu atau beberapa tugas yang dinilai sendiri setiap minggu (dengan seleksi acak), dan mungkin mengubah nilainya. Untuk mengetahui cara mengakses tugas, klik file tugas. Bacaan terdiri dari catatan kelas. Ada dua belas set. pdf mencatat file. Satu untuk masing-masing topik kursus. File. pdf ini dapat diakses melalui Web. Informasi lebih lanjut tentang berbagai topik yang dibahas dalam kursus dapat ditemukan melalui referensi yang tercantum di akhir setiap bab catatan kelas. Kelas didasarkan sepenuhnya pada kinerja pada tugas, masing-masing bernilai 10 poin. Tidak ada ujian. Jumlah total poin yang mungkin untuk 12 topik adalah 12 x 10 120. Nilai A yang dibutuhkan 90-100 persen dari poin yang mungkin. Nilai B membutuhkan 80-90 persen. Nilai C membutuhkan 70-80 persen, dan sebagainya. Nilai diberikan dengan penilaian diri yang dipandu oleh rubrik yang disajikan di kelas. Jumlah poin yang diterima harus ditandai di bagian atas setiap tugas bergradasi. Markup penugasan Anda harus menyertakan anotasi dari setiap penurunan harga dengan mengacu pada rubrik yang diilustrasikan di kelas (misalnya -0,5, rp3 menunjukkan pengurangan sebesar -0,5 karena kesalahan yang terkait dengan rubrik poin 3) Tugas, diberikan di kelas pada hari Kamis, akan Karena (diunggah ke D2L oleh Anda) sebelum memulai kelas pada hari Selasa berikutnya. Setengah jam pertama periode pertemuan hari Selasa akan didedikasikan untuk presentasi rubrik penilaian, penilaian sendiri atas penugasan yang telah selesai, dan pengunggahan tugas yang dinilai sendiri ke D2L. Jadwal ini memberi Anda waktu 4 hari untuk menyelesaikan dan mengunggah tugas ke D2L sebelum pukul 09:00 hari Selasa. D2L melacak waktu penugasan diupload, dan tidak ada hukuman yang dinilai selama diunggah sebelum pukul 09:00 pada hari Selasa tanggal jatuh tempo. Jika Anda memiliki beberapa jadwal yang harus jauh dari kelas (misalnya, kehadiran di sebuah konferensi), Anda bertanggung jawab untuk mengunggah tugas sebelum pukul 09:00 hari Selasa karena waktunya, dan untuk mengupload versi self-graded pada pukul 10:15 pagi. hari yang sama. Dengan kata lain, jadwalnya sama dengan siswa yang berada di kelas. Jika keadaan darurat muncul (misalnya Anda terkena flu) dan tidak dapat melakukan tugas atau penilaian sesuai jadwal, kirimkan saya email dan kami akan sampai di akomodasi. Jika tidak, denda 5 poin (setengah dari total poin yang tersedia untuk latihan) akan dinilai. Pengenalan data pengorganisasian rangkaian waktu untuk analisis Suatu deret waktu didefinisikan secara luas sebagai serangkaian pengukuran yang dilakukan pada waktu yang berbeda. Beberapa kategori deskriptif dasar deret waktu adalah 1) panjang vs pendek, 2) bahkan langkah waktu vs langkah waktu yang tidak rata, 3) diskrit vs kontinyu, 4) periodik vs aperiodik, 5) stasioner vs nonstasioner, dan 6) univariat vs multivariat . Sifat-sifat ini dan juga tumpang tindih temporal dari beberapa seri, harus dipertimbangkan dalam memilih kumpulan data untuk analisis dalam kursus ini. Anda akan menganalisis rangkaian waktu Anda sendiri di kursus. Langkah pertama adalah memilih seri tersebut dan menyimpannya dalam struktur di file tikar. Keseragaman dalam penyimpanan pada awalnya sangat sesuai untuk kelas ini sehingga perhatian kemudian dapat difokuskan pada pemahaman metode deret waktu, bukan debug kode komputer untuk menyiapkan data untuk analisis. Struktur adalah variabel Matlab yang mirip dengan database sehingga isinya diakses oleh penanda lapangan tekstual. Struktur dapat menyimpan data dari berbagai bentuk. Sebagai contoh, satu bidang mungkin merupakan matriks deret waktu numerik, yang lain mungkin berupa teks yang menjelaskan sumber data, dsb. Dalam tugas pertama Anda akan menjalankan skrip Matlab yang membaca rangkaian waktu dan metadata Anda dari file teks ascii yang Anda siapkan sebelumnya dan Menyimpan data di struktur Matlab dalam file matrik tunggal. Dalam tugas selanjutnya Anda akan menerapkan metode time series ke data dengan menjalankan skrip dan fungsi Matlab yang memuat file mat dan mengoperasikan struktur tersebut. Pilih data sampel yang akan digunakan untuk tugas selama kursus Baca: (1) Notes1.pdf, (2) Persiapan, dapat diakses dari menu bantuan MATLAB Jawab: Jalankan skrip geosa1.m dan jawablah pertanyaan yang tercantum dalam file di a1.pdf Bagaimana membedakan kategori deret waktu Bagaimana cara memulai dan berhenti MATLAB Bagaimana cara memasukkan perintah MATLAB pada command prompt Bagaimana membuat angka di jendela gambar Bagaimana cara mengekspor tokoh ke pengolah kata Anda Perbedaan antara skrip dan fungsi MATLAB Bagaimana cara menjalankan skrip dan fungsi Bentuk variabel struktur MATLAB Bagaimana menerapkan skrip geosa1.m untuk mendapatkan serangkaian rangkaian waktu dan metadata ke dalam struktur MATLAB Distribusi probabilitas deret waktu menggambarkan probabilitas bahwa pengamatan masuk ke dalam kisaran nilai tertentu. Distribusi probabilitas empiris untuk rangkaian waktu dapat dicapai dengan memilah dan memberi peringkat nilai dari seri. Quantiles dan persentil adalah statistik yang berguna yang dapat diambil secara langsung dari distribusi probabilitas empiris. Banyak uji statistik parametrik mengasumsikan deret waktu adalah sampel dari populasi dengan distribusi probabilitas populasi tertentu. Seringkali penduduk dianggap normal. Bab ini menyajikan beberapa definisi dasar, statistik dan plot yang terkait dengan distribusi probabilitas. Sebagai tambahan, sebuah tes (uji Lilliefors) diperkenalkan untuk menguji apakah sampel berasal dari distribusi normal dengan mean dan varians yang tidak ditentukan. Jawaban: Jalankan skrip geosa2.m dan jawab pertanyaan yang tercantum dalam file di a2.pdf Definisi istilah: deret waktu, stasioneritas, kepadatan probabilitas, fungsi distribisi, quantile, spread, lokasi, mean, standar deviasi, dan condong Bagaimana menafsirkan Grafik paling berharga dalam analisis deret waktu - deret seri waktu Bagaimana menafsirkan kotak petak, histogram dan plot probabilitas normal Parameter dan bentuk dari distribusi normal Uji Lilliefors untuk normalitas: deskripsi grafis, asumsi, hipotesis nol dan alternatif Peringatan pada interpretasi Tingkat signifikansi uji statistik ketika deret waktu tidak acak dalam waktu Bagaimana menerapkan geosa2.m untuk memeriksa properti distribusi dari deret waktu dan menguji seri untuk normalitas Autokorelasi mengacu pada korelasi deret waktu dengan nilai masa lalu dan masa depannya sendiri. Autokorelasi juga kadang disebut korelasi tertinggal atau korelasi serial. Yang mengacu pada korelasi antara anggota dari serangkaian angka yang disusun pada waktunya. Autokorelasi positif bisa dianggap sebagai bentuk ketekunan yang spesifik. Kecenderungan sebuah sistem untuk tetap berada dalam keadaan yang sama dari satu pengamatan ke pengamatan berikutnya. Misalnya, kemungkinan besok hujan lebih besar jika hari ini hujan daripada jika hari ini kering. Seri waktu geofisika sering kali autokorelasi karena proses inersia atau carryover dalam sistem fisik. Misalnya, sistem tekanan rendah yang berkembang perlahan dan bergerak di atmosfer bisa memberi ketekunan pada curah hujan harian. Atau drainase yang lambat dari cadangan air tanah mungkin memberi korelasi dengan arus tahunan sungai yang berturut-turut. Atau fotosintat yang tersimpan mungkin memberi korelasi dengan nilai tahunan indeks cincin-pohon berturut-turut. Autokorelasi mempersulit penerapan uji statistik dengan mengurangi jumlah pengamatan independen. Autokorelasi juga dapat mempersulit identifikasi kovariansi signifikan atau korelasi antara deret waktu (misalnya presipitasi dengan deret pohon). Autokorelasi dapat dieksploitasi untuk prediksi: rangkaian waktu autokorelasi dapat diprediksi, probabilistik, karena nilai masa depan tergantung pada nilai arus dan masa lalu. Tiga alat untuk menilai autokorelasi deret waktu adalah (1) rangkaian deret waktu, (2) scatterplot yang tertinggal, dan (3) fungsi autokorelasi. Jawaban: Jalankan skrip geosa3.m dan jawab pertanyaan yang tercantum dalam file di a3.pdf Definisi: autokorelasi, ketekunan, korelasi serial, fungsi autokorelasi (acf), fungsi autocovariance (acvf), ukuran sampel efektif Bagaimana mengenali autokorelasi dalam deret waktu Plot Bagaimana menggunakan scatterplots yang tertinggal untuk menilai autokorelasi Bagaimana menafsirkan acf diplot Bagaimana menyesuaikan ukuran sampel untuk autokorelasi Definisi matematis dari fungsi autokorelasi Persyaratan yang mempengaruhi lebar pita kepercayaan dihitung dari acf Perbedaan antara satu sisi dan dua - dari uji autokorelasi lag-1 yang signifikan Bagaimana menerapkan geos3.m untuk mempelajari autokorelasi deret waktu Spektrum deret waktu adalah distribusi varians rangkaian sebagai fungsi frekuensi. Objek analisis spektral adalah untuk memperkirakan dan mempelajari spektrum. Spektrum tidak mengandung informasi baru selain fungsi autocovariance (acvf), dan kenyataannya spektrumnya dapat dihitung secara matematis dengan transformasi acvf. Tapi spektrum dan acvf menyajikan informasi tentang varians deret waktu dari sudut pandang komplementer. Acf merangkum informasi dalam domain waktu dan spektrum dalam domain frekuensi. Jawaban: Jalankan skrip geosa4.m dan jawab pertanyaan yang tercantum dalam file di a4.pdf Definisi: frekuensi, periode, panjang gelombang, spektrum, frekuensi Nyquist, frekuensi Fourier, bandwidth Alasan untuk menganalisis spektrum Bagaimana menafsirkan spektrum diplot dalam hal distribusi Varians Perbedaan antara spektrum dan spektrum normal Definisi jendela lag seperti yang digunakan dalam memperkirakan spektrum dengan metode Blackman-Tukey Bagaimana pilihan jendela lag mempengaruhi bandwidth dan varians spektrum perkiraan Bagaimana menentukan spektrum suara putih Dan spektrum autoregresif Bagaimana membuat sketsa beberapa bentuk spektral yang khas: white noise, autoregressive, quasi-periodic, frekuensi rendah, frekuensi tinggi Bagaimana cara menerapkan geosa4.m untuk menganalisis spektrum deret waktu dengan metode Blackman-Tukey Autoregressive-Moving Model rata-rata (ARMA) model Autoregressive-moving-average (ARMA) adalah model matematis dari ketekunan, atau autokorelasi, dalam deret waktu. Model ARMA banyak digunakan dalam hidrologi, dendrochronologi, ekonometri, dan bidang lainnya. Ada beberapa kemungkinan alasan pemasangan model ARMA pada data. Pemodelan dapat berkontribusi untuk memahami sistem fisik dengan mengungkapkan sesuatu tentang proses fisik yang membangun ketekunan ke dalam rangkaian. Sebagai contoh, model keseimbangan air fisik sederhana yang terdiri dari istilah untuk input presipitasi, penguapan, infiltrasi, dan penyimpanan air tanah dapat ditunjukkan untuk menghasilkan rangkaian aliran arus yang mengikuti bentuk model ARMA tertentu. Model ARMA juga bisa digunakan untuk memprediksi perilaku deret waktu dari nilai masa lalu saja. Prediksi tersebut dapat digunakan sebagai dasar untuk mengevaluasi kemungkinan kemungkinan variabel lain terhadap sistem. Model ARMA banyak digunakan untuk prediksi deret waktu ekonomi dan industri. Model ARMA juga bisa digunakan untuk menghilangkan ketekunan. Dalam dendrochronology, misalnya, pemodelan ARMA diterapkan secara rutin untuk menghasilkan kronik waktu residu indeks ring-width tanpa ketergantungan pada nilai masa lalu. Operasi ini, yang disebut prewhitening, dimaksudkan untuk menghilangkan kegigihan yang terkait secara biologis dari rangkaian sehingga residu lebih sesuai untuk mempelajari pengaruh iklim dan faktor lingkungan luar lainnya terhadap pertumbuhan pohon. Jawaban: Jalankan skrip geosa5.m dan jawablah pertanyaan yang tercantum dalam file di a5.pdf Bentuk fungsional model AR dan ARMA yang paling sederhana Mengapa model seperti itu disebut sebagai autoregressive atau moving average Tiga langkah dalam pemodelan ARMA Pola diagnostik dari Autokorelasi dan fungsi autokorelasi parsial untuk rangkaian waktu AR (1) Definisi kesalahan prediksi akhir (FPE) dan bagaimana FPE digunakan untuk memilih model ARMA terbaik Definisi statistik Portmanteau, dan bagaimana dan residu residu dapat Digunakan untuk menilai apakah model ARMA secara efektif memodelkan ketekunan dalam rangkaian Bagaimana prinsip parsimoni diterapkan dalam pemodelan ARMA Definisi prewhitening Bagaimana pengaruh sebelum perang mempengaruhi (1) kemunculan deret waktu, dan (2) spektrum deret waktu Bagaimana menerapkan geosa5.m ke ARMA-model rangkaian waktu Analisis spektral - Metode periodogram merapikan Ada banyak metode yang tersedia untuk memperkirakan spektrum deret waktu. Dalam pelajaran 4 kita melihat metode Blackman-Tukey, yang didasarkan pada transformasi Fourier dari fungsi autocovariance yang merapikan dan dipotong. Metode periodogram merapikan mengeliminasi transformasi acf dengan transformasi Fourier langsung dari deret waktu dan perhitungan periodogram mentah, sebuah fungsi yang pertama kali diperkenalkan pada tahun 1800 untuk mempelajari deret waktu. Periodogram mentah diratakan dengan menerapkan kombinasi atau rentang satu atau lebih filter untuk menghasilkan spektrum yang diperkirakan. Kelancaran, resolusi dan varians perkiraan spektral dikendalikan oleh pilihan filter. Pemulusan periodogram baku yang lebih ditekankan menghasilkan spektrum yang bervariasi, atau kontinum null yang mendasari, yang dengannya puncak spektral dapat diuji signifikansinya. Pendekatan ini adalah alternatif dari spesifikasi bentuk fungsional dari kontinum null (misalnya spektrum AR). Jawaban: Jalankan skrip geosa6.m dan jawab pertanyaan yang tercantum dalam file di a6.pdf Definisi: periodogram mentah, filter Daniell, rentang filter, kelancaran kontinuitas null, stabilitas dan resolusi spektrum meruncing, padding, kebocoran Empat langkah utama dalam memperkirakan Spektrum oleh periodogram yang merapikan Bagaimana pengaruh pilihan bentang filter pada kelancaran, stabilitas dan resolusi spektrum Bagaimana kontinum null digunakan dalam pengujian untuk kepentingan puncak spektral Bagaimana menerapkan geosa6.m untuk memperkirakan spektrum suatu waktu Seri dengan metode periodogram merapikan dan uji periodisitas pada frekuensi tertentu Tren dalam deret waktu adalah perubahan bertahap dan lambat dalam beberapa properti seri selama keseluruhan interval yang sedang diselidiki. Trend kadang-kadang didefinisikan secara longgar sebagai perubahan jangka panjang dalam mean (Gambar 7.1), namun juga dapat merujuk pada perubahan pada sifat statistik lainnya. Misalnya, rangkaian cincin pohon dari lebar cincin yang diukur sering memiliki kecenderungan yang berbeda dan juga mean (Gambar 7.2). Dalam analisis deret waktu tradisional, deret waktu didekomposisi menjadi tren, komponen musiman atau periodik, dan fluktuasi yang tidak teratur, dan berbagai bagian dipelajari secara terpisah. Teknik analisis modern sering memperlakukan seri tanpa dekomposisi rutin seperti itu, namun pertimbangan tren yang terpisah masih sering dibutuhkan. Detrending adalah operasi statistik atau matematis untuk menghilangkan tren dari rangkaian. Detrending sering diterapkan untuk menghilangkan fitur yang diduga mendistorsi atau mengaburkan hubungan yang diminati. Dalam klimatologi, misalnya, tren suhu akibat pemanasan kota mungkin mengaburkan hubungan antara keruh dan suhu udara. Detrending juga kadang-kadang digunakan sebagai langkah preprocessing untuk mempersiapkan time series untuk analisis dengan metode yang mengasumsikan stationarity. Banyak metode alternatif tersedia untuk detrending. Tren linier sederhana dalam mean dapat dihapus dengan mengurangkan garis lurus kuadrat terkecil. Tren yang lebih rumit mungkin memerlukan prosedur yang berbeda. Sebagai contoh, spline smoothing kubik biasanya digunakan dalam dendrochronology agar sesuai dan menghilangkan tren ring-width yang mungkin tidak linier, atau bahkan tidak meningkat secara monoton atau menurun seiring berjalannya waktu. Dalam mempelajari dan menghilangkan kecenderungan, penting untuk memahami efek detrending pada sifat spektral dari deret waktu. Efek ini dapat diringkas dengan respon frekuensi fungsi detrending. Jawaban: Jalankan skrip geosa7.m dan jawab pertanyaan yang tercantum dalam file di a7.pdf Definisi: respons frekuensi, spline, spline smoothing kubik Pro dan kontra rasio vs perbedaan detrending Interpretasi istilah dalam persamaan untuk parameter spline Bagaimana memilih Spline secara interaktif dari respons frekuensi yang diinginkan Bagaimana spektrum dipengaruhi oleh detrending Bagaimana mengukur pentingnya komponen tren dalam deret waktu Bagaimana menerapkan geosa7.m untuk secara interaktif memilih fungsi detrending spline dan detrend time series Perkiraan spektrum suatu waktu Series memberikan distribusi varians sebagai fungsi frekuensi. Bergantung pada tujuan analisis, beberapa frekuensi mungkin lebih menarik daripada yang lain, dan ini mungkin berguna untuk mengurangi amplitudo variasi pada frekuensi lain dengan menyaringnya secara statistik sebelum melihat dan menganalisis rangkaian. Misalnya, variasi frekuensi tinggi (dari tahun ke tahun) dalam catatan debit terukur dari daerah aliran sungai mungkin relatif tidak penting untuk persediaan air di baskom dengan waduk besar yang dapat menyimpan beberapa tahun limpasan rata-rata tahunan. Bila variasi frekuensi rendah menjadi perhatian utama, diharapkan untuk memperlancar catatan debit untuk menghilangkan atau mengurangi fluktuasi periode pendek sebelum menggunakan catatan debit untuk mempelajari pentingnya variasi iklim terhadap persediaan air. Smoothing adalah bentuk penyaringan yang menghasilkan deret waktu di mana pentingnya komponen spektral pada frekuensi tinggi berkurang. Insinyur listrik menyebut filter jenis filter low-pass ini, karena variasi frekuensi rendah diperbolehkan melewati filter. Pada filter low-pass, frekuensi rendah (periode lama) ombak hampir tidak terpengaruh oleh smoothing. Hal ini juga memungkinkan untuk menyaring rangkaian sedemikian rupa sehingga variasi frekuensi rendah berkurang dan variasi frekuensi tinggi tidak terpengaruh. Filter jenis ini disebut filter high-pass. Detrending adalah bentuk high-pass filtering: garis tren yang dipasang melacak frekuensi terendah, dan residu dari garis tren memiliki frekuensi rendah yang dilepaskan. Jenis penyaringan ketiga, yang disebut penyaringan band-pass, mengurangi atau menyaring frekuensi tinggi dan rendah, dan meninggalkan beberapa pita frekuensi menengah yang relatif tidak terpengaruh. Dalam pelajaran ini, kita membahas beberapa metode perataan, atau penyaringan low-pass. Kita sudah membahas bagaimana spline smoothing kubik mungkin berguna untuk tujuan ini. Empat jenis filter lainnya dibahas di sini: 1) simple moving average, 2) binomial, 3) Gaussian, dan 4) windowing (metode Hamming). Pertimbangan dalam memilih jenis filter low-pass adalah respons frekuensi yang diinginkan dan rentang, atau lebar filter. Jawaban: Jalankan skrip geosa8.m dan jawab pertanyaan yang tercantum dalam file di a8.pdf Definisi: filter, bobot filter, rentang filter, filter low-pass, filter high-pass, tanggapan frekuensi filter band-pass filter Bagaimana Gaussian Filter berhubungan dengan distribusi Gaussian Bagaimana membangun sebuah filter binomial sederhana secara manual (tanpa komputer) Bagaimana mendeskripsikan fungsi respons frekuensi dalam hal sistem dengan input dan output sinusoidal Bagaimana menerapkan geosa8.m untuk merancang secara interaktif Gaussian, binomial Atau filter lowpass Hamming-window untuk rangkaian waktu Koefisien korelasi product moment Pearson mungkin merupakan statistik tunggal yang paling banyak digunakan untuk meringkas hubungan antara dua variabel. Signifikansi statistik dan peringatan interpretasi koefisien korelasi sebagaimana diterapkan pada deret waktu adalah topik pelajaran ini. Dengan asumsi tertentu, signifikansi statistik dari koefisien korelasi bergantung hanya pada ukuran sampel, yang didefinisikan sebagai jumlah observasi independen. Jika deret waktu diautokorelasi, ukuran sampel efektif, lebih rendah dari ukuran sampel sebenarnya, harus digunakan saat mengevaluasi signifikansi. Hubungan transien atau palsu dapat menghasilkan korelasi yang signifikan untuk beberapa periode dan bukan untuk yang lain. Variasi waktu kekuatan korelasi linier dapat diperiksa dengan plot korelasi yang dihitung untuk jendela geser. Tetapi jika banyak koefisien korelasi dievaluasi secara bersamaan, interval kepercayaan harus disesuaikan (penyesuaian Bonferroni) untuk mengkompensasi kemungkinan peningkatan pengamatan beberapa korelasi tinggi dimana tidak ada hubungan. Interpretasi korelasi geser juga dapat dipersulit oleh variasi waktu mean dan varians dari seri, karena korelasi geser mencerminkan kovariat dalam hal penyampaian standar dari mean pada jendela waktu yang diminati, yang mungkin berbeda dari mean jangka panjang. Akhirnya, perlu ditekankan bahwa koefisien korelasi Pearson mengukur kekuatan hubungan linier. Scatterplots berguna untuk memeriksa apakah hubungan itu linier. Jawaban: Jalankan skrip geosa9.m ​​dan jawab pertanyaan yang tercantum dalam file di a9.pdf Definisi matematika dari koefisien korelasi Asumsi dan hipotesis untuk uji signifikansi koefisien korelasi Bagaimana menghitung tingkat signifikansi koefisien korelasi dan untuk menyesuaikan tingkat signifikansi untuk autokorelasi dalam Seri waktu individu Peringatan untuk interpretasi koefisien korelasi Bonferroni penyesuaian terhadap tingkat signficance korelasi di bawah beberapa perbandingan Inflasi varians estimasi koefisien korelasi ketika deret waktu autokorelasi Kemungkinan efek transformasi data pada korelasi Bagaimana menafsirkan plot dari korelasi geser Bagaimana cara menerapkan geosa9. M untuk menganalisis korelasi dan korelasi geser antara pasangan deret waktu Hubungan tertinggal adalah karakteristik dari banyak sistem fisik alami. Korelasi tertunda mengacu pada korelasi antara dua deret waktu yang bergeser dalam waktu relatif terhadap satu sama lain. Korelasi tertunda penting dalam mempelajari hubungan antara deret waktu karena dua alasan. Pertama, satu seri mungkin memiliki respons tertunda terhadap seri lainnya, atau mungkin respons tertunda terhadap stimulus umum yang mempengaruhi kedua seri. Kedua, respon dari satu seri ke rangkaian lainnya atau stimulus luar dapat diolesi pada waktunya, sehingga stimulus yang dibatasi pada satu pengamatan menghasilkan respons pada beberapa pengamatan. Misalnya, karena penyimpanan di waduk, gletser, dan lain-lain, pelepasan volume sungai dalam satu tahun mungkin bergantung pada curah hujan pada beberapa tahun sebelumnya. Atau karena perubahan pada kepadatan mahkota dan penyimpanan fotosintat, lebar cincin pohon dalam satu tahun mungkin bergantung pada iklim beberapa tahun sebelumnya. Koefisien korelasi sederhana antara kedua seri yang selaras dengan tepat waktu tidak mencukupi untuk mengkarakterisasi hubungan dalam situasi seperti itu. Fungsi yang berguna yang akan kita bahas sebagai alternatif untuk koefisien korelasi sederhana adalah fungsi cross-correlation dan fungsi respon impuls. Fungsi cross-correlation adalah korelasi antara seri bergeser satu sama lain sebagai fungsi dari jumlah pengamatan offset. Jika seri individu diautokorelasi, fungsi korelasi silang diperkirakan dapat terdistorsi dan menyesatkan sebagai ukuran hubungan yang tertinggal. Kita akan melihat dua pendekatan untuk mengklarifikasi pola korelasi silang. Salah satunya adalah secara individual menghilangkan kegigihan dari, atau prewhiten, seri sebelum estimasi korelasi silang. Dalam pendekatan ini, kedua seri pada dasarnya dianggap setara. Alternatifnya adalah pendekatan sistem: lihat seri sebagai sistem linier dinamis - satu seri input dan output lainnya - dan perkirakan fungsi respons impuls. Fungsi respon impuls adalah respon output pada saat ini dan masa depan terhadap pulsa input hipotetis yang dibatasi pada waktu saat ini. Answer: Run script geosa10.m and answer questions listed in the file in a10.pdf Definitions: cross-covariance function, cross-correlation function, impulse response function, lagged correlation, causal, linear How autocorrelation can distort the pattern of cross-correlations and how prewhitening is used to clarify the pattern The distinction between the equal footing and systems approaches to lagged bivariate relationships Which types of situations the impulse response function (irf) is an appropriate tool How to represent the causal system treated by the irf in a flow diagram How to apply geos10.m to analyze the lagged cross-correlation structure of a a pair of time series Multiple linear regression Multiple linear regression (MLR) is a method used to model the linear relationship between a dependent variable and one or more independent variables. The dependent variable is sometimes also called the predictand, and the independent variables the predictors. MLR is based on least squares: the model is fit such that the sum-of-squares of differences of observed and predicted values is minimized. MLR is probably the most widely used method in dendroclimatology for developing models to reconstruct climate variables from tree-ring series. Typically, a climatic variable is defined as the predictand and tree-ring variables from one or more sites are defined as predictors. The model is fit to a period -- the calibration period -- for which climatic and tree-ring data overlap. In the process of fitting, or estimating, the model, statistics are computed that summarize the accuracy of the regression model for the calibration period. The performance of the model on data not used to fit the model is usually checked in some way by a process called validation. Finally, tree-ring data from before the calibration period are substituted into the prediction equation to get a reconstruction of the predictand. The reconstruction is a prediction in the sense that the regression model is applied to generate estimates of the predictand variable outside the period used to fit the data. The uncertainty in the reconstruction is summarized by confidence intervals, which can be computed by various alternative ways. Answer: Run script geosa11.m (Part 1) and answer questions listed in the file in a11.pdf The equation for the MLR model Assumptions for the MLR model Definitions of MLR statistics: coefficient of determination, sums-of-squares terms, overall-F for the regression equation, standard error of the estimate, adjusted R-squared, pool of potential predictors The steps in an analysis of residuals How to apply geosa11.m (part 1) to fit a MLR regression model to predict one variable from a set of several predictor variables Validating the regression model Regression R-squared, even if adjusted for loss of degrees of freedom due to the number of predictors in the model, can give a misleading, overly optimistic view of accuracy of prediction when the model is applied outside the calibration period. Application outside the calibration period is the rule rather than the exception in dendroclimatology. The calibration-period statistics are typically biased because the model is tuned for maximum agreement in the calibration period. Sometimes too large a pool of potential predictors is used in automated procedures to select final predictors. Another possible problem is that the calibration period itself may be anomalous in terms of the relationships between the variables: modeled relationships may hold up for some periods of time but not for others. It is advisable therefore to validate the regression model by testing the model on data not used to fit the model. Several approaches to validation are available. Among these are cross-validation and split-sample validation. In cross-validation, a series of regression models is fit, each time deleting a different observation from the calibration set and using the model to predict the predictand for the deleted observation. The merged series of predictions for deleted observations is then checked for accuracy against the observed data. In split-sample calibration, the model is fit to some portion of the data (say, the second half), and accuracy is measured on the predictions for the other half of the data. The calibration and validation periods are then exchanged and the process repeated. In any regression problem it is also important to keep in mind that modeled relationships may not be valid for periods when the predictors are outside their ranges for the calibration period: the multivariate distribution of the predictors for some observations outside the calibration period may have no analog in the calibration period. The distinction of predictions as extrapolations versus interpolations is useful in flagging such occurrences. Answer: Run script geosa11.m (Part 2) and answer questions listed in the file in a12.pdf Definitions: validation, cross-validation, split-sample validation, mean square error (MSE), root-mean-square error (RMSE) standard error of prediction, PRESS statistic, hat matrix, extrapolation vs interpolation Advantages of cross-validation over alternative validation methods How to apply geosa11.m (part 2) for cross-validated MLR modeling of the relationship between a predictand and predictors, including generation of a reconstruction and confidence bands Downloading Files -- tsfiles. zip The Matlab class scripts and user-written functions are zipped in a file called tsfiles. zip. To get the files, first create an empty directory on your computer. This is where you will store all functions, scripts and data used in the course. Go to D2L, or click on tsfiles. zip to download the zip file to that directory and unzip it there. When you run matlab, be sure that directory is your current matlab working directory. Powerpoint lecture outlines miscellaneous files. Downloadable file other. zip has miscellaneous files used in lectures. Included are Matlab demo scripts, sample data files, user-written functions used by demo scripts, and powerpoint presentations, as pdfs (lect1a. pdf, lect1b. pdf, etc.) used in on-campus lectures. I update other. zip over the semester, and add the presentation for the current lecture within a couple of days after that lecture is given. To run the Matlab scripts for the assignments, you must have your data, the class scripts, and the user-written Matlab functions called by the scripts in a single directory on your computer. The name of this directory is unimportant. Under Windows, it might be something like C:geos585a. The functions and scripts provided for the course should not require any tailoring, but some changes can be made for convenience. For example, scripts and functions will typically prompt you for the name of your input data file and present Spring17 as the default. That is because Ive stored the sample data in Spring17.mat. If you want to avoid having to type over Spring17 with the name of your own data file each time you run the script, edit the matlab script with the Matlab editordebugger to change one line. In the editor, search for the string Spring17 and replace it with the name of your. mat storage file (e. g. Smith2017), then be sure to re-save the edited script.

No comments:

Post a Comment