求matlab绘图大神,隐函数绘图请问以下程序错在哪了?为啥第2个图作不出来?仔细讲解一下for循环那段,clc;clear all;w=[];z=[];figure;h=ezplot('w^(z/500)*0.01*(pi*sin(z/500*pi/2)+2*cos(z/500*pi/2)*log(w))+2*log(w)',[0,10^6

来源:学生作业帮助网 编辑:作业帮 时间:2024/10/06 16:45:21
求matlab绘图大神,隐函数绘图请问以下程序错在哪了?为啥第2个图作不出来?仔细讲解一下for循环那段,clc;clear all;w=[];z=[];figure;h=ezplot('w^(z/500)*0.01*(pi*sin(z/500*pi/2)+2*cos(z/500*pi/2)*log(w))+2*log(w)',[0,10^6
xNW_eEe{b~Auf/ekRU!P !$..eoѳ_SW9s33C|.I(M z{h5T69eʼn d\9mϽ1N.H1)FŤ L&Tl5>\֒w0-SA|~߫vt~M/fr2|őFI@`:4.d -,Сa/m$TJ.[ZH٦yo^y8jK~dLOrjY}Q£7WZ_ :\B g]fon/9h,Wͭ_SYѯ0 s ,Tw5WQY%g~crZME=x^2Qfx,~/P^o!0]=IlQb* 92:ż"A9oU({"~b%eQ$;U)rOvvWݼaݎQ:ɔcl7lKןw^M|7Uc{ќ{U7\?ܣp.ƅΉ*D.SWL\X Iى< SZW3^-56Ck1Mۖu]å)ɤcsN 9-&ln[418pcq7U\%f]u5Dom^} k`tb[x5Kyk Y|lya,l**[<7*+x^TvenQe| wXs,qjCbO0jz=[Gp횁IItH_`^gwdEf^rFӻM>ܩ1]!ǟzeP"(J$T-grl@̤rB`*$ 20() {RrL?E

求matlab绘图大神,隐函数绘图请问以下程序错在哪了?为啥第2个图作不出来?仔细讲解一下for循环那段,clc;clear all;w=[];z=[];figure;h=ezplot('w^(z/500)*0.01*(pi*sin(z/500*pi/2)+2*cos(z/500*pi/2)*log(w))+2*log(w)',[0,10^6
求matlab绘图大神,隐函数绘图
请问以下程序错在哪了?为啥第2个图作不出来?仔细讲解一下for循环那段,
clc;
clear all;
w=[];
z=[];
figure;
h=ezplot('w^(z/500)*0.01*(pi*sin(z/500*pi/2)+2*cos(z/500*pi/2)*log(w))+2*log(w)',[0,10^6,0,1000]);
x=get(h,'Xdata');
y=get(h,'Ydata');
[m,n]=size(x); %从这里开始到for循环结束没看懂,x指的是什么?从哪来的?
for i=1:m
for j=1:n
w=[w x{i,j}(1:end)];%这是什么?
z=[z y{i,j}(1:end)];
end
end
f=1000*(1+1./((w.^z*0.01).^2)+2*cos(z*pi/2)./(w.^z*0.01)).^0.5;
figure;
plot(w,f,'--r');
xlabel('w');
ylabel('f');

求matlab绘图大神,隐函数绘图请问以下程序错在哪了?为啥第2个图作不出来?仔细讲解一下for循环那段,clc;clear all;w=[];z=[];figure;h=ezplot('w^(z/500)*0.01*(pi*sin(z/500*pi/2)+2*cos(z/500*pi/2)*log(w))+2*log(w)',[0,10^6

改一下程序

clc;
clear all;
w=[];
z=[];
figure;
h=ezplot('w^(z/500)*0.01*(pi*sin(z/500*pi/2)+2*cos(z/500*pi/2)*log(w))+2*log(w)',[0.01,10,-10,10]);
x=get(h,'Xdata');% 获取上面画出的图所有点的横坐标
y=get(h,'Ydata'); %获取上面画出的图所有点的纵坐标
if strcmp(class(x),'double')&strcmp(class(y),'double') %判断是否是double类型如果是就转为cell类型
    x=num2cell(x);
    y=num2cell(y);
end
[m,n]=size(x); %求x大小,和y一样
for i=1:m
for j=1:n
w=[w x{i,j}(1:end)];)];%取出x的第i行j列的元胞,它是一个矩阵或者向量,不管怎么样,通过(1:end)转化为一个行向量放到w向量中.每次执行一次,w就增加一个行向量x{i,j}(1:end)

z=[z y{i,j}(1:end)];
end
end
f=1000*(1+1./((w.^z*0.01).^2)+2*cos(z*pi/2)./(w.^z*0.01)).^0.5;
figure;
plot(w,f,'--r');
xlabel('w');
ylabel('f');