急急急!用matlab解矩阵方程组 -1.05t(1)+0.81t(2)-1.032=0 t(i-1)+t(i+1)-2t(i)=0 0.81t(10)-0.91t(11)
w=[2000,3000,1000,500,4000,2000,1000];
t=[48.7,51.9,61.2,71.8,48.7,51.9,64.0];
f=-[t,t]'; %求最大值,化成求负数的最小值
A=[
w,w;
t,t;
0,0,0,0,t(5),t(6),t(7),0,0,0,0,t(5),t(6),t(7); %t5*n51+t5*n52+t6*n61+t6*n62+t7*n71+t7*n72<=302.7
1,0,0,0,0,0,0,1,0,0,0,0,0,0; %n11+n12<=8
0,1,0,0,0,0,0,0,1,0,0,0,0,0; %n21+n22<=7
0,0,1,0,0,0,0,0,0,1,0,0,0,0; %n31+n32<=9
0,0,0,1,0,0,0,0,0,0,1,0,0,0; %n41+n42<=6
0,0,0,0,1,0,0,0,0,0,0,1,0,0; %n51+n52<=6
0,0,0,0,0,1,0,0,0,0,0,0,1,0; %n61+n62<=4
0,0,0,0,0,0,1,0,0,0,0,0,0,1; %n71+n72<=8
];
B=[80000;2040;302.7;8;7;9;6;6;4;8];
intlist=ones(14,1); %n全部都是整数
xmin=zeros(14,1); %n变量最小值都是0
xmax=inf*ones(14,1); %最大值不设置
ctype=-1*ones(10,1); %所有约束都是小于等于约束,这里有10个条件
[n,how]=ipslv_mex(f,A,B,intlist,xmax,xmin,ctype)
%如果how=0的话,说明找到最优解
fmax=-f'*n %这里还原负值
%%%%%%%%%%%%
ipslv_mex得去网上找
分析:pQ=0说明Q是不满秩的(否则p只有全0解与下一个条件矛盾),p一定存在于Q的核空间中。p*ones(n)=1说明p是归一化的,这个最后处理一下就行。
首先用rank(Q)<n确认一下Q不满秩,再用null(Q)命令求出核空间的基,从中任取一个归一化即可。
程序:
if rank(Q)<n
ps = null(Q); % 或者 ps=null(Q,'r');更精确些
p = ps(:,1);
p = p/sum(p);
end
clc
clear all
A=zeros(11);
B=zeros(11,1);
A(1,1)=-1.05;
A(1,2)=0.81;
A(11,10)=0.81;
A(11,11)=-0.91;
B(1)=1.032;
for ii=2:10
A(ii,ii-1)=1;
A(ii,ii+1)=1;
A(ii,ii)=-2;
end
t=A\B
最后一个方程没给全吧,我就当最后是个零来给程序。
clc
clear all
A=zeros(11);
B=,1);
A(1,1)=-1.05;
A(1,2)=0.81;
A(11,10)=0.81;
A(11,11)=-0.91;
B(1)=1.032;
for ii=2:10
A(ii,ii-1)=1;
A(ii,ii+1)=1;
A(ii,ii)=-2;
end
t=A\B 就行了鐧惧害鍦板浘
用matlab做下题急急急!
结果:哈桑不用担心,黄金还有盈余的 代码:clear;半椭球相关系数值 a=30.6;b=29.6;c=30;xmin=0;xmax=pi\/2;ymin=0;ymax=pi\/2;被积函数,在极坐标系下进行积分 integellipse=@(x,y)4*a*b*cos(x).*((sin(x).^2)+...c^2\/a^2*(cos(x)).^2.*(cos(y)).^2+...c^2\/...
加急!用matlab画一个半径为5的球
x0=2;y0=3;z0=4;%球心 r=5;%半径 下面开始画 [x,y,z]=sphere;mesh(r*x+x0,r*y+y0,r*z+z0)axis equal === 或者:clear;clc;[x,y]=meshgrid(-1:0.01:1);z=sqrt(1-x.^2-y.^2);z1=real(z);mesh(x,y,z1);hold on;z=-sqrt(1-x.^2-y.^2);z2=real(z);m...
Matlab问题 急!!!求!
题主的问题可以这样来解决。根据题意,输入x为矩阵,而输出y也为矩阵。对于任意矩阵进行计算,软件系统要求用点运算符进行运算,所以我们自定义函数应这样来写,即 function y=fx(x)y=1.\/((x-2).^2+0.1)+1.\/((x-3).^4+0.01);end 然后,在命令窗口中或主程序中,输入 x=rand(5,3)...
matlab 急急急急急急!!!
可以在matlab中建立一个m文件,将以下function内容复制到m文件即可。使用该m文件,只需在matlab的命令行中敲入[n_couples,i_month]=population(1,24)即可。 function内容如下:---function内容开始--- function [n_couples,i_month]=population(mature_couples,months)该程序是用来计算每个月的兔子对数 ...
急求!!! 大学数学,用matlab解决问题,题目是一维抛物型偏微分方程差分解...
在MATLAB中编写脚本文件:A = 0.5; %方程系数 it0 = inline('sin(pi*x)','x'); %初始条件 bx0 = inline('0'); bxf = inline('0'); %边界条件 xf = 2; M = 80; T = 0.1; N = 100;[u1,x,t] = EF_Euler(A,xf,T,it0,bx0,bxf,M,N);figure(1),clf,mesh(u1)x...
怎样用matlab读取一个文件夹下的多个子文件夹中的多个jpg图片,急用...
这样:imgPath = 'E:\/imageData\/'; % 图像库路径 imgDir = dir([imgPath '*.jpg']); % 遍历所有jpg格式文件 for i = 1:length(imgDir) % 遍历结构体就可以一一处理图片了 img = imread([imgPath imgDir(i).name]); %读取每张图片 end ...
急急急!!!高分求用MATlAB解答 在线等
function S = int2sum(X)INT2SUM 把一个正整数分解为若干连续正整数之和。S = INT2SUM(X)输入: X 为正整数,小数则去掉小数部分,负数取绝对值 返回: S 为列数为2的矩阵,每行表示一种分解的情况,分解出的连续整数以第一列为首项,第二列为尾项。如果无法分解成连续正整数之和,...
用matlab求解多元一次方程组 求大神指导 急!!急!!急!! 救命啊!!
3、就题主现在的问题而言,就是无解的,这可以用线性规划的方法验证:Ae=[2 3 5 6; 1 1 1 1];be=[89; 20];lb=[5 0 0 0];ub=[10 5 5 5];x=linprog([1 0 0 0],[],[],Ae,be,lb,ub)得到结果是无法找到可行解。上述结论也可以通过Lingo进行验证:max=x1;2*x1+3*x2+5*...
急!如何用MATLAB构建理想低通滤波器,巴特沃斯低通滤波器,指数低通滤 ...
说明:buttord函数使用阻带指标计算3dB截止频率,这样阻带会刚好满足要求,而通带会有富余。2.buttap(N)[z0,p0,k0]=buttap(N)用于计算N阶巴特沃斯归一化(3dB截止频率Ωc=1)模拟低通原型滤波器系统函数的零、极点和增益因子。说明:如果要从零、极点模型得到系统函数的分子、分母多项式系数向量ba、...
急急!!MATLAB高手来
可能我没明白你的意思,不过用:ezplot('5*0.4476*x.^4*50*pi\/9-4*2.1740*x.^3*50*pi\/9+3*3.4935*50*pi\/9*x.^2 -2*1.4556*50*pi\/9*x.^1+0.2330*50*pi\/9')ezplot('5*0.4476*x.^4*50*pi\/9-4*2.1740*x.^3*50*pi\/9+3*3.4935*50*pi\/9*x.^2 -2*1.4556...
路佩西乐: 先将xP=0两边取转置,得P'x'=0,求出x'再转置回来求出x. 当矩阵方程P'x'=0,rank(P)=r<n时,即P是奇异的,这时方程有无穷多个解,用MATLAB可求出它的一个基本解系, 基本解系的线性组合就是它的通解(一般解) 而求基本解系用matlab 中的命令 x=null(P', r )即可.其中:r=rank(P).
新密市19624657714: 急求matlab编程求矩阵方程解 - ?
路佩西乐: %楼主说的不可能,下面是我修改了一下,A 是4行6列1 2 3 4A=[4,3,2,1;3,4,1,2;1,2,4,3;6,5,6,9;2,3,4,5;12,2,1,6]'; b=[1,2,3,4]'; x=b\A; disp(x) 计算结果: 1 2 3 4 5>> A=[4,3,2,1;3,4,1,2;1,2,4,3;6,5,6,9;2,3,4,5;12,2,1,6]'; b=[1,2,3,4]'; x=b\A; disp(x)0.66670.73330.96672.33331.33331.4333
新密市19624657714: matlab解矩阵方程求教啊 - ?
路佩西乐: 直接弄个高斯列主元 消去法就是了e68a84e8a2ade799bee5baa6e79fa5e9819331333335316435 function X=Gauss(A,B) n=length(B); X=zeros(n,1); c=zeros(1,n); d1=0 for i=1:n-1max=abs(A(i,i));m=i;for j=i+1:nif max<abs(A(j,i))max=abs(A(j,i...
新密市19624657714: 用matlab求解矩阵方程AX=B - 2X - ?
路佩西乐: 程序:展开全部 A=[1 2 3;1 3 1; 0 1 2];B=[1 0; 0 1;-1 0];X=(A+2)\B%化简矩阵方程.AX=B对应X=A\B.XA=B对应X=B/A结果:X =6.0000 -0.2500-3.0000 0.5000-1.0000 -0.2500
新密市19624657714: matlab求矩阵方程组 - ?
路佩西乐: 重新整理一下你的方程.你的第一个方程是 (pi0,pi1)=(pi0,pi1)*[A00,A01;A10,R*A2] 取个转置,就是 (pi0,pi1)(T) = [A00,A10;A01,R*A2]*(pi0,pi1)(T) 相当于求解线性方程组 Ax = 0 A = [A00 - 1,A10;A01,R*A2 - 1] x = (pi0,pi1)(T) 你还有...
新密市19624657714: matlab如何解矩阵方程AX=b,A为3*3矩阵,b为3*1矩阵,X为3*1矩阵, - ?
路佩西乐:[答案] X=b乘以A的逆,A的逆可以用inv(A)解出
新密市19624657714: matlab求解矩阵方程组矩阵A(8行9列)矩阵X(8行1列)矩阵B(8行1列)每个值都为0AX=B,就是AX=0,矩阵A已知,求X用matlab怎么解 - ?
路佩西乐:[答案] 先将xP=0两边取转置,得P'x'=0,求出x'再转置回来求出x. 当矩阵方程P'x'=0,rank(P)=r
新密市19624657714: Matlab如何求解矩阵微分方程组的符号解例如:dX/dt=HX,X为列向量,H为矩阵H=[2*c,A,0;0,g,0;A,0,a],c,a,g,A为符号常量,初值X(0)=[0,1,0]',求方程组的解X? - ?
路佩西乐:[答案] 使用dsolve函数 dsolve('Dx=2*c*x+A*y','Dy=g*y','Dz=A*x+A*z','x(0)=0','y(0)=1','z(0)=0') 就是把方程组写出三个微分方程进行计算 祝你学习愉快!
新密市19624657714: matlab 矩阵方程求解 - ?
路佩西乐: 好吧,再修改成 |K-w^2*M|=0 实际上就是个三次多项式的求解 程序:syms v M=[2 1 0 0;1 4 1 0;0 1 4 1;0 0 1 2] K=[1 -1 0 0;-1 2 -1 0;0 -1 2 -1;0 0 -1 1] f=det(K-v*M) solve(f) 求出了v,w就是v的开方,你应该会求吧
新密市19624657714: matlab求矩阵方程 - ?
路佩西乐: It may not be the best answer.a=[0 1;1 -1]; q=[1 0;0 1]; syms p1 p2 p3 p4 p; p=[p1 p2;p3 p4]; f=a'*p+p.'*a+q; pp=solve(f(1),f(2),f(3),f(4),'p1','p2','p3','p4'); [pp.p1 pp.p2;pp.p3 pp.p4] ans = [ -1/2-p4, -1/2+p4] [ -1/2, p4]