MatLab - рисование графика

[Cobert]

Здравствуйте. Мне потребовалось нарисовать 2 графика функции:
$ y=ax^b $ $ y=a^{bx} $ , где a и b - произвольные коэффициенты
Я стал действовать так:

Используется синтаксис Matlab M

x=-20:20
y=2*x^4
plot(x,y)

но программа выдала ошибку о том, что матрица должна быть квадратной и аргумент функции не должен содержать степень. Что мне нужно сделать, что бы избежать такой ошибки?

P.S. Я понимаю, что вопрос глупый, но я только вчера первый раз сел за MatLab, соответственно, ничего в нем не знаю.

Добавлено спустя 46 минут 12 секунд:

Все, тему можите закрывать. Я сам разобрался.

[photon]

для тех, кто вдруг столкнется с подобной проблемой, правильно записать:

y=2*x.^4; 

наличие точки обеспечит поэлементное возведение в степень, а не попытку умножить матрицу саму на себя, а точки с запятой - отключит вывод результата вычисления на экран (это медленная операция, в данном примере всего 51 точка - не страшно, но когда ворочаются массивы в десятки мегабайт, то их вывод на экран займет почти все время, которое требуется для выполнения программы )

[LLenich]

Здравствуйте!
Наверное это покажется простым, но тем не менее вопрос такой: скажите пожалуйста как построить график функции 3D, в случае, если данные заданы в файле в столбцах, каждая строка - определенная точка поверхности, какие команды я должен использовать, чтобы считать данные и построить поверхность?
Заранее спасибо тем, кто поможет.

[photon]

посмотрите help и примеры.

Для считывания из файла функция load(), для построения поверхности mesh(), surf(), contour()

[ann-mari]

$$

помогите,пожалуйста.Как на графике отметить корни и экстремумы?

[photon]

Вычислите положение нулей и экстремумов, соберите их в отдельный массив и постройте на том же графике какими-нибудь выделяющимися символами.

[ewert]

Типичный код:

Используется синтаксис Matlab M

t=a:(b-a)/300:b;
z=[];
for k=1:length(t)-1
    if t(i)*t(i+1)<=0
        z=[z, fzero('f', t(i))];       % очередной корень
    end       
end

plot(t,f(t), z,z*0,'o'),    grid

(Некоторые корни могут случайно сдублироваться, но для графика это неважно.)

[photon]

ewert, лучше избегать циклов в матлабе, а работать с массивами. Я бы использовал конструкцию типа

Используется синтаксис Matlab M

t(1:end-1).*t(2:end)

[ewert]

photon в сообщении #256630 писал(а):

ewert, лучше избегать циклов в матлабе, а работать с массивами.

Я в курсе; но что Вы собираетесь делать с этим замечательным массивом дальше? Всё равно ведь потом придётся его прореживать.

Цикл дешевле обойдётся. Другое дело, что (в целях оптимизации) лучше применить его к примерно Вашему массиву. Точнее, к примерно такому:

Используется синтаксис Matlab M

y=f(t);
find=(y([1:length(y)-1])*y([2:length(y)]<=0);
for i=1:length(find)
    if find(i)
    . . . . . . . . .

[photon]

ewert, зачем же так циклить? А уж применять в качестве имен переменных названия стандартных функций совсем нехорошо.

Используется синтаксис Matlab M

x=1:0.01:20;
y=sin(x);
plot(x,y)
hold on
roots_pos=find(sign(y(1:end-1).*y(2:end))-1);
plot(x(roots_pos),y(roots_pos),'o', 'Color', 'red'), grid on

ну можно дополнить использованием FZero, если хотите

[ewert]

Я просто никогда с командой find не имел дела.

Но -- тоже выражу некоторое недоумение.

Во-первых: зачем же холдить, если можно вывести график сразу?

Во-вторых: гораздо проще и понятнее

Используется синтаксис Matlab M

roots_pos=find(y(1:end-1).*y(2:end)<=0)

[photon]

ewert, насчет упрощения с find() - согласен, насчет выводить вместе - наглядность теряется, а так, конечно, можно. Почему взял sign? - по идее это очень быстрая функция

[ewert]

Ну, операция сравнения (по идее) ещё эффективнее. Хотя бы потому, что она одна, а не две.

А hold я не люблю без необходимости. Вечно приходится не забывать его отключить.

[wtf]

Всех приветствую. Нужна помощь в посторении графика в Матлабе с помощью plot.

F=3x1+x2 -> max
при ограничениях
х1 + 3х2 <= 21
2х1+ х2 <= 12
х1 <= 5
x1>=0, x2 >=0
НУжен тест построения графика в матлабе. Заранее спасибо.

[ewert]

Графика чего?...