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