木虫 (职业作家)
中国特色 ![]()
|
【答案】应助回帖
cenwanglai: 计算强帖+1, 感谢耐心回复 2012-04-29 20:35:57
不必客气啊。换成下面的代码即可:
clc;
xita=15*pi/180:5*pi/180:85*pi/180;
nxita=length(xita);
x=-30:0.01:30;
nx=length(x);
A=zeros(nxita,nx);
k=61.01*pi/180.;
z=-tan(xita)./tan(k);
for ix=1:nx
A(:,ix)=tan(xita).*(x(ix)-0.5);
end
x1=-0.5.*(z+1)./(1-z);% 求到的根
y1= -tan(k).*(x1+0.5);
Root=zeros(nxita,2);
for ix=1:nxita
Root(ix,1)=x1(ix);
Root(ix,2)=y1(ix);
end
disp(Root);
plot(x,A);
hold on;
y=-tan(k).*(x+0.5);
plot(x,y,'LineWidth',4,'Color','r'); |
|