Exame Modelo de AIPE

Olá a todos!

Já está on-line o exame modelo de AIPE. Juntamente com o exame coloquei na woc um vídeo de apoio, relativo à última questão de Labview.

Bom estudo,

Eu já estive a ver o exame modelo e o exercicio que me está a dar a volta á cabeça é o 1º...
se Ahem=2pirH e Aci=2pirh então o custo é (300+400)*2*pi*r*h ??
o que eu fiz foi:
syms r h;
h=solve('500000=(pi*h*r^2*(1+(2/3)*r))',h); %pr não passar a capacidade
for i=3:0.04:7
a=subs(h,r,i);
if(a<hmin)hmin=a;
end
if(a>hmax)hmax=a;
end
end

mensagem=['a altura minima é ' num2str(hmin) ' e a maxima é ' num2str(hmax) ' com o raio de 3 a 7'];
disp(mensagem);
----
Não é bem isto que se pede,mas pronto é uma experiência

Que tal uma abordagem diferente: arranjar uma função custo que tenha apenas uma variável independente (a altura e o raio não são independentes pois o depósito terá de ter o volume especificado). Depois trata-se de achar um mínimo dessa função! Como é que faziam isso em Análise Matemática? Achando o zero da primeira derivada....

Olá!
Reparei agora que no exame modelo houve um pequeno lapso na questão 1, nas fórmulas relativas à parte hemisférica do depósito. Já coloquei o update na woc!

Em relação á primeira pergunta: Eu já pus uma questão idêntica no forum próprio do MatLab mas visto que estão directamente relacionadas...Como representar C graficamente uma vez que é uma variavel simbólica??

É que me dá um erro estranho, logo não consigo resolve-lo

Já agora a minha visão em relação à resolução da questão:

syms Vt C R H Ac Ah pi % tudo vars simbolicas

R=[3 4 5 6 7] % R a variar entre 3 e 7

Vt=500000 % Volume total

Ac=2*pi*H*R % Area cilindro
Ah=2*pi*H*R % Area hemisferio
Vt=pi*H*R.^2+((2*pi*H*R.^3)/3) % Relaçao entre H e R atraves de volume

H=Vt/((pi*R.^2)+((2*pi*H*R.^3)/3)) % H em ordem a R

C=(300+400)*2*pi*H*R % C em ordem a H e R

plot(C,R,H) %grafico de C (ineficaz ???)


Qualquer ajuda é bem vinda

abraços

Faz o download da nova versão do exame modelo pois na versão inicial havia um pequeno erro nas fórmulas fornecidas na Questão 1.

syms h r 'positive';
R=3:0.1:7;%vou analisar 41 pontos
clc;
F1='500000=(pi*h*r^2)+(2/3)*pi*r^3';
F2='c=300*(2*pi*h*r)+400*(2*pi*r^2)';
x=1:41;
x=x';
y=1:41;
y=y';
n=0;
s=solve(F1,F2);
for i=3:0.1:7
n=n+1;
g1=subs(F1,'r',i);
g2=subs(F2,'r',i);
s=solve(g1,g2);

x(n,1)=double(s.c);
y(n,1)=double(s.h);
end

subplot(2,2,1);
plot(y,x);
xlabel('Altura');
ylabel('Custo');
subplot(2,2,2);
plot(R',x);
xlabel('Raio');
ylabel('Custo');
subplot(2,2,3);
plot3(R',y,x);
xlabel('Raio');
ylabel('Altura');
zlabel('custo');
subplot(2,2,4);
plot(R',y);
xlabel('Raio');
ylabel('Altura');

pause;close;

Eu ainda não sei bem se isto se faz assim...Ainda falta achar o minimo mas tal como o professor disse é uma questão de Diff...Com este codigo apenas analiso os graficos do custo

Vejam lá se descobrem o que isto faz:

syms H R real
V=2/3*pi*R^3+pi*R^2*H
H=solve(V-500,H)
C=2*pi*H*R*300+2*pi*R^2*400
f=diff(C)
solucao=solve(f)
for i=1:3
if conj(solucao(i))==solucao(i)
resultado=solucao(i);
break
else
continue
end
end
RAIO=double(resultado)

Na parte do 'for' penso que o codigo está mal por lapso do professor..mas a solução está bem pois a solução 1 é o que se pretende...
Obrigado professor por apresentar a resolução visto que a gente andava ás aranhas com este problema

'if conj(solucao(1))==solucao(1)' isto não havia de percorrer as 3 soluções?
Não haveira de ser assim 'if conj(solucao(i))==solucao(i)' ?

Sim, já corrigi! Lapsos destes nem o matlab os detecta!

Obrigado pela atenção (e participação no fórum!),

Espero que este não tenha muitos bugs...Se tiver agradeço que se manifestem

fich=fopen('c:\ascii.bin','\t')

A=fread(fich,inf,'char');
xlabel(A(1,1));
ylabel(A(2,1));
B=str2num(A(1,2:end));
C=str2num(A(2,2:end));
disp(B);
disp(C);
plot(B,C)


Grato por qualquer observação

a=1:10;b=2:11;
tudo1=['a';a'];
tudo2=['b';b'];
tudo=[tudo1 tudo2];
save -ascii 'c:\teste.txt' tudo
clear all;
load -ascii 'c:\teste.txt'
x=input('insira o nome da 1ª ','s');
y=input('insira o nome da 2ª ','s');
[m,n]=size(teste);

conta=0;
for i=1:n
if(strcmp(char(teste(1,i)),x)==1 & conta==0)
variavel1=char(teste(1,i));
a=teste(2:end,i);
conta=conta+1;
end
if(strcmp(char(teste(1,i)),y)==1 & conta==1)
variavel2=char(teste(1,i));
b=teste(2:end,i);
conta=conta+1;

end
end
if conta==0
disp('nada desse tipo encontrado');
clear all;

elseif conta==1
disp('apenas 1 lida');
clear all;
else
plot(a,b);
xlabel(variavel1);
ylabel(variavel2);
pause;
close;

end




A parte inicial foi para eu testar se dava

Markoni em [m n]=size(teste) podes por apenas m=size(teste), m fica com duas dimensões...não é importante é apenas optimização de código.

Já agora queria agradecer a todos a disponibilidade e a celeridade com que se vêm tratando estas questões ( é importante, pois o exame está próximo )

PS: sim, já vi que a minha solução está erradissima

Olá a todos!

Uma pequena ajuda nesta questão, que também está a dar que falar pelo que vejo aqui no fórum.

Vamos admitir que criámos um ficheiro no Excel onde colocámos os labels das colunas na primeira linha ('Sérgio' e 'Valor') e os dados numéricos nas linhas seguintes. Gravámos os dados em formato txt (Save as: Text(tab delimited)(*.txt)) em c:\teste.txt

Agora queremos abrir este ficheiro em matlab. Fiz o que vem a seguir:

>> var=importdata('c:\teste.txt')
var =
data: [2x2 double]
textdata: {'Sérgio' 'Valor'}
colheaders: {'Sérgio' 'Valor'}

>> var.data
ans =
1 4
2 5

>> var.textdata
ans =
'Sérgio' 'Valor'

Espero que isto ajude a iluminar as vossas mentes.

Bom estudo,

hum...com o importdata ele separa logo..pois mas no enunciado diz que está tudo em ascii logo é só numericos...mas sempre é util visto que carrega tudo para uma variavel só e depois é só manipular
Eu não gosto é do facto de ao ver as strings que carrega aparecerem com plicas

Não fui eu que fiz o enunciado dessa questão mas a ideia subjacente ao tipo de ficheiro é aquela que eu procurei transmitir na mensagem anterior. Este tipo de ficheiros é muitas vezes usado em sistemas de aquisição de dados.

Quanto às plicas, repara que se pedires ao utilizador uma string com o comando input e se ela ficar armazenada numa variável, as plicas também estão lá! É uma característica das strings

por acaso o professor não tem por ai mais exercicios desses(questão 1 e/ou 2 do exame modelo)?..para a gente estar mais á vontade com este tipo de exercicios

neste ultimo problema que o professor resolveu como é que eu posso aceder as colunas separadamente!
por exemplo termos a tabela

X Y
1 2
3 4

como é que eu posso aceder separamente ao vector X ou Y??

é simples..acedes ao indice
var.data(1:end,1)->o teu x
var.data(1:end,2)->o teu y

Boas!

Tentei resolver a questão 2 do exame modelo e fiz este pedaço de código. O problema é que não funciona e prtt decidi pô-lo aqui pra ouvir as vossas opiniões para perceber o que está mal e como corrigir.

Aqui está:

-admitindo que teste.txt é um ficheiro de dados em formato ASCII guardado na directoria c:\ (c:\teste.txt):

var=importdata('c:\teste.txt');

A=input('Introduza o número das duas colunas que deseja descarregar do ficheiro para adicionar ao gráfico no formato [num_col_1 num_col_2]: ');

col1=A(1);
col2=A(2);

vector1=var.data( :,col1);
vector2=var.data( :,col2);

plot(vector1,vector2);

vector1=vector1( : );
vector2=vector2( : );


xlabel(var.textdata(1));
ylabel(var.textdata(2));

fprintf('Dados descarregados do ficheiro de texto:\n%s\t%s\n%d\t%d ',var.textdata(1),var.textdata(2),vector1,vector2);

-------------------------------------------------------------------

Notas:

-Logo no inicio ao inicializar a variavel "var" dá 1 erro que é o seguinte: "??? Error using ==> importdata
Unable to open file."

-O fprintf não funciona correctamente. Já não contando que aparecem as plicas aparece também 1 warning:"The argument for the %s format specifier must be of type char (a string).".

Gostava de saber se la volta a dar a este pedaço de código...

"-Logo no inicio ao inicializar a variavel "var" dá 1 erro que é o seguinte: "??? Error using ==> importdata
Unable to open file."...
Isto acontece porque não deves ter o ficheiro criado;
vector1=vector1( : );
vector2=vector2( : );
isto não está a fazer nada...
e quanto ao fprint,este não apresenta variaveis do tipo 'campo',mas se fizeres um type cast...funciona,Já testei:o

obrigado.

vector1=vector1( : );
vector2=vector2( : );

Isto serve para que no fprintf a informação apareça em colunas pois vector1 e vector2 sao dois vectores linha, ou seja, transforma vector1 e vector2 em dois vectores coluna.

abraço

Tá bem visto...

No enunciado diz:"Sabe apenas que este ficheiro tem sempre um número par de vectores
(coluna) numéricos e que o primeiro elemento de cada vector tem uma string com o nome da variável."
Se está em ascii o ficheiro apenas tens numericos,ou seja, tens de descobrir qual a variavel que está na 1ªposição,ainda falta isto no teu codigo;além disso no enunciado diz que está tudo em colunas, por isso é que eu disse que aquela parte do codigo é totalmente inutil...
Eu no codigo que postei encima peço ao user as variaveis que ele quer carregar (penso que esse é o objectivo) e percorro as colunas do ficheiro fazendo um string compare...
Não sei se já resolves-te o problema que tinhas no fprintf mas aqui está o que eu faria:

fprintf('Dados descarregados do ficheiro de texto:\n%s\t%s\n',char(var.textdata(1)),char(var.textdata(2)));
fprintf('%d\t%d\n',vector1,vector2);


Espero ter ajudado