matlab解矩阵函数方程
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/12 20:27:55
I1=Z1\(V1-Va)每个都代表是对应的向量和矩阵你要自己赋值注意是左除号,=旁边那个因为没给具体数据,没法给你列出来你也可以用这个[I1;I2;I3]=[Z11Z12Z13;Z21Z22Z23;
这样的。>>x=solve('a*x^2=1','x')x=1/a^(1/2)-1/a^(1/2)>>x=subs(x,'a',2)x=0.7071-0.7071
如果XA=Bmatlab键入x=B/A或者x=mrdivide(B,AX=-0.33330.33331.33330.66670.33330.33330.66670.83331.3333
n=3;A=reshape(mod(1:n*n,n+1),n,n);B=tril(A.^2-A+1)自己可以设计那种函数再问:3和7是我随便取的,我是想做成这样的:100a110a2a11再答:cle
L=ldl(A)[L,D]=ldl(A)[L,D,P]=ldl(A)[L,D,p]=ldl(A,'vector')[U,D,P]=ldl(A,'upper')[U,D,p]=ldl(A,'upper'
functionI=test()forn=1:4 fork=1:4 I(n,k)=cos((2*(n-1)+1)*(k-1)*pi)/8; &nb
你还是利用高代知识把表达式化成p=.,然后就好算了.
这个问题俺也感兴趣,提供一种方法,仅供参考.clc;clearP=sym('[p11p12p13;p21p22p23;p31p32p33]')A=rand(3,3)At=A'Q=diag(diag(A
你的方程好像有问题,我求解之后得[00-20].没有用及x1+x2+x3+x4=1,如果运用那个条件就会与解发生矛盾.
一元二次方程求出来有两个解,分别放在xx的第一行和第二行:x=solve('1-0.2*x-0.8*x^2-y','x')y=0:0.1:5;xx=eval(x);
随便找本矩阵论的书就可以找的到解线性方程组的方法,有很多.再问:人工解法知道啊,我想知道怎么用matlab实现?再答:见第七,第八页http://wenku.baidu.com/link?url=sO
symskmxw2K=[.]M=[.]X=[.]delta=K-w2*M)*X;s=solve(delta(1),delta(2),delta(3),delta(4),'w2','x');s.w2%w
R=B\A因为行数小于列数,得到的R是基本解
%不用for直接就可以操作A=rand(10)%10*10随即矩阵B=cos(A)%对应元素求余弦C=A^3%矩阵运算D=A.^3%数组运算,对应元素3次取幂E=exp(A)%e的A次方%可以算出来,
X=b乘以A的逆,A的逆可以用inv(A)解出
dsolve('Dx=g*x+b')ans=-b/g+exp(g*t)*C1然后把你的数值带入求解.向量也好,矩阵也好都一样带入,只是乘除号用点运算号.*,./
好吧,再修改成|K-w^2*M|=0实际上就是个三次多项式的求解程序:symsvM=[2100;1410;0141;0012]K=[1-100;-12-10;0-12-1;00-11]f=det(K-
你这个有点像有限元里面的东西.下面的K是我乱写的.你如果给定的不是前三个位移后三个应力的话,你自己得适当地调整,实际上就是调整下面的K1和K2,你一个命令地看,应当可以OK的.K=[2,1,0,0,0
一般来讲指数函数用expm(A)会比funm(A,@exp)快一点,其它函数则无所谓.