Question

求大神解决一个matlab问题

百度不给过只能把程序截图了，如果需要，我可以把这段程序发过去，麻烦各位了。
t=20:0.1:180;r=10;s=50;a=480;b=480;o=71/1852;d=60*180/pi*asin(sin(pi*s/(60*180))/sin(pi*t/180));
v=sqrt(a^2+b^2-2*a*b*cos(t));
l=1/2.996*v*r*sqrt(2)/a/sin(t);
e=v/a/sin(t);
h=l/500.*(exp(e*o./l)-exp(-e*o./l)).*(exp(-2.*d/l)-exp(-(2.*d+250)./l))*2*0.39*19/1652;semilogy(t,h)
??? Error using ==> mrdivide
Matrix dimensions must agree

Answer 1

clear all;clc;
t=20:0.1:180;
r=10;
s=50;
a=480;
b=480;
o=71/1852;
d=60*180/pi*asin(sin(pi*s/(60*180))./sin(pi*t/180));
v=sqrt(a^2+b^2-2*a*b*cos(t));
l=1/2.996*v*r*sqrt(2)/a./sin(t);
e=v/a./sin(t);
h=l/500.*(exp(e*o./l)-exp(-e*o./l)).*(exp(-2.*d/l)-exp(-(2.*d+250)./l))*2*0.39*19/1652;
semilogy(t,abs(h));

追问

这个做出来的图和原文章的那个图不一样。

追答

那你要看看你的算法有没有问题
另外，这句改一下：
h=l/500.*(exp(e*o./l)-exp(-e*o./l)).*(exp(-2*d./l)-exp(-(2*d+250)./l))*2*0.39*19/1652;

追问

你再看这个，其实是要做三个图，都是一个算式呀，这个就是以r为变量，然后做出来的图和书上是一样的
r=4:0.1:16;t=45;s=50;a=480;b=480;o=71/1852;d=60*180/pi*asin(sin(pi*s/(60*180))/sin(pi*t/180));
v=sqrt(a^2+b^2-2*a*b*cos(t));
l=1/2.996*v*r*sqrt(2)/a/sin(t);
e=v/a/sin(t);
h=l/500.*(exp(e*o./l)-exp(-e*o./l)).*(exp(-2*d./l)-exp(-(2*d+250)./l));semilogy(r,h)

追答

clear all;clc;

r=4:0.1:16;

t=45;s=50;

a=480;b=480;

o=71/1852;

d=60*180/pi*asin(sin(pi*s/(60*180))./sin(pi*t/180));

v=sqrt(a^2+b^2-2*a*b*cos(t));

l=1/2.996*v*r*sqrt(2)/a./sin(t);

e=v/a./sin(t);

h=l/500.*(exp(e*o./l)-exp(-e*o./l)).*(exp(-2*d./l)-exp(-(2*d+250)./l));

semilogy(r,abs(h));

追问

这个图我是做出来了，我的意思是换了个变量，变量t变成数组后，就做不出来了，我是问第一个程序，第二个我就是拿来说可能算数的公式是对的

追答

r=20:0.1:180;
t=10;
--------------你这两句搞反了
semilogy(r,abs(h));
---------这句也要跟着变

追问

不是不是这两个是两个图，不是搞反了，你可能没有懂我的意思，这个里面t代表的是角度，r代表的是一个另外的参数（RNP）,我是论文里要做两个图，一个是那个h随着角度的变化，一个是h随着RNP的变化，我做出来一个图，就是r在变的那个图做出来了，另一个做不出来

追答

我想我明白你的意思了，再给你看看，一会补充给你

追问

太谢谢啦，麻烦你了，我是要做三个图才要用matlab的，昨天刚刚开始学这个，学的不好，我在想是不是要用循环来做，可是我又不会编那个程序

追答

发现问题了，主要是：以theta为变量时，要不用弧度，要不用角度
用弧度，直接用sin()、cos()即可
用角度，则要用sind()、cosd()等
你的t=20:0.1:180，你看180弧度都多少个周期了
改为：t=pi/6:pi/50:pi;这样的就可以了