Question

matlab编程高手进！我需要一个求方程在一定范围内的解的程序，方程不难

syms w y;
u=2.65e+10;
c1b=5.057e+10;
cv=2596.15;
cs=3128.89;
cp=3117.47;
b=11.6;
us01=0;
us12=20;
h=0.000000001;

f=w-sin(y*(b/h)*sqrt(cp^2/cv^2-1))*((u*(b/h)*sqrt(1-cp^2/cs^2)+(b/h)^2*us12)/(c1b*(b/h)*sqrt(cp^2/cv^2-1)))+cos(y*(b/h)*sqrt(cp^2/cv^2-1));
方程是f=0 ，他的图是类似cos函数的曲线，一个w对应一个y值，我想要y在区间【0，1nm】内w的值，y按0.01nm变化，1nm=0.000000001m，小数点后面有9位数，千万不能弄错。
程序算出来应该有100组数据，算出的数据可以显示在命令窗口，也可以存到txt文件里，路径要说清楚就行，求大神把程序贴上来或者发我私信，重赏啊重赏~~~

Answer 1

clear;
clc;
syms w y;
u=2.65e+10;
c1b=5.057e+10;
cv=2596.15;
cs=3128.89;
cp=3117.47;
b=11.6;
us01=0;
us12=20;
h=0.000000001;
h=ezplot(w-sin(y*(b/h)*sqrt(cp^2/cv^2-1))*((u*(b/h)*sqrt(1-cp^2/cs^2) ...
+(b/h)^2*us12)/(c1b*(b/h)*sqrt(cp^2/cv^2-1)))+cos(y*(b/h)*sqrt(cp^2/cv^2-1)),[-20 20 0 1e-9]);
title('f=0');
set(h,'LineStyle','--','color','r') ;
grid on;
view(90,-90);
C=get(h,'contourMatrix');
wData=C(1,2:end);
yData=C(2,2:end);

save ydata.txt yData -ascii
save wdata.txt wData -ascii

其实ezplot是无法设置步长的，只能用默认的
储存数据的路径就在当前matlab的工作路径

追问

可以不用ezplot来求数据？你求出来的数据太少啦，我之所以要带有步长的数据，是因为要带到origin来绘图的，数据少了会失真，图也会变得难看，你能不能再帮我看看？