作业帮 > 综合 > 作业

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

来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/05/17 21:50:27
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文件里,路径要说清楚就行,求大神把程序贴上来或者发我私信,重赏啊重赏~
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来绘图的,数据少了会失真,图也会变得难看,你能不能再帮我看看?