正在学习中，希望跟论坛一起成长！

给我一个比比吧！！！

偶是刚刚上来的，给点希望吧！！！

向专家的目标进发

新人前来报道，希望不断努力不段进取？

申请做基础问题版块版主。帮助他人同时也可提高自己，共同学习，共同进步。请管理员核查。

欢迎你加入我们：

请让我们用一周时间来评估一下，一周以内给你答复。

谢谢你的支持。。。

你好，也许你觉得我是新注册，没这个资格，我是matlab的爱好者，呵呵，我想申请虚拟现实的版主，我的研究方向是智能机器人，我要研究的是遥操作平台系统，所以我对虚拟现实方面接触算比较多吧，所以我想跟论坛一起成长，刚好版主的位置缺了，可以让我先坐一下不？如果以后有更合适的人选我愿意让位的。呵呵，支持论坛。

顶个！
现在没有满足条件哦

以后有的是机会，

继续努力中、、、、、、、、、、、

程序及运行
function function_2_1
n0=input('enter n0=');
n1=input('enter n1=');%纤维的折光系数
r=input('enter r=');%异形截面纤维轮廓中心基圆的半径
r1=input('enter r1=');%异形截面纤维轮廓顶端基圆的半径
R=input('enter R=');%异形截面纤维轮廓外接基圆的半径
I0=input('enter I0=');%入射光的光强
n=input('enter n=');%theta0的等分个数
m=input('enter m=');%thetat的等分个数
t=asin((r-r1)/(R-r1));%AB直线部分与Y轴正向的夹角
thetat=(pi/2)-t;%AB直线部分光线的法线与Y轴正向的夹角
lAB=sqrt((R-r1).^2-(r-r1).^2);%AB直线部分长度
theta0=0;deltatheta0=((pi/2)-0)/m;%入射光的初始角赋值和theta0的等分角度
while (theta0<=(pi/2))
    k=0:m;
    theta0=theta0+k*(deltatheta0);
    thetaR2=thetat-theta0;
    feR2=asin((n0*sin(thetaR2))/n1);
    Rr2=(1/2)*(((tan(thetaR2-feR2).^2)/(tan(thetaR2+feR2).^2))+((sin(thetaR2-feR2).^2)/(sin(thetaR2+feR2).^2)));
    Ir2=(I0*Rr2*sin(thetaR2))/lAB;
    k=k+1;
end
x1=0:0.01pi/2);
y1=0:0.01(pi/2)-t);
[x,y]=meshgrid(x1,y1);
z=Ir2;
surf(x,y,z)
colormap,colorbar
>> function_2_1()
enter n0=1
enter n1=1.5
enter r=2
enter r1=1
enter R=5
enter I0=1
enter n=0
enter m=157
Warning: Matrix dimensions must agree, not rendering mesh.
多谢了