Tài liệu matlap toàn tập_4

Tham khảo tài liệu tài liệu matlap toàn tập_4, công nghệ thông tin, kỹ thuật lập trình phục vụ nhu cầu học tập, nghiên cứu và làm việc hiệu quả
Nội dung trích xuất từ tài liệu:
Tài liệu matlap toàn tập_4 76disp([(1:M) B Ip Pr])format short gCh¹y ch−¬ng tr×nh nµy th× kÕt qu¶ nh− sau:>> Amount=10000Interest Rate=8.9Number of months = 36Payment =317.5321 Amortization Schedule Payment Balance Interest Principle 1.00 9756.63 74.17 243.37 2.00 9511.46 72.36 245.17 3.00 9264.48 70.54 246.99 4.00 9015.65 68.71 248.82 5.00 8764.99 66.87 250.67 6.00 8512.46 65.01 252.53 7.00 8258.07 63.13 254.40 8.00 8001.78 61.25 256.28 9.00 7743.60 59.35 258.19 10.00 7483.49 57.43 260.10 11.00 7221.47 55.50 262.03 12.00 6957.49 53.56 263.97 13.00 6691.56 51.60 265.93 14.00 6423.66 49.63 267.90 15.00 6153.77 47.64 269.89 VÝ dô nµy minh ho¹ cÊu tróc lÆp for vµ if-else-end. Nã còng minh ho¹ viÖc sö dông script M_file.§Ó tÝnh to¸n mét kho¶n cho vay bÊt kú b¹n chØ cÇn thay ®æi d− liÖu vµo ë phÇn ®Çu cña ch−¬ng tr×nhvµ b¹n ch¹y l¹i nã.VÝ dô: Chuçi lªn xuèng VÊn ®Ò: cho x0 lµ mét sè nguyªn bÊt kú. Gi¶ sö chuçi xk ®−îc ®Þnh nghÜa nh− sau: xk+1 = xk/ 2 nÕu xk lµ ch½n vµ xk+1 = 3xk+ 1 nÕu xk lµ lÎChuçi nµy cã thuéc tÝnh g× nÕu chuçi sè dõng l¹i khi xk =1, chuçi ph©n kú hay héi tô vÒ 1. Gi¶i ph¸p: Chóng ta chØ cÇn vßng lÆp while ®Ó xÐt xem khi nµo xk= 1 vµ sö dông cÊu trócif-else-end ®Ó thùc hiÖn viÖc tÝnh to¸n d·y xk. Trong MATLAB th× ch−¬ng tr×nh nh− sau:function up_down% up_down.m script file for up/down sequence proplemx=zeros(500,1); %preallocate storage for x(k)x(1)=round(abs(input(Enter a number> )));k=1;while (x(k)>1)&(k 77 end k=k+1; % increment sequence counterendx=x(x>0) % keep values generated only and dispay themM=0:499;plot(M,x) KÕt qu¶ cña ch−¬ng tr×nh nµy kh¸ thó vÞ, vÝ dô víi x=2m , trong ®ã m lµ mét sè nguyªn th× chuçisÏ rÊt ng¾n (t¹i sao?), h¬n n÷a bÊt cø khi nµo gi¸ trÞ cña mét sè h¹ng trong chuçi lµ luü thõa cña 2 th×chuçi sÏ nhanh chãng dõng l¹i, nh−ng ®èi víi nh÷ng sè x t−¬ng ®èi nhá th× kÕt qu¶ lµ mét chuçi kh¸thó vÞ. VÝ dô x1=27. HÇu nh− tÊt c¶ c¸c gi¸ trÞ ban ®Çu ®Òu sinh ra mét chuçi cã gi¸ trÞ rÊt ngÉu nhiªnnh− h×nh vÏ d−íi ®©y víi x(1)=837799. LiÖu b¹n cã d¸m kÕt luËn chuçi nµy héi tô hay kh«ng!§å thÞ kÕt qu¶ cña ch−¬ng tr×nh víi x(1)=837799 lµ: H×nh 11.1 --------------------oOo-------------------ch−¬ng 12 HµM M_FILE Khi b¹n sö dông c¸c hµm MATLAB nh− inv, abs, angle, vµ sqrt, MATLAB nhËn gi¸ trÞ mµ b¹ntruyÒn vµo, dùa vµo kÕt qu¶ ®ã, tÝnh to¸n kÕt qu¶ cña hµm vµ tr¶ l¹i cho b¹n kÕt qu¶ tÝnh to¸n. C¸clÖnh tÝnh to¸n b»ng hµm còng nh− c¸c biÕn trung gian ®−îc t¹o ra bëi c¸c lÖnh nµy b¹n ®Òu kh«ngnh×n thÊy, tÊt c¶ nh÷ng g× b¹n tr«ng thÊy chØ lµ c¸c gi¸ trÞ nhËp vµo vµ c¸c gi¸ trÞ ®−a ra, v× vËy cã thÓcoi mét hµm nh− mét c¸i hép ®en. C¸c thuéc tÝnh nµy lµm cho hµm trë lªn rÊt h÷u dông ®èi víi c¸clÖnh tÝnh to¸n mµ ph¶i dïng ®Õn c¸c hµm to¸n häc phøc t¹p th−êng xuÊt hiÖn khi b¹n gi¶i quyÕt 78nh÷ng vÊn ®Ò lín. Dùa vµo −u ®iÓm nµy, MATLAB cung cÊp mét cÊu tróc ®Ó b¹n cã thÓ tù t¹o méthµm cho m×nh d−íi d¹ng mét M_file. Hµm flipup d−íi ®©y lµ mét vÝ dô vÒ viÖc dïng hµm M_file:function y=flipup(x) % FLIPUP Flip matrix in up/down directiopn. % FLIPUP(x) return x with columns preserved and rows flipped % in the up/down direction. For example. % % x=1 4 becomes 36 % 25 25 % 36 14 % % See also FLIPLR, ROT90, FLIPDIM. % Copyright (c) 1984-96 by the MathWork, Inc. % $Revision: 5.3 $ $Date: 1996/10/24 18: 41: 14 $if ndim(x)~=2 error( ‘X must be a 2-D matrix.’);end[m, n] = size(x);y = x(m: -1: 1, :); Mét hµm M_file cã vÎ rÊt gièng ...