aqyugu
2012-03-18, 11:31
我在编写光反射增益时，遇到一些问题，求解惑，谢谢！简化表述方式，就是
.......................................................................
P=∫∫p（a，b，x，y）dadc, a∈[0,5]c∈[0,3],x∈[0,5]y∈[0,5].画出P在（x，y）上
的图。
.......................................................................
我编了个程序，大致是
-------------------------------------------
for x=0:0.5:5;
for y=0:0.5:5;
P=dblquad(@(a,c)p(a,c,x,y),0,5,0,3);
end
end
mesh(x,y,P)
-------------------------------------------
我遇到的问题就是
1：p（a,c,x,y）的表达式太长了，但有@（a，c），我不清楚怎样把它变短（比如p=cos（
a）*cos（c）*sin（x）*sin（y），记s=cos（a）*cos（c），t=sin（x）*sin（y），即
p=s*t，表达式变短了）；
2：这里画图时用了mesh不知道可不可以？plot3可以吗？
3：我这样编写的程序是不是有整体错误，直接实现不了题目的目的。有别的函数可以用吗
？或者更好的编写方法？
求解惑，谢谢啦！
.......................................................................
P=∫∫p（a，b，x，y）dadc, a∈[0,5]c∈[0,3],x∈[0,5]y∈[0,5].画出P在（x，y）上
的图。
.......................................................................
我编了个程序，大致是
-------------------------------------------
for x=0:0.5:5;
for y=0:0.5:5;
P=dblquad(@(a,c)p(a,c,x,y),0,5,0,3);
end
end
mesh(x,y,P)
-------------------------------------------
我遇到的问题就是
1：p（a,c,x,y）的表达式太长了，但有@（a，c），我不清楚怎样把它变短（比如p=cos（
a）*cos（c）*sin（x）*sin（y），记s=cos（a）*cos（c），t=sin（x）*sin（y），即
p=s*t，表达式变短了）；
2：这里画图时用了mesh不知道可不可以？plot3可以吗？
3：我这样编写的程序是不是有整体错误，直接实现不了题目的目的。有别的函数可以用吗
？或者更好的编写方法？
求解惑，谢谢啦！