edongliang
2011-10-23, 18:43
m文件命名为q4.m
function jieguo = q4(x4,z4)
%QQ Summary of this function goes here
% Detailed explanation goes here
syms x1 x2 x3 X4 y1 y2 y3 y4 z1 z2 z3 z4;
x1=0;
y1=0;
z1=0;
x2=sqrt(3)/2;
y2=0.5;
z2=0;
x3=sqrt(3)/2;
y3=-0.5;
z3=0;
y4=0;
W=[1,0.5,0.5;0,sqrt(3)/2,sqrt(3)/6;0,0,sqrt(6)/3];
W1=inv(W);
A=[x2-x1,x3-x1,x4-x1;y2-y1,y3-y1,y4-y1;z2-z1,z3-z1,z4-z1];
T=A*W1;%求亮矩阵的积
T1=inv(T);%求矩阵的逆矩阵
CC=det(T)*det(T1);%求两矩阵行列式相乘
jieguo=CC;
end
matlab 中输入
x4=(0.1:0.1:2);
y4=(0.1:0.1:2);
[X,Y]=meshgrid(x4,y4);
for i=1:20
for j=1:20
QQ(i,j)=q4(X(i,j),Y(i,j));
end
end
我计算出来怎么总是1啊
function jieguo = q4(x4,z4)
%QQ Summary of this function goes here
% Detailed explanation goes here
syms x1 x2 x3 X4 y1 y2 y3 y4 z1 z2 z3 z4;
x1=0;
y1=0;
z1=0;
x2=sqrt(3)/2;
y2=0.5;
z2=0;
x3=sqrt(3)/2;
y3=-0.5;
z3=0;
y4=0;
W=[1,0.5,0.5;0,sqrt(3)/2,sqrt(3)/6;0,0,sqrt(6)/3];
W1=inv(W);
A=[x2-x1,x3-x1,x4-x1;y2-y1,y3-y1,y4-y1;z2-z1,z3-z1,z4-z1];
T=A*W1;%求亮矩阵的积
T1=inv(T);%求矩阵的逆矩阵
CC=det(T)*det(T1);%求两矩阵行列式相乘
jieguo=CC;
end
matlab 中输入
x4=(0.1:0.1:2);
y4=(0.1:0.1:2);
[X,Y]=meshgrid(x4,y4);
for i=1:20
for j=1:20
QQ(i,j)=q4(X(i,j),Y(i,j));
end
end
我计算出来怎么总是1啊