数学建模社区-数学中国

标题: 求助matlab高手 [打印本页]

---

作者: 1134864661    时间: 2014-12-31 16:32
标题: 求助matlab高手
任给一矩阵A，若A为4×m（m≥4）,利用初等行变化为最简阶梯阵（不能用rref函数），并给出矩阵的秩。否则A为方阵，给出矩阵的转置和矩阵的秩。如果不是方阵用rref函数给出矩阵的行最简阶梯阵及矩阵的秩。怎么用matlab编写

---

作者: QYY1994    时间: 2015-1-12 20:30

clc
clear all
M1=[1 2 3 4 5 6 7 8 9 10
   3 4 5 6 7 8 9 10 11 12
   5 7 7 8 9 10 11 12 13 14
   7 8 9 10 11 12 13 14 15 16];
M2=[1 2 3 4
    3 4 5 6
    5 6 7 8
    7 8 9 10];
M3=[1 2 3
    2 3 4
    3 4 5
    4 5 6];
A=M1;%A=M2;%A=M3;
if size(A,2)>4
for i=2:4
A(i,:)=A(i,:)-A(i,1)/A(1,1)*A(1,:);
end
k=find(A(2:4,2));
k=k(1)+1;
t=A(2,:);
A(2,:)=A(k,:);
A(k,:)=t;
for i=3:4
A(i,:)=A(i,:)-A(i,2)/A(2,2)*A(2,:);
end
k=find(A(3:4,3));
c=size(k);
if (c(1)~=0)
k=k(1)+2;
t=A(3,:);
A(3,:)=A(k,:);
A(k,:)=t;
A(4,:)=A(4,:)-A(4,3)/A(3,3)*A(3,:)
else
k=find(A(3:4,4));
k=k(1)+2;
t=A(3,:);
A(3,:)=A(k,:);
A(k,:)=t;
A(4,:)=A(4,:)-A(4,4)/A(3,4)*A(3,:)
end
r1=rank(A)
else if(size(A,2)==4)
        T=A'
        r2=rank(A)
    else B=rref(A)
        r3=rank(A)
    end
end

这几天在准备期末考试，没有及时回复，抱歉。。。

---

作者: 1134864661    时间: 2015-1-15 14:34
谢谢。再问一个题，1.作曲线y=x(1-x)在[0,1]上的转动切线从几何上说明水平切线的存在性。2.在【0,1】上取100个点，生成数据点z=sin(y)，用这些进行最近点插值，求出x=0.9999处的值同时对数据点进行三次多项式拟合怎么编写？

---

作者: 1134864661    时间: 2015-1-15 14:35
1.作曲线y=x(1-x)在[0,1]上的转动切线从几何上说明水平切线的存在性。2.在【0,1】上取100个点，生成数据点z=sin(y)，用这些进行最近点插值，求出x=0.9999处的值同时对数据点进行三次多项式拟合怎么编写？

---

欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)

Powered by Discuz! X2.5