MATLAB - Laboratorium2

Programy z pętlą i warunkiem. Wykresy XY. Zapis do pliku i odczyt z pliku.

Zadania do wykonania w laboratorium:

Zad1. Pętla WHILE z instrukcją warunkową IF

Cena=145.5;
while 1
  x=input('Zgadnij cene:');
  if x>Cena
     fprintf('Za duzo!\n');
  else
     fprintf('Za malo! \n');
  end
end

Przypomnienie:
Petla typu WHILE:

	while wyrazenie
	  . . . . . instrukcje
	end

B) Dokonaj następujących przeróbek programu:

• Dodaj na pocz. komendy czyszczenia pamięci i ekranu.
• Program ma wczytywać z klawiatury przy pomocy funkcji input dowolną wartość zmiennej Cena (a nie 145.5)
  
• Zmień petlę WHILE tak, aby po wprowadzeniu wartości x równej 0 zakończyło się działanie programu (help while - wyświetli objaśnienia).
  
• Przerobić program aby wykrywał odgadnięcie Ceny (czyli przypadek, że Cena jest równa x ) i pisał jeden z trzech odpowiednich komunikatów.

Zad.2. Pętla FOR i zapis do pliku

A) Napisz program (pr3.m) z pętlą FOR...END, który wyświetla na ekranie liczby nieparzyste od 1 do 9 oraz ich pierwiastki w następującej postaci:

Liczba Pierwiastek 1 1.000 3 1.732 5 2.236 7 2.646 9 3.000

(help elfun - wyświetli objaśnienia funkcji)

B) Zmodyfikuj ten program tak, aby oprócz wyświetlenia tabelki także wyprowadził ją do pliku o nazwie: PIERWIASTKI.TXT.

Możesz wyświetlić zawartość tego pliku poleceniem:
type 'Pierwiastki.txt'

Zad.3. Wykresy dwuwymiarowe - użycie: subplot, ezplot, ezpolar

Funkcje ezplot(...) oraz ezpolar(...) mają kilka postaci. Wyświetl ich opisy przy pomocy help. Wykorzystaj funkcję subplot(lw,lk,nr) aby w jednym oknie graficznym wyświetlić 4 wykresy: ezplot('sin(x)', 'cos(5*x)') ezplot('sin(x)', 'cos(x)') ezpolar('sin(t)',[0 pi]) ezpolar('sin(2*t).*cos(3*t)',[0 pi])

Zad.4. Wykresy dwuwymiarowe - funkcja plot, atrybuty linii oraz opis wykresu

Przy pomocy funkcji plot, dla x z przedziału od 0 do 2*pi narysuj na jednym wykresie trzy przebiegi funkcji: a) sin(x) - czerwoną linią ciągłą b) sin2(x) - niebieską linią kreskowaną c) cos2(x) - czarną linią typu kropka-kreska Dodaj legendę (opis wykresów) Objaśnienia masz w help plot oraz TUTAJ i tu UWAGA: Jeśli x jest ciągiem czyli wektorem czyli macierzą wierszową to x*x albo x^2 są niewykonalne bo to mnożenie macierzowe (wiersz razy kolumna). Poprawne są natomiast tzw. działania tablicowe: x.*x albo x.^2 dotyczące poszczególnych elementów

Zad.5. Interpolacja i aproksymacja

Pomiary co 100 mm ugięcia płaskownika zamocowanego poziomo na jednym końcu i obciążonego pionową siłą na drugim zamieszczono w tabeli:

Sporządź wykres według tej tabeli (linią kropkowaną z zaznaczeniem punktów z tabeli) a następnie zastosuj (zgodnie z opisem w Kompendium):

a = polyfit(x, y, n)
gdzie:	a - szukany wektor współczynników wielomianu,
	x , y - współrzędne danych punktów (np.: dane z pomiarów),
	n - stopień dopasowywanego wielomianu

interpolację "splajnami" przy użyciu funkcji Ys=interp1(x,y,X,'spline') - dodaj do wykresu jako linię czerwoną przerywaną. Zwiększ grubość linii i rozmiar markerów (edycja wykresu) oraz dodaj tytuł i legendę jak poniżej:

Która z tych operacji jest prawidłowa w tym przypadku jeśli wiadomo, że dane były obarczone błędami pomiaru?

Zad.6. Program z pętlą, warunkami, zapisem do pliku i generowaniem wykresu XY.

Dla zginanej belki (jak na rysunku) napisz program generujący wykresy siły tnącej T(x) oraz momentu gnącego Mg(x) w funkcji odległości x od początku belki.
Dane są: a, P. Program ma wczytać dane z klawiatury. Tabela wartości x, T(x), Mg(x) ma być generowana do pliku o nazwie BELKA1.TXT.
Podpowiedzi: użyj instrukcji FOR oraz IF. Wykresy powinny być w osobnych układach współrzędnych.

Zadanie domowe nr 2: