Question

matlab数学实验

哪错了？

a=-6:0.5:6;c=a;
[C A]=meshgrid(a,c)
Y=abs((a.^2-a*c*402^(2/3)-2*c.^2*402^(4/3))/(a.^2+c.^2*402^(4/3)+2*a*c*402^(2/3)));
mesh(C,A,Y)

Answer 1

把Y=abs()中的小写a,c都换成拆迹大写的兄悉，羡御乎就可以了。对了，还有矩阵乘*换成点乘.*

追问

Y=abs((a.^2-a*c*402^(2/3)-2*c.^2*402^(4/3))/(a.^2+c.^2*402^(4/3)+2*a*c*402^(2/3)));
这是一个二元方程，我其实是想找到Y<1时的a和c的值 ，应该怎么办啊？？
急急急！！！
解决了这个问题一定给好评！

追答

clear all
clc
a=-6:0.5:6;
c=a;
[A,C]=meshgrid(a,c);
Y=abs((A.^2-A.*C*402^(2/3)-2*C.^2*402^(4/3))/(A.^2+C.^2*402^(4/3)+2*A.*C*402^(2/3)));
subplot(1,2,1);
mesh(A,C,Y)
xlabel('a');ylabel('c');zlabel('Y');
subplot(1,2,2);
plot(A(Y<1),C(Y<1),'*');
xlabel('a');ylabel('c');

用A(Y<1),C(Y<1)来取出Y<1所对应的a,c值

追问

>> syms A C
>> solve('abs((A.^2-A.*C*402^(2/3)-2*C.^2*402^(4/3))/(A.^2+C.^2*402^(4/3)+2*A.*C*402^(2/3)))<1','x','y')

这样做为什么不对？

Answer 2

应该是a.*c两处修改一下，再试试！

追问

Y=abs((a.^2-a*c*402^(2/3)-2*c.^2*402^(4/3))/(a.^2+c.^2*402^(4/3)+2*a*c*402^(2/3)));
这是一个二元方程，我其实是想找到Y<1时的a和c的值 ，应该怎么办啊？？
急急急！！！