21.510.50-0.5-1-1.5-2-2x-1.5-1-0.50Y0.511.52

%以上方法可以绘出正常的横坐标为y纵坐标为x的图像，但发现在y=0处x延伸至正负无穷。

h1=ezplot(y1,[-2 2 -2 2],1); h2=ezplot(y2,[-2 2 -2 2],1); -(2 - 3 x2)1/221.510.50-0.5-1-1.5-2-2x-1.5-1-0.50x0.511.52 reset(symengine); syms x y S;

%以上方法绘出的图像存在一个空隙，且默认为y-x图像。

S = dsolve('Dy*y/5+x/4=0','x');

ezplot(subs(y^2-(S(1))^2, 'C3', 1),[-2,2 -2,2],2); grid on;

%用椭圆方程绘图不产生间隙

(5 x2)/4 + y2 - 221.510.50-0.5-1-1.5-2-2y-1.5-1-0.50x0.511.52

24 一阶微分方程

syms a b;

ys=dsolve('Dy-a*x^2-b*x=0','y(0)=2',x)

ys =

(x^2*(3*b + 2*a*x))/6 + 2

25 边值问题

fs=dsolve('Df-3*f=4*g,Dg+4*f=3*g','f(0)=0,g(0)=1')

fs =

g: [1x1 sym] f: [1x1 sym]

fs.g fs.f

ans =

cos(4*t)*exp(3*t) ans =

sin(4*t)*exp(3*t)

第三章

3.行下标列下标

rng('default'); A=rand(3,5); L=A>0.5

L =

1 1 0 1 1 1 1 1 0 0 0 0 1 1 1 [a,b]=find(L==1)

IND=sub2ind(size(A),a,b)

IND = 1 2 4 5 8 9 10 12 13 15

4.循环运算、数组运算

t=0:0.1:10; N=length(t);

y1=zeros(size(t)); for k=1:N

y1(k)=1-exp(-0.5*t(k))*cos(2*t(k)); end

plot(t,y1); xlabel('t'); ylabel('y1');

1.51y10.50012345t678910 y2=1-exp(-0.5*t).*cos(2*t); plot(t,y2); xlabel('t'); ylabel('y2');

1.51y20.50012345t678910 5.回答问题

clear all;

A=magic(3);B=rand(3); A*B B*A

ans =

5.4072 11.5771 3.0037 6.3884 10.3215 4.9680 2.7058 7.5337 4.8496 ans =

2.5916 3.8303 5.2097 3.4833 5.6313 3.6800 10.9646 9.0086 12.3554

相同，对于矩阵而言对位相乘无差异

不相同，点乘与矩阵乘法进行的不是同一种运算。 不相同，左乘右乘运算不同。

相同，A左点除B等同于B右点除A，均是对位计算。 不相同，左除右除运算亦不相同。

A*A\\B-B A*(A\\B)-B

A*(A*inv(B))-B