Giao trinh matlab v5.2 P14

Muốn cho chương trình thi hành được thì phải có phương pháp xác định được các vùng trống của bộ nhớ để làm kho chứa dữ liệu. Phương pháp được biết nhiều nhất là thông qua tên của các biến.
Nội dung trích xuất từ tài liệu:
Giao trinh matlab v5.2 P14PhÇn bμi tËp vÝ dô vμ lêi gi¶ielseif u(l+j)==u(l+1)N(l,j)=(U-u(l))*N(l,j-1)/(u(l+j-1)-u(l));elseN(l,j)=(U-u(l))*N(l,j-1)/(u(l+j-1)-u(l)) +(u(l+j)-U)*N(l+1,j-1)/(u(l+j)-u(l+1));endendendendif j==3for l=1:(n+k-3)if u(l+j-1)==u(l)if u(l+j)==u(l+1)N(l,j)=0;elseN(l,j)=(u(l+j)-U)*N(l+1,j-1)/(u(l+j)-u(l+1));endelseif u(l+j)==u(l+1)N(l,j)=(U-u(l))*N(l,j-1)/(u(l+j-1)-u(l));elseN(l,j)=(U-u(l))*N(l,j-1)/(u(l+j-1)-u(l)) +(u(l+j)-U)*N(l+1,j-1)/(u(l+j)-u(l+1));endendendendif j==4X(t)=0;Y(t)=0;Z(t)=0;for l=1:(n+k-4)if u(l+j-1)==u(l)if u(l+j)==u(l+1)N(l,j)=0;elseN(l,j)=(u(l+j)-U)*N(l+1,j-1)/(u(l+j)-u(l+1));endelseif u(l+j)==u(l+1)N(l,j)=(U-u(l))*N(l,j-1)/(u(l+j-1)-u(l));else 133PhÇn 1 - C¬ sëPhÇn bμi tËp vÝ dô vμ lêi gi¶iN(l,j)=(U-u(j))*N(l,j-1)/(u(l+j-1)-u(l)) +(u(l+j)-U)*N(l+1,j-1)/(u(l+j)-u(l+1));endendendendendfor l=1:nX(t)=X(t)+x(l)*N(l,k);Y(t)=Y(t)+y(l)*N(l,k);Z(t)=Z(t)+z(l)*N(l,k);endend %U%hold on;if s1==1subplot(2,1,1);plot3(X,Y,Z);line(x,y,z);hold onsubplot(2,1,2);hold onplot(X,Y,M);line(x,y);hold on endif s1==2subplot(2,1,1);plot3(X,Y,Z);line(x,y,z);hold onsubplot(2,1,2);plot(X,Y,G);line(x,y);hold on endif s1==3subplot(2,1,1);plot3(X,Y,Z);line(x,y,z);hold onsubplot(2,1,2); 134PhÇn 1 - C¬ sëPhÇn bμi tËp vÝ dô vμ lêi gi¶iplot(X,Y,R);line(x,y);hold on endendelseif strcmp(action,Plane);n=input(n=);m=input(m=);k=input(Vao bac k=);h=input(vao bac h=);q=y;while q==yq=input(Neu tiep tuc thi danh y nguoc lai la n =);for i=1:(n+k)if inu(i)=n-k+1;elseu(i)=i-k;endendfor i=1:(m+h)if imw(i)=m-h+1;elsew(i)=i-h;endend%x=input(Nhap vao n.m toa do Px=);%y=input(Nhap vao n.m toa do Py=);%z=input(Nhap vao n.m toa do Pz=);x(1,1)=-3;x(1,2)=-3;x(1,3)=-3;x(1,4)=-3;x(2,1)=-1;x(2,2)=-1;x(2,3)=-1;x(2,4)=-1;x(3,1)=1;x(3,2)=1;x(3,3)=1;x(3,4)=1;x(4,1)=3;x(4,2)=3;x(4,3)=3;x(4,4)=3;y(1,1)=0;y(1,2)=3;y(1,3)=3;y(1,4)=0;y(2,1)=3;y(2,2)=5;y(2,3)=5;y(2,4)=3; 135PhÇn 1 - C¬ sëPhÇn bμi tËp vÝ dô vμ lêi gi¶iy(3,1)=3;y(3,2)=5;y(3,3)=5;y(3,4)=3;y(4,1)=5;y(4,2)=5;y(4,3)=5;y(4,4)=5;z(1,1)=5;z(1,2)=3;z(1,3)=-3;z(1,4)=-5;z(2,1)=5;z(2,2)=3;z(2,3)=-3;z(2,4)=-5;z(3,1)=5;z(3,2)=3;z(3,3)=-3;z(3,4)=-5;z(4,1)=5;z(4,2)=3;z(4,3)=-3;z(4,4)=-5;v=input(vao khoang can ve cua u(1,2..n)=);g=input(vao khoang can ve cua w(1,2..n)=);for i=1:(n+k-1)if u(i)< u(i+1)if u(i)==v-1N(i,1)=1;elseN(i,1)=0;endelseN(i,1)=0;endendfor i=1:(n+k-1)t=N(i,1);endfor i=1:(m+h-1)if w(i)< w(i+1)if w(i)==g-1M(i,1)=1;elseM(i,1)=0;endelseM(i,1)=0;endendfor i=1:(m+h-1)t1=M(i,1);endX1=[];Y1=[];Z1=[];for U=(v-1):0.1:(v-0.1)t=0;for W=(g-1):0.1:(g-0.1)t=t+1; 136PhÇn 1 - C¬ sëPhÇn bμi tËp vÝ dô vμ lêi gi¶ifor i=2:hif i==2for l=1:(m+h-2)if w(l+i-1)==w(l)if w(l+i)==w(l+1)M(l,i)=0;elseM(l,i)=(w(l+i)-W)*M(l+1,i-1)/(w(l+i)-w(l+1));endelseif w(l+i)==w(l+1)M(l,i)=(W-w(l))*M(l,i-1)/(w(l+i-1)-w(l));elseM(l,i)=(W-w(l))*M(l,i-1)/(w(l+i-1)-w(l)) +(w(l+i)-W)*M(l+1,i-1)/(w(l+i)-w(l+1));endendendendif i==3for l=1:(m+h-3)if w(l+i-1)==w(l)if w(l+i)==w(l+1)M(l,i)=0;elseM(l,i)=(w(l+i)-W)*M(l+1,i-1)/(w(l+i)-w(l+1));endelseif w(l+i)==w(l+1)M(l,i)=(W-w(l))*M(l,i-1)/(w(l+i-1)-w(l));elseM(l,i)=(W-w(l))*M(l,i-1)/(w(l+i-1)-w(l)) +(w(l+i)-W)*M(l+1,i-1)/(w(l+i)-w(l+1));endendendendendfor j=2:kif j==2for l=1:(n+k-2) 137PhÇn 1 - C¬ sëPhÇn bμi tËp vÝ dô vμ lêi gi¶iif u(l+j-1)==u(l)if u(l+j)==u(l+1)N(l,j)=0;elseN(l,j)=(u(l+j)-U)*N(l+1,j-1)/(u(l+j)-u(l+1));endelseif u(l+j)==u(l+1)N(l,j)=(U-u(l))*N(l,j-1)/(u(l+j-1)-u(l));elseN(l,j)=(U-u(l))*N(l,j-1)/(u(l+j-1)-u(l)) +(u(l+j)-U)*N(l+1,j-1)/(u(l+j)-u(l+1));endendendendif j==3for l=1:(n+k-3)if u(l+j-1)==u(l)if u(l+j)==u(l+1)N(l,j)=0;elseN(l,j)=(u(l+j)-U)*N(l+1,j-1)/(u(l+j)-u(l+1));endelseif u(l+j)==u(l+1)N(l,j)=(U-u(l))*N(l,j-1)/(u(l+j-1)-u(l));elseN(l,j)=(U-u(l))*N(l,j-1)/(u(l+j-1)-u(l)) +(u(l+j)-U)*N(l+1,j-1)/(u(l+j)-u(l+1));endendendendif j==4for l=1:(n+k-4)if u(l+j-1)==u(l)if u(l+j)==u(l+1)N(l,j)=0;elseN(l,j)=(u(l+j)-U)*N(l+1,j-1)/(u(l+j)-u(l+1));endelse 138PhÇn 1 - C¬ sëPhÇn bμi tËp vÝ dô vμ lêi gi¶iif u(l+j)==u(l+1)N(l,j)=(U-u(l))*N(l,j-1)/(u(l+j-1)-u(l));elseN(l,j)=(U-u(j))*N(l,j-1)/(u(l+j-1)-u(l)) +(u(l+j)-U)*N(l+1,j-1)/(u(l+j)-u(l+1));endendendendend % kt for j=2:kX1(t)=0;Y1(t)=0;Z1(t)=0;for l=1:nfor i=1:mX1(t)=X1(t)+x(l,i).*N(l,k).*M(i,h);Y1(t)=Y1(t)+y(l,i).*N(l,k).*M(i,h);Z1(t)=Z1(t)+z(l,i).*N(l,k).*M(i,h);endendplot3(X1,Y1,Z1);hold onview([3 3 6])endendview([3 3 3])end 139PhÇn 1 - C¬ sëToolBox - Digital Signal ProcessingPhÇn 2 ...