matlab解方程小数,MatLab入门及解方程_梦里香巴拉的空间

% MATLAB 6.1 实用指南

help elmat

help eye

root=solve('a*x^2+b*x+c=0','x')

% 续行 ...

% LOOKFOR Search all M-files for keyword.

lookfor integral

% pi ，inf无穷大 i(j)复数中的虚数单位

MATLAB的输入与输出语句--------------------------------------------------------------------------------

输入语句

输入数值

?x=input('please input a number:')

please input a number:22

x = 22

输入字符串

?x=input('please input a string:','s')

please input a string:this is a string

x = this is a string

输出语句

输出显示命令

自由格式 (disp)

?disp(23+454-29*4)

361

?disp([11 22 33; 44 55 66; 77 88 99])

11 22 33

44 55 66

77 88 99

?disp('this is a string')

this is a string

格式化输出 (fprintf)。

fprintf('The area is %8.5f\n', area) % 注意输出格式前须有%符号，

%跳行符号须有\符号

The area is 12.56637 % 输出值为8位数含5位小数

在这里你如果学过c语言就能很好的理解了。

错误消息显示命令

?error('this is an error')

??? this is an error

%%%%%%%%%%%%%%%     数组           %%%%%%%%%%%%%%%%%%%%%%%%% 创建数组x＝(0:0.5:2) first:increment:last 方法2： linspace(first_value,last_value,number)

%     logspace(first_value,last_value,number)

A=[2,5,6,3;1,2,3,4]

B=[-7,9,-2,0;1,1,1,1]

c=A+B

c=c+2

c=c/2

%%    随机数[10,50] to generate a 5-by-5 array of uniformly distributed random numbers on the interval [10,50]

a = 10; b = 50;

x = a + (b-a) * rand(5)

% 求行列式大小 det

a=[1 2 0;2 5 -1; 4 10 -1]

b=det(a)

%求逆及转置

c=inv(a)

b=det(c)

B=C'

%%

A(:,j)is the j-th column of A

A(i,:)is the i-th row of A

A(:,:)is the equivalent two-dimensional array.

For matrices this is the same as A.

A(j:k)is A(j), A(j+1),...,A(k)

A(:,j:k)is A(:,j), A(:,j+1),...,A(:,k)

A(:,:,k)is the kth page of three-dimensional array A.A(i,j,k,:)is a vector in four-dimensional array A. The vector includes A(i,j,k,1), A(i,j,k,2), A(i,j,k,3), and so on.A(:)is all the elements of A, regarded as a single column. On the left side of an assignment statement, A(:) fills A, preserving its shape from before. In this case, the right side must contain the same number of elements as A.

%%%

A = pascal(3);

B = magic(3);

m = 3; n = 3;

for i = 1:m

for j = 1:n

C(i,j) = A(i,:)*B(:,j);

end

end

%矩阵幂运算 ^    不同于数组的幂运算 .^

A=[1,1,1;1,1,1;1,1,1]

B=A.^2

C=A^2

%单位阵eye(n,m)

a=eye(4,3)

%元素全1矩阵

a=ones(4,3)

%零矩阵 zeros(n,m)

a=zeros(4,3)

%魔方矩阵 magic(n) 随机阵rand(n,m)

magic(4)

%伴随矩阵 compan(A) 注意A为向量，不能是矩阵

u=[1 0 -7 6];

A=compan(u)

%%(x+y)n杨辉三角形表组成的矩阵称为帕斯卡(Pascal)矩阵求(x+y)5系数,对角线值

pascal(6)

% 函数rot90(A,k)将矩阵A旋转90?的k倍，当k为1时可省略

% 求方阵A的逆矩阵可调用函数inv(A)

% 特征值和特征向量的函数是eig(A)

% 有AX=B 则X=?

A=[21,34,20;5 78 20;21 14 17; 34 31 38];

B=[10 20 30 40]';

X=B\A

% diag(A)函数用于提取矩阵A主对角线元素，产生一个具有min(m,n)个元素的列向量。

%%%%%%%%%%%%%%     % 多项式运算           %%%%%%%%%%%%%%%%%%%%%%%%%%

T=[2 5 0 4 1 4]

poly2sym(T)

%f=2x^5+5x^4+4x^2+x+4

T=[2 5 0 4 1 4]

r=roots(T)

poly(r)

x=-2.7709

f=2x^5+5x^4+4x^2+x+4

%%

%多项式乘积 C=Conv(a,b) eg. (x+1)(x2+x+3)

a=[0,1,1]

b=[1,1,3]

mul=conv(a,b)

% 多项式除 (s4+7s3+16s2+18s+8)/(s+3)

a=[1,3]

b=[1,7,16,18,8]

[q,r]=deconv(b,a)

%求多项式值：polyval与polyvalm，它们的输入参数均为多项式系数向量P和自变量x。

% 两者的区别在于前者是代数多项式求值，而后者是矩阵多项式求值。

%多项式＋ 乘conv(T1,T2);除 [T r]=deconv(t1,t2)

T1=[2 5 0 4 1 4];

T2=[0 0 5 1 3 2];

T=conv(T1,T2);

% 求导 polyder

T1=[2 5 0 4 1 4];

h=polyder(T1);

%多项式求值 求x=[3 4]处的值 h=polyval(T1,[3 4])

T1=[2 5 0 4 1 4];

s=[1,4,5]

h=polyval(T1,s)

%求根 h=roots(T1)

h=roots(T1)

% 求多项式x4+8x3-10的根。命令如下：

A=[1,8,0,0,-10];

X=roots(A)

G=poly(X)

% 若已知多项式的全部根，则可以用poly函数建立起该多项式，其调用格式为：P=poly(x)

% 若x为具有n个元素的向量，则poly(x)建立以x为其根的多项式，且将该多项式的系数赋给向量P。

%多项式的拟合用 [p,s]=polyfit(x,y,n) 最小2乘法

x=linspace(0,2*pi,100);

y=cos(x)

T=polyfit(x,y,6)

y1=polyval(T,x)

plot(x,y,'go',x,y1,'b--')

%%%    关系与逻辑操作符    测试函数    %%%%%%%%

%     Eq(A,B) Lt(A,B) Gt(A,B) Le(A,B) Ge(A,B)

%     And(A,B) Or(A,B) Not(A) Xor(A,B) Any(A) All(A)

A=[1 0 2 1]

B=[0 0 1 1]

C=And(A,B)

%%%%% 通用数据分析函数 %%%%%%%

%    Max Min Mean Median Std Sum Prod求积 hist trapz Sort Sumsum CumProd

a=[1 2 5; 4 5 3]

[y,I]=max(a)

[y,I]=max(a,[],2)

y=[0 1 2]

trapz(y)

% diff 差分 Gradient(F) 梯度

% 对多项式求导数的函数是：

p=polyder(P) %：求多项式P的导函数

p=polyder(P,Q)    %：求P·Q的导函数

[p,q]=polyder(P,Q) %：求P/Q的导函数，导函数的分子存入p，分母存入q。

%%%%%%%%%%%%%%%%%%%%%%         符号运算         %%%%%%%

% (1)利用单引号来生成符号表达式。

% (2)用sym函数建立符号表达式。

% (3) 使用已经定义的符号变量组成符号表达式

% syms    符号变量名1 符号变量名2 …

f='a*x^2+bx+c=0'

f=sym('ax+b=0') % sym生成符号表达式

symvar('a*x+y')

diff('x^n') %求导

diff('x^n','n') %对n求导

% 符号基本代数运算 symadd,symsub,symmul,symdiv,sympow

f='4*x+5';

g='2*x^2+5*x+6';

symadd(f,g)

eval(t) %其中t为字符串。它的作用是把字符串的内容作为对应的MATLAB语句来执行。